Stability in the Stefan problem with surface tension (II) 



Mahir Hadzic and Yan Guo 
Division of Applied Mathematics 
Brown University, Providence, RI 02912, U.S.A. 

May 22, 2009 

Abstract 

Continuing our study of the Stefan problem with surface tension effect, in this paper, we estabhsh 
sharp nonhnear stabihty and instabihty of steady circles. Our nonlinear stability proof relies on an 
energy method along the moving domain, and the discovery of a new 'momentum conservation law'. Our 
nonlinear instability proof relies on a variational framework which leads to the sharp growth rate estimate 
for the linearized problem, as well as a bootstrap framework to overcome the nonlinear perturbation with 
severe high-order derivatives. 
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1 Introduction 

The Stefan problem is one of the best known parabolic two-phase free boundary problems. It is a simple 
model of phase transitions in liquid-solid systems. 

Let 51 C denote a domain that contains a liquid and a solid separated by an interface F. As the melting 
or cooling take place the boundary moves and we are naturally led to a free boundary problem. Define the 
solid phase il.~{t) as a region encircled by r{t) and define the liquid phase r2+(t) The unknowns 

are the location of the interface |r(i);t>o| and the temperature function v: [0,T] x fi^M. Let Fq be the 

initial position of the free boundary and tio : ^ M be the initial temperature. We denote the normal velocity 
of F by and normalize it to be positive if F is locally expanding (t) . Furthermore, we denote the mean 
curvature of F by k. With these notations, (w,F) satisfies the following free boundary value problem: 

dtv-Av^O in n\T. (1.1) 

v = aK. on F, (1.2) 

F-K]t on F, (1.3) 

Vn^O on dn, (1.4) 

v{0,-) = vo; F(0)=Fo (1.5) 

where i7>0 is the surface tension coefficient. Given v, we write and for the restriction of v to 57^ (t) 
and 0~(t), respectively. With this notation [vn]lL stands for the jump of the normal derivatives across the 
interface T{t), namely [vn]- •=^n j where n stands for the unit normal on the hypcrsurfacc T{t) with 
respect to il~{t). 

Steady states and the perturbation {u,f). It is well known (cf. [55]) that there are steady circles 
for the problem (|l.ip - (|1.5p and they are parametrized by triples {R,xo,yo), where ^>0 and xo,yo are real 
numbers such that the ball Bj^(xo,yo) with radius R centered at (xo^yo) belongs to Q. To each such triple, 
we associate a time-independent solution given by (v,r) = (a/ R,Sji{xo,yo)), where Sji{xo,yo) = dBj^{xo,yo). 

Coordinates. We shall assume that $7 = 5/}, (0) CM^ is a ball of radius i?* > 1. We shall parametrize F 
as a 'graph' over the unit circle centered at (a;o,?/o) with ^/x^ + y^ + 1 < i?*: 

r{t) = {xen\x= {R{t, 6) cose + xn, R{t, 9) sind + yo) } , (1-6) 

where Here i?: [0,oo[xS^ is a sufficiently smooth function such that Uo<t<T-'^(^) T>0 

and additionally Fo = {a;efi| a;= (_Ro(^)cos0-|-xo,i?o(6')sin^?-|-j/o)} for some i?o>0. We introduce a map 
0: [0,c»[x§i^ [0,oo[xF, 

0(t, 9) (i?(i, 9) cose* + a;o, 6i) sin 61 + 2/0) ■ (1-7) 

The map induces a metric on F, whose line element is denoted by \g\d9 and in local coordinates \g\ = 
R^ + Rg. With these notations, we can express the unit tangent r, the unit normal n, the mean curvature 
K and the normal velocity V in local coordinates, i.e. as a function of 6. It follows that 

(Re cos9 — Rsm9, Rg sinfl + i?cos6') (Re sm9 + Rcosd , —Re cos9 + Rsm9) 

; — ; , n — 



1-91 (1.8) 
,o0 = if=_--y)^ and Vo^^^,n^^^. 

In the rest of the paper we shall normalize a=l and analyze the stability of the steady state (w,r) = 
{l,Si{xo,yo))- Hence, it is natural to introduce the perturbation / : [0,oo[x§^ from the steady circle 
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Si{xo,yo) by setting f:=R—l. Similarly we set u:=v — l. Linearizing the curvature H and the Jacobian 
\g\ around 1, we obtain 

Hif) = l-f-fge + N{f); |g| = 1 + ^ + (1.9) 

where N{f) denotes the nonlinear remainder in the expansion of H{f) and = 0(|/|^ + j/el^). If we set 
^^0 = 1^0"! and /o = -Ro~lj we can formulate the problem in terms of the perturbation {u,f): 

dtu-Au = in n\r. (1.10) 

on r, (1.11) 

y = on r, (1.12) 

u„ = on dn, (1.13) 

u{0,-)^uo; /(0) = /o (1.14) 

Motivation. To motivate our results let us first investigate the eigenvalue problem for the linearized 
Stefan problem around a steady circle (centered at (0,0)): 

\v~Av^O onn\E,\ (1.15) 

v^-f-fee on§\ (1.16) 

[vn]t = -Xf onS\ (1.17) 

Vn = ondn; [v]t = on §\ (1.18) 

with A^O. Observe that RHS of the equation (|1.16p simply is the linearization around 1 of the curvature 
operator H. Note that by integrating ()1.15|) over ri\S^ and using the boundary condition ()1.17p . we obtain 

the first important identity In'^'^ Is^ Multiplying (jl.isp by v, integrating over il\S^ and using 

the boundary conditions (|1.16p and (jl.lTp . we are led to the identity 

{fl-n)+l |VH^-0. (1.19) 

Setting J('i;,/) := j^^v"^ ^ k^{fe-P] ^nd denoting Pw = t; - |^ J^u, P/ /- |^ Jgi /, we note that 

where we note that X{v,f)= J^\Pv\'^ + /gi {/| — (P/)^} > by Wirtinger's inequality and JfiV — — Jsif- 
Hence, if 

1 1 

then we easily see that for A > 0, equation (|1.19p can not have non-trivial solutions and thus, all the eigenvalues 
are non-positive. On the other hand, if C < 0, one can show that there exists a pair (v,f) such that X(w, /) < 0. 
One can formulate a variational problem so that there exists Aq > 



^0 '^.oVw^ J 



where the constraint set C is given by 

C :={(«,/) I V e H\nsi), f e H'^/^iS'), V = ~f-fee on S\ fe'" = 0, ^ \Wv\^ ^ o}. 

The minimizer {vo,zq) is the eigenvector associated to Aq and leads to the fastest growing mode e^^{vQ,zo) 
for the linear Stefan problem. For a complete analysis of the eigenvalue problem via a different method, 
see EH. 
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The above analysis suggests that in the stable case C > 0, we should expect the solutions to exist globally 
close to the steady sphere Contrary to that, in the unstable case C<0, an instability should develop. 

Notation. The normal velocity of a moving surface E is denoted by Vs and the index S is dropped 
when it is clear what surface we are referring to. We also set ils :=f7\S]. For a given function f E L'^{E>^), 
we introduce the Fourier decomposition 

oo 

f{9) — Qfc cos k6 + bk sin k9. 

k=0 

For any n>0, we define the projection operators P„ and P„+: 

oo 

P„/:=a„cosn0 + fe„sinn6', P„+/ := afccosfc0 + 6fesinfc6'. (1-20) 

k—n 

For the sake of simplicity, we denote P/:=Pi+/ = / — Pq/. If a function w is prescribed on fl then 'Pqw:— 
JqW and Pwi^w — PqW. Furthermore, we will typically leave out the domain when referring to a 

Sobolev norm of functions g:S^— >R, i.e. H^H/f := ||5||ir=(si)- We shall often employ the following notation: 
d^g :—de<^dtbg. For a given set AcM.'^, we denote the closure of A by cl(A). For any te]0,oo] introduce 

C°°{[0,t[- r!±) [u:[Q,t] x^^R; ulci(o±) e C°°([0,i[xcl(r!±))}, 

and analogously define C°°{n^) by suppressing the time dependence. For any fcGNU{0} set 

i/''(f}=^) = |u:f7^M; ulci(o±) e-H"''(cl(f^^))|; ||u||ff'=(o±) := ||wlci(n+)||//'=(ci(n+)) + ||u1c1(o-)||h'^(c1(o-)), 

(1.21) 

where H'^ stands for the standard W'^''^ Sobolev space. 

Conservation laws. Assume for the moment u and / are classical solutions to (ll.lOp - (I1.12p . The 
first 'conservation-of-mass' law arises by simply integrating p.lOp over Qy, using the Stokes' formula and 
the boundary conditions: 

dt f u + dt f {f + ^}^0. (1.22) 

The second law states the fundamental energy identity; it arises by multiplying p.ip by (1 + u), integrating 
over n and making use of the fact that Jj,V k dS = dt\T{t)\ (where |r(t)| = Jgi \g\d9 stands for the surface 
volume of r(t)): 

f {l + u)' + dt f J{l + f ? + pB+ I |Vup = 0. (1.23) 

In addition to these two well-known identities, we discover a new 'conservation-of-momentum' law. It 

arises by multiplying the equation (11.101) by the harmonic functions pa(x,y) — x+ -prf and ph(x,y) = y + -prr 

I I I y I 

respectively. After integrating over fJpj using the Green's identity and the boundary conditions, we arrive 
at: 

(1.24) 



dt f upa = dt f {Fa{RS-Fa{l,e)}; dt f upb = dt f {Fb(i?,0)-Ffc(l,0)}. 
Jn Jsi Jn Jsi 



Here 



R 



2 rcos^ + XQ 



Fa{R,9):= {rcose + xo + Ri-, . ., )dr (1.25) 

^ (rcostZ + Xo)^ + (rsmt^ + yo) 

FtiR,9):= [""{rsme + yo + Rl-. J^^i^^+^o ^^^^^^ 

' ^ (rcos0-|-a;o)^-|-(rsmt'-|-2/o) 



The details of the proof are presented in Lemma 13^ Let us denote 



™o= /_"o+ /_ {/o + ^}, (1.27) 



n 



noPa- / {Fa{R{O),0)~Fa{l,0)}; mb:= / uoPb~ / {Fb{R{0),e) - F,{l,e)} . (1.28) 
n Jn Jsi 

These conservation laws play an important role in our proof of the stability. 



4 



1.1 Main theorems in the stable case 

Motivated by (|1.23p . we naturally define the zero-tli order energy and dissipation as 



{uj):^l f fu+f + + I?(o)(Zi,c^)= / (1.29) 

^ Jn Jn jsi * Jn 



Let / S N and I > 3. By we shall denote the energy quantity involving ^-derivatives 



I „ /-I 



£:(+)(z^,^) = i^ / uf, + lYl I \^^tA^ + \Yl I (1-30) 

In the same vain, is defined through: 

/ Uf.+Y, I |VW,.|2 + ^ / Ll^^^l]. (1.31) 

Of course, from these definitions it is neither clear that 5(q) defines a positive definite functional, nor it is 
clear that and 'D^j^-^ are positive definite, due to the presence of the term —J^-^ui'^f, in the definitions 
above. However, under the appropriate smallness assumptions the positive definiteness will become apparent 
later. Summing (|1.29p . (|1.30p and (|1.3ip . we introduce the temporal energy and the temporal dissipation: 

£{U,uj):=£i^o){^,uj)+£^+){U,uj), (1.32) 

V{U,uj):^V^n){U,uj)+Vi^+^{U,uj). (1.33) 

In addition to the (temporal) energy £ and the (temporal) dissipation V, we introduce the higher order 
space-time energy quantities £ and T) that incorporate mixed space-time derivatives. For any 0<v<l, we 
set: 

i-i i-i 
<&^{uJ):=Y.^'>^,{uJ)- D''(u,/):=^;.«D,(u,/), (1.34) 

g=0 (3=0 

where, for Q<q<l~2 



^ JS^ \9\ 

(1.35) 

where P2+ is given in (|1.20p . and 

£z_i(7.,/):=i||Vii||^.„_„(^±). (1.36) 

Observe that 2;^ is positive definite under the smallness assumption on /. Namely, for any l<q<l — 2 set 
p = dfl\^qf ■ Then the integral over §^ in the definition of £g takes the form 

The first inequality above simply follows from Wirtinger's inequality, and 1 1 \ g \ lq+'i ~ 1| < C*! |/| is 
small. Similarly, for any < q < ^ ^ 1 

2),(u,/):=||Vwt,-i-,|p^.,+i(,,±) +229-1^^ {|/0.,+i,,-,|2-|/,„,,_,|2|>o, (1.37) 

by the Wirtinger's inequality. When dealing with the case i^=l, we denote the quantities (B^{u,f) and 
D^{u,f) by t{u,f) and D{u,f) respectively. Finally, for any 0<q;, 0<t^<l we define the total energy and 
total dissipation as: 
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Theorem 1.1 Assume C>0. Then there exist 0<a,i^<l, a constant C* and a sufficiently small constant 
M* > such that if 

Ea .U 

iuQjo) + C^iml + ml + ml)< M* , 
then there exists a unique solution to the Stefan problem il.lO]) - ^Ll^ {u,f), satisfying the global bound 

1 /■* 

Eo,,,{u,f){t) + C^{ml + ml + ml) + - / D^^Au,f){T)dT<Eo,,,{uJ){s) + C,{ml + ml + ml), t>s>0. 

(1.38) 

The second theorem improves the stabihty to an asymptotic stabihty statement, and states that the solution 
converges to a 'nearby' steady sphere: 

Theorem 1.2 Under the assumptions from the previous theorem, there exists (a;o,yo) close to (a;o,?/o) o,nd 
R close to 1 such that the solution {v,T) to the Stefan problem converges exponentially fast to the steady 
state {l/R,Sji{xo,yo)). 

Remark. The analogous statement holds for the steady states {v,T) = {a/ R,Sji{xo,yo)) of the original Stefan 
problem (jl.ip - ()1.3|) . Here, the stability condition reads (^{a,R) :— — \g}\fi2 > 0, where IS*/;.! = 2TrR is the 

length of the circle of radius R. 

Remark. The role of the weight v in the definition (|1.34p of <B'^ and is technical in nature. It is 
introduced naturally in order to close the energy estimates in Subsection 13.2.51 

1.2 Main theorem in the unstable case 

As it is shown in Subsection [13 with the help of an appropriate change of variables, the Stefan problem can 
be formulated as a problem on a fixed domain flsj^^xo.yo) ^'^^ ^^^^ unknowns w:flsi{xo.vo) ^md /:S^— >M 
and it takes the form 

dt{wJ)^C{wJ) + U{wJ) (1.39) 

where the linear operator C and the nonlinearity U are defined explicitly by (|4.400p and (|4.424p respectively. 
In order to formulate the instability statement, we need to specify the notion of the initial perturbation. 
Lemma 231 based on [35], under the instability assumption C<0, shows that the linearized operator for the 
Stefan problem possesses countably many real eigenvalues Ao,Ai,A2,... with finite multiplicity. The leading 
eigenvalue Aq is positive and simple, Ai = and all the remaining eigenvalues are strictly negative. If m{i) 
stands for the multiplicity of the eigenvalue Ai, we associate the eigenvectors Ci^k to the eigenvalue A^, where 
l<fc<TO(z). The set U°io ^^^=1 forms an orthonormal basis with respect to the inner product (•,•)/ 

defined by (|4.409p . We refer to 

oo ni(i) 

{wo,fo) = ^^Ci^ket,k; co,iy^O (1.40) 

i=0 fc=l 

as the generic profile if there exists a non-trivial contribution along the direction of the growing eigenvector 
(vqjZq), i.e. co,i =/=0. In the following theorem we assert the local existence for small initial data and nonlinear 
instability of the steady spheres. 

Theorem 1.3 For any smooth generic profile such that ||wo||^2(-q-) + ||/o|li2 = 1, there exists a sufficiently 
small number 6o>0, such that for any S>0, there exists a unique solution {w^,f^) to U.39\) with initial 
conditions {iSwo,Sfo)} ^ which on the time interval [0,T^ :— ^\og^] satisfies 

\\{w'j'){t)~c„,iSe^°\vo,zo)\\L-<C6^e^^"' + CS, 

where C is independent of S. Furthermore, there exists a constant Ko>0 independent of S, such that 

\\^'iT')\\Un) + \\fHT')\\h>Ko. (1.41) 

Remark. Our theorem shows that the nonlinear unstable dynamics is characterized by the growing mode 
e^°'^{vQ,zo) over the time scale [0,T'^]. If we perturb away from 5i(0,0)=§"'^ we show that vq is spherically 
symmetric ((wo,zo) is the eigenvector associated to Aq), hence our instability does NOT contain morphological 
changes, but only spreading or shrinking along the radial direction. 
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1.3 Previous work 



The Stefan problem has been studied in a variety of mathematical literature over the past century (see for 
instance [31]). If we set cr = in (|1.2p . the resulting problem is called the classical Stefan problem. It has 
been known that the classical Stefan problem admits unique global weak solutions in several dimensions 
([13], [13] and [13]). The references to the regularity of weak solutions of the two-phase classical Stefan 
problem are, among others, f2^, f3j, [4], [9]. Local classical solutions are established in [22] and [27] . 

If the diffusion equation (II. ip is replaced by the elliptic equation Au = 0, then the resulting problem 
is called the quasi-stationary Stefan problem with surface tension (also known as Hele-Shaw or MuUins- 
Sekerka problem). Global existence for the two-phase quasi-stationary Stefan problem close to a sphere in 
two dimensions has been obtained in [5] and , and in arbitrary dimensions in [12| . Global stability for the 
one-phase quasi-stationary Stefan problem is established in [I6] . Local-in-time solutions in parabolic Holder 
spaces in arbitrary dimensions are established in f6]. 

As to the Stefan problem with surface tension (also known as the Stefan problem with Gibbs-Thomson 
correction), global weak existence theory (without uniqueness) is developed in [T], [5S] and [30] ■ In [H] the 
authors consider the Stefan problem with small surface tension i.e. cr^ 1 in (|1.2p . The local existence for the 
Stefan problem is studied in [55] . In [TU] the authors prove a local existence and uniqueness result in suitable 
Besov spaces, relying on the L^-regularity theory. Linear stability and instability results for spheres (which 
are equilibria for the Stefan problem with surface tension) are contained in |28j . Some of the references for 
the Stefan problem with surface tension and kinetic undercooling effects are [7], [29], [31], [33] . 

1.4 Strategy of the proofs and basic ideas 

In [U], the authors constructed the first global smooth solution near a flat surface. This article is a continu- 
ation of such a study. However, there are several major mathematical novelties in our current work. Let us 
first note that just from the definition (|1.29p of £(o) it is not even clear that this quantity is positive definite. 
An important step is to Taylor-expand the term \/{i + f)^ + fg— 1. We obtain 

^(0) = ^ / /^+ / f+l I fo+ I 0{\f\^ + \fe?) 

To extract a positive quadratic contribution, we exploit the conservation-of-mass law (|1.22p and note that 
JqU + /gi / = "j/gi Using this observation, by a calculation similar to the one in the equation follow- 
ing (|1.19p . we can conclude 

£(0)=^ / ip«p+k( / uY+l f {/i-ip/n+o(i/p+i/.n 

Note that Jg^ {/e >0 by the Wirtinger's inequality. Hence, if C>0 it does seem that £(o) is non- 

negative, modulo the third order remainder. Unfortunately, note that the L^-norm of the first mode Pi/ 
(recall (fOO]) ) is certainly not controlled by the expression /gi [fj - \ff\^] = /gi {IP2+/0P - IP2+/P}. The 
analysis of the first modes turns out to be impossible without the new conservation-of- momentum law (|1.24p . 
More specifically, the problem lies in ensuring the positivity of £(0) due to the non-trivial presence of the 
first modes in the third order remainder. The key to circumvent this severe difficulty is to exploit (|1.24p . so 
to write: 

{l + Rl)ai^ f upa+ f 0(|/|2 + |/2|), {l + Rl)b,= f uph+ f 0{\f\' + \f^\); 
Jn Jsi Jo Jsi 

and thus, we can bound the L^-norm of Pi/ by the L^-norms of u and (Pp-f P2+)/, which are in turn 
bounded by the energy f (0) . See Lemma 13.101 

Energy. Unlike our previous work [21], the higher-order energy estimate is considerably more delicate 
due to the geometric structure of our problem. In order to preserve the boundary condition, one needs to 
take tangential derivatives. Unfortunately, usual spatial derivatives d^i could lead to linear growth in time 
when we interchange from one local chart to another, which presents a difficulty in constructing the global 
solution. Instead, we divide our main energy estimate in two steps. We note that dt is the only globally 
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defined vector field wfiicfi is 'almost' tangential with respect to r{t), so we start by deriving the estimate for 
the temporal energy, defined as in (|1.32[) and (|1.33[) : 



sup £{s)+ [ V{s)<£{Q) + C sup (Vf + VS) / V{s)ds. (1.42) 

0<s<t Jo 0<s<t Jo 

Instead of deriving the similar energy estimates for the other tangential derivatives, we make use of the 
temporal energy estimate (|1.42p and the heat equation in the bulk to bound £ by f and V. We refer to this 
as the reduction step. The estimate we obtain for the space-time energy is the following: 

sup (£''(s)+ / D''(s)<(£''(0) + C sup Ve' [ 'D''{s)ds + C* [ V{s)ds. (1.43) 

0<s<t Jo 0<s<t Jo Jo 

Based on the previous two estimates, we are able to obtain the a-priori bound for the solution of the original 
Stefan problem 

1 /■* 

sup Ea,^{s) + - / D^_^{s)ds<Ea_^{0), 

0<s<t ^ Jo 

which is the central estimate that will allow us to extend the solution globally. We chose to define our energy 
on a moving domain, with the aim of preserving the transparency of the geometric structure of the problem. 
We prove the temporal energy estimates in Subsection 13.2.31 and with the space-time energy estimates in 
Subsection 13.2.51 

Regular izat ion and the iterative procedure. In order to establish the local existence for the Stefan 
problem, we set-up an iteration scheme, which generates a sequence of iterates {(m™, /"*)}„£«■ As in [2T| . 
such an iteration is well defined, but it breaks the natural energy setting due to the lack of exact cancellations 
in the presence of cross-terms. We design the elliptic regularization 

ftR feH r 1+ , 

m m 

to overcome this difficulty. For a fixed e, we use it to prove that {(w'",/'")}meN is a Cauchy sequence in 
the energy space. However, note that the expression m™"*"^ — does not a-priori make sense, since the 
iterates u™ and u'"^^ are defined on different domains 17™"^ and fi™, respectively. Thus, we transform the 
regularized Stefan problem to a problem on a fixed domain by applying a suitable change of variables and 
then prove the Cauchy property. 

Stability for C > 0. After establishing the existence of the global solution under the positivity assumption 
on C, we further prove the exponential decay of the solution to a nearby steady circle which is uniquely 
determined via conservation laws (|1.22p and (|1.24p . Just like in pP, we have to prove that the energy E^^u 
is bounded by the dissipation Dq. j^. However, the situation is considerably more complicated in comparison 
to the Stefan problem close to flat steady surfaces. The reason is that the Dirichlet condition (|1.2p can 
not be directly used to estimate the L^-norm of / against the L^-norm of Vu, since Jp(K— 1)7^0. Again, 
we crucially exploit the conservation laws (|1.22p and (|1.24p together with the condition C>0, to prove this 
estimate. This is carried out in Subsection 13.41 

Instability for C < 0. The passage from the linear to nonlinear instability is a delicate problem in 
general. In Section 3] we start by proving the a-priori bounds and the local existence under the assumption 
C<0 (Lemma 14.11 and Theorem 14. 2p . analogously to the stable case. However, to prove the instability, we 
resort to a bootstrap technique as developed by Guo and Strauss in |19]. The main idea is to show that on 
the time scale of instability, some stronger norm 1 1 1 - 1 1 1 (given by the high-order energy) is actually controlled 
by the weaker norm || • || (which is simply i^-norm in our case). The details of this general procedure are 
described in Lemma 14.101 However, in order to apply Lemma 14.101 a detailed understanding of the linear 
problem and the growing eigenmode is necessary. In Subsection l4.31 we provide a variational characterization 
of the positive eigenvalue (Lemma l4.5p and establish sharp growth estimate for the linear problem in suitable 
norms (Lemma 14. 9p . To finally prove Theorem 11.31 (Subsection I4.4p , we have to reformulate the Stefan 
problem so that the Duhamel principle can be applied. Again, wc transform the problem to a problem 
on the fixed domain with the severe nonlinearity N{f) coming from the curvature term, contained in the 
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Dirichlet boundary condition. To ensure the linearity of the boundary condition, we further modify the 
unknowns and as a resuh obtain the problem (|4.420p - (|4.423p . to which we apply Lemma [4. 101 

This work is a continuation of the program of developing a robust energy method to investigate and 
characterize morphological stabilities/instabilities in the setting of free boundary problems in PDE. In par- 
ticular, in a forthcoming paper |2pi, the nonlinear stability and instability of steady spheres in arbitrary 
dimensions will be analyzed. 

Due to the rotational invariance, in the rest of the paper we shall, without loss of generality, analyze the 
stability of the steady states with yo = 0, i.e. of the form (?},f ) = (l,5'i(xo,0)). 

2 Energy identities 
2.1 Regularization 

In order to prove Theorem II. 1[ we shall regularize the Stefan problem by adding viscosity to the 'jump' 
relation (jl.3p . Namely, for e>0, we consider the following equation: 



V + eA=[un]t onT, (2.44) 
where Ao(j)= Thus, expressed as an equation on the equation (|2.44p reads 

ftR fen r 1+ , 

where (f> is given by (|1.7p {yo — Q). We shall refer to the problem (|l.ip . (II. 2p and (|2.44p as the regularized 
Stefan problem. In the rest of this subsection we define a number of energy quantities that will be necessary 
for our subsequent energy estimates. 

The stable case. We introduce the additional e-dependent energy terms that incorporate the viscosity: 

S,{U,u:):=8{U,u:)+S\lo)- f^M:=-^e/ (2-45) 

i-i . 

V,{U,u;):^V{U,Lo)+V'{uj); := Ve / -cc-^.^.+i}. (2.46) 

k=o -^s^ 

Accordingly, with regard to the higher-order space-time norms € and S), for any < g < / — 1 we define: 

€l{f) := e22«^^ {|/,.,+4,,-i-, p - |/,.,^3,,-i-, p} - £2^* |P2+/,2<,+4,,-i-, p ( ^^^^^^^3 - l) , (2.47) 



||/e29+3ti-g 1^ — |/e2,+2ti-<j (2.48) 
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We define the total higher order energies in the natural way: 

i-i i-i 
r^'^^^i/'^e^; €^ = £"4- r^''; S)"'"^ = ^i^«2)^; = S"^ -f 2)^^^ (2.49) 

q=0 q=0 

Again, if iy=l, we simply drop the index 1 in the notation introduced in (|2.49p . 

The unstable case. In order to prove the local existence theorem in the case C < 0, we introduce the 
appropriate energy quantities. Using ()1.20p . we define the energy A4 by 

M{ll,u;):=^Y. [ Uf. + lY. I |VW*^I' + |E / " 1^2+^.,. p}; (2.50) 
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and the dissipation N is defined through: 



^{U,u;):^Y. / ^''=+E / |VZ^t.p + E / HV-|P2+^,.p}. (2.51) 



The quantities and Af are the analogues of £ and I? defined by p.32p and (|1.33p . In analogy to (|1.34p . 
for any < 1 we define the space-time energy quantities: 

/-I i-i 
Tl" = ^i^''Tlg, OT' = ^z^««n, (2.52) 

g=0 q=0 

where, for any 0<q<l — l, we define 

OT,(w,/):=||Vu,,-l-,|||,.,(^)+229 f {|/e.,+2t,-l_,|2-|P2+/e.,+i,,_i_,|2}; 

(2.53) 

ai,(?/,/):=||Vut,-i-,||^.,+i(^)+22'J / {|/e.,+i,,-,|2-|P2+/,,,,,_,|2}. 



§1 



Finally, for any 0</3, 0<i^< 1, the toial energy and the toia^ dissipation for the unstable case are defined 

by 

As to the energy quantities that arise from the regularization (|2.44p . analogously to (|2.45p - (|2.49p . we define 

/-I » 

M,iU,cu):=M{U,iu) + M'{u;); M'{uj):=Y,^ {uj^gst, -\V2+^02t''\^} (2.54) 

fe=o "'s' 

and 

Afe{ll,Lj):^J\f{U,uj)+Af'{uj);Af'{uj):^S2e / {c^^^^.+i - |P2+c<Je2t<e+i P}. (2.55) 

Furthermore, for any 0<q<l — l we define: 

fm^(/):=e229 / {|/e„+4,,-i-,|2-|P2+/e.,+3,,-i-,|2};^^(/):=e22« / {|/,.,+3t,-,|2- |P2+ W.^^-^j^}, 

(2.56) 

and hence the total higher order energies are given by 

;-i i-i 

= an^ = aH"+fm^'"; ^^-^^^^z/'fri^; o^.^^DT'+fri^^". (2.57) 

13=0 q=0 

2.2 Trace 

Let 7 stand for the trace operator on the curve F. By 9^* we shall denote the time differentiation operator 
along the hypersurface F defined by 

a;ito0 = 9t(uo0). (2.58) 
With these notations, applying the chain rule in the above formula, we arrive at the following decomposition: 

d;u^-f[ut] + Vun + ViiUs, (2.59) 

where o (/) = (cos 6*, sin 6) ■t = is the tangential velocity, V stands for the normal velocity of the curve 
C and Us — dgU is the tangential derivative of u (recall (|1.8p ). Note that the pull-back d* of the operator ds, 
with respect to the map </) : ^ F, is given by i9* — j^^dg (in short, in local coordinates ds = j^c^e)- Iterating 
the formula (|2.59p k times and using the product rule, we arrive at the following decomposition formula for 
the operator dfk : 

j[utk]=d*kU-['-/,dtk]{u), (2.60) 
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where the commutator [7,9(fc] (u) is defined through 

k-l 

[jA'']{u):^Y1 E CP'''''d:.-.-AV'V[)d:d?,ut., (2.61) 

p— q+r—k—p 

and C^-i-^ are some positive real constants. The following trace lemma will be of crucial importance for 
estimating terms appearing in (|2.6ip : 

Lemma 2.1 Let vEH'^^^{il), /eiJ*''(S^) and q,r£N such that r + q = k. Then the following inequality 
holds: 

\\^:^?,v\\LH^)<C{^ + \\f\\H^)\\v\\H.^^^n)■ 
Proof. Using the formulas dn — Uidxi and ds —Tidx^^ we may write 

r 

Here Fj (n) is a function involving at most j — 1 9* -derivatives of n and thus at most j c^e-derivatives of /, since 
n = l + 0{fg). On the other hand Qm{u) consists of terms containing m Cartesian derivatives of m. The idea is 
to use the Sobolev estimates: if j < we estimate < C'll-Pjllffi by the Sobolev inequality, whereas 

Qu+j-i is estimated in L^-norm. We then use the trace inequality, to conclude ||(3fe4.j_i||i2 < C||u||//fc+j(Q). 
In the case j >-^^, we estimate Fj in i^-norm and Qk+j-i in L°°-norm together with the Sobolev and 
trace inequalities. Combining the two cases the claim of the lemma easily follows. Note that we only use the 
trace inequality for the Cartesian derivatives on T. Observe that in the inequality ||f ||l2(p) < C||i'||fl'i(n), 
the constant C depends only on the C^-norm of /. However, for fc> 2, the C^-norm of / is bounded by the 
H^+^-noYTn of /. □ 



2.3 Identities 

Let I=[Q,q] for some {)<q<oo. The derivation of the energy identities crucially depends on the following 
model problem: 



Ut-M^^T on fix I 


(2.62) 


U^^U+U^ on F, 


(2.63) 




(2.64) 


[dM]t^V + V on r. 


(2.65) 




(2.66) 



We should think of Z^^, V, G'('(/;,x) s-nd lT'{ip,x) as certain kind of 'error terms' that will not contribute to the 
extraction of the energy. We first define the stable case energies and (for the model problem) through: 

£,{U,cu) := i ^ [U' + \VU\'} + i {iu^g-u;'} + J^^ {col^-ul, } (2.67) 



and 



V, {U,u) ■.= J^[u! + \VU\'}+J^^ [u;l - } + f {^ht - ^Ih } ■ (2.68) 
As to the unstable case, we define the energies M.^ and by 

M,{U,u) :- [U^ + \VU\^] + i£ {loI ~ |P2+cc>P} + {ul, - \^2+u:gA^). (2.69) 

KiU,Lo):^ f {Uf + \VU\^}+ f {Lol-\¥2+LOt\^}+e f H% - |P2+c^e=tP}. (2.70) 
Remark. Note that the first modes of a; are included in the energy quantities A^e and A^. 
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Lemma 2.2 The following identities hold: 

(a) ^ 

+ ^ P,(s)ds = £,(0) + J {U+Ut)T + j [O + P^+J^ ^ {Q + Sy, (2.71) 

(b) ^ 

Meit) + K{s)ds^MM + j^{l^+l^t):F + j {0'' + P"^+J^^ ^ |q" + 5"}, (2.72) 
where (recalling U.8\) ) 

= 0'^ = 0(U):^U^U^+UV + V{U^f, (2.73) 
p = p« = p(^).= (5*^±_[^^a^]^±)^±+5*^l)_^V"|V^±|2^ (2.74) 

Q^Q{x^^,i^^F)^ -ujt{^ - x) +i^et {^e - Xe) - £^92^ (^92 - xe^ ) + ewes^ (w^a - xe^ ) 

+X^t{y^ - l)+xe-^t{j^ - l) + exioeH{^ - 1) + exe^^e-t(^^^^ - 1) (2.75) 

+hUodp\g\-G{uJtT^ + err^) I5I ; - Q + (Po + Pi )^(Po + POcc-t + enL0g2f>iiug2t . 

\g\F \g\F' 

S^S{x,uj,ip,F) ■.^-ujt{uJt-Xt)+(^et{i^et-Xet)-£i^eH{^eH-XeH) + <^^eH{'^9H-XeH) 
+XtUJt{j^ - ^) + XeHtot{j^ - l) + eXtUJeH{^ - 1) + exe2tWe*t(^;^^ - 1) 

-{~x{F~')t-Xe<^), + Gt)Vocl>\g\ + {f-^+-^)h\g[,^^ 



2 



(2.76) 

Proof. Clearly, part (b) follows easily from part (a). Hence, in the rest of the proof we only focus on 
the identity in the part (a) of Lemma 12.21 The first energy identity arises by multiplying p.62p by U and 
integrating over fi. We obtain: 

]-dt lu^~ I V{U^f+ [ \VU\^- [ U^Ui^ [ UT. 
2 Jii Jr Jq Jr Jn 

Using and dUS]): 

- ju^U^^~ juV- j^{U^Ut+UV) (2.77) 
We shall evaluate — j-^UV by transforming it to the integral over S"^^ using (|2.64p and (|2.66|) : 

We proceed toward the extraction of the energy contribution from the first integral on RHS of (|2.78p . Note 
that 



-e / e 



,.'F ' i;\g\F'' '\g\F ' \g\F"'' ' F^ U \g\F^ UF^ U ^F' 

-dt [ {Ljl-u;^} + ledt [ {ujI.-ujI,} 

f f f (2-79) 

+ / u;t{uj~x)~ ^et{<^e-Xe) + e. I ^e^ti'^e'^-Xe^)-'^ I ^eni^e^-Xe^) 



-/^Xc..(^-l)-^^X.-*(^-l)-./^X-.^*(^-l)-.^^X.^c....(^-l) 

The identities (P?77)) - ^(TT^ give the formulas (PTTH]) and ^^^7^ for O and Q, respectively. The second 
energy identity arises by multiplying (j2.62p by lAt and integrating over fJ. We obtain: 

Jn ^ Jn Jr Jr Jn 
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Using and (PTBS)) 

-jutut ^-j^ {d;u + d;u^ - [^A]u^)ut = - ^ d;uv - ^ { {d;u^ - [7, dt]u^)ut + d;uv] (2.8O) 

We shall evaluate — j-^dlUV by transforming it to tlic integral over using (|2.64p and (|2.66[) . By (|2.64p 
Thus, we have 

-/.r«v.-Xa.(«.*)Vo*,^-X(-f-^-x(l,.-,».(-^)..o.)(-^-;^..)M 

■fA [ XeHt^t f Xt(^BH f XeH^eH , [ / A. ( \ ^ n \^, m \ 



The identities (PTHD]) and (PIM]) imply the formulas (^7^)1 and ^{^J^ for P and 5* respectively. 



(2.81) 



3 Stable case ( > 
3.1 Iteration scheme 

In this section we set up an iteration scheme for proving the local existence of solutions to the Stefan 
problem (jl.lOp - (|1.14p . We want to find a sequence of iterates {u^ ,T^)kefi^ where each curve F'^ is expressed 
as a graph over via the mapping [O,oo[x§i ^ T'^, (j)'' (t,9) = {R''{t,9)cose + xo, R''{t,e)sm9). 

The iteration is defined as follows: given (u"^,R'^) for some meN and smooth initial data u''^^^{Q,-) — 
uo(-), find by solving the equation 

atU™+i-Aw"+i==0, on f^r-, (3.82) 

with the given Dirichlet boundary condition 

After we have obtained u™^-^, for given smooth initial data i?™+^(0, •) = i?o(-)i find by solving the 

following equation on F™: 

+eA™+i = [a„u"+i]t , (3.84) 
where 1/"+^ : F" -> M is given by o = -i^l^ and A^+^ : F™ M is defined through A"+i o = 

-^fpry- Here \g"'\ = ^(P'")2 + (i?™)2. The relation ([STSi]) can be expressed as an equation on S^: 

rm+l r,m rm+1 

- I -e4^ = [g„u"'+i]+o0'" on S\ (3.85) 

Our goal is to prove that the sequence of iterates (w™,/™) converges to the solution of the regularized Stefan 
problem. Applying the temporal differentiation operator dfk (0<fc<Z — 1) to the equations p.82p . (|3.83p 
and (|3.85l) . we obtain 

u™+i-A^ri = on O™, (3.86) 
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m+1 -m I ~,m 



uT+uT on r-, 



where 

and thus, (|3.88p can be rewritten as 



■ + G^ on 



where 

We can thus rewrite 
where 



[5„w;^+i]t=<+< on r'^ 



/■m+1 nrn f 
J fk + l -fi- J/ 



m+1 



|5" 



|5" 



+ on §\ 



k-l 



1 



m+1 



9=0 



k-l 



um _ im ,f,m . im _ \ ^ /-m+lo 



q=0 



fe-i 

g=0 



1 



m+1 

e*t'J+U|^m|/t'=— J' 



(3.87) 
(3.88) 

(3.89) 
(3.90) 
(3.91) 

(3.92) 

(3.93) 
(3.94) 

(3.95) 
(3.96) 



In (|3.95p above, N*{f) is the quadratic remainder in the decomposition iJ(/) = 1 — /— -^j^ + iV*(/), i.e. 
iV*(/) = — (1 — /) — (jij)^^. As a first step toward the introduction of the energy norm for the iterates, 
we state the version of the conservation law (|1.22p which takes into account the iteration. 



Lemma 3.1 The following conservation law holds: 



§1 



ft 



m+1 



m rm+l 



rn 



(3.97) 



Proof. Upon integrating p.82[) over fi™ and using integration by parts we arrive at 



,m+l 



[dn 



,m+ll+ _ _ 



§1 



m + l /-m+l 



|5" 



|5" 



/■m+1 



rm /"m+l 
/ Jt ■> 



which imphes the claim. 

From p.97p . we immediately see 



For any fc e N define: 
where 



(3.98) 
(3.99) 



£e,+ (Z^,w)=f(+)(W,^) + ^ Ve / {u!g3tk~LJg2tk}; V^^+{U,uj)=Vt^+){U,uj) + y^e ( {^034^+1-^92(^+1} 
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Setting 

U^dt,u"'+\ u}^K'^+\ x^Rl", ij^RJ", U^u-^. U^u"^, V = Vr, V = 5^, G = G;", h^K^ , F=l , J^ = 
the identity (|2.7ip implies: 



(3.100) 



£ll+\t)+ f V''^+\s)ds = £,AO)+ f I [0'^ + P"'}+ f / {Q™ + 5"'j. (3.101) 



Here O™ := Et'i 0^^,P^-^ ^1=0 = EtJi QT and 5" = ELo ^fc": whereby 

O™ :=0(a,fcM"'+i), P™:=P(9t.w™+i) 

and 

Qm return pm+1 pm\ cm C f jym jym+l r)m\ 

On the level of iterates, for any j G N we define the 0-th order energy by 



(3.102) 
(3.103) 



^(0) " 2 



(Pu^y + I 



^(0) 



Remark. Note that £^ is slightly different from f (q) . 
3.1.1 The 0-th order energy identity 



To derive the 0-th order energy identity we multiply the equation (|3.82p by u"^'^^ and integrate over fi" 
We obtain 







/ 




Jo™ 





A m+1 



m+l|2 _ 



0. 



(3.104) 



We apply the e-dependent part of the formulas (P?75)) and f^^7^ to evaluate -e /p„K"'A™+^ With x^f"\ 
w = /™+\ F=l, V = i?", 3 = 5™, G' = A^*(/'") and /i = 0, we obtain 



where 



/■m+1/ /•m+1 -fm-i , /•m+l/zm+l rmN , ^/•mi-m+i/ 



m rm+1 1 



(3.105) 



(3.106) 



where the quadratic remainder term N*{f) is defined in the line after (|3.96p . Using (|3.98p and (|3.105p . we 
obtain from (|3.104p : 



ldt{f {u-^'f-f {r^'r}-f fr'+\dte[ {\fU'-\fU'}^l iv.^+^i^ 



^mym+1^771 _|_ / yr/i+l^ym+1 _j_ / 

We now make the fundamental calculation: 



(3.107) 



j {u'^+^ f- j {r+^ f = j {Pu"'+^)^dx+\n"'\{Pou"'+y - J (P/"+i)2-|§i|(Po/"+i) 



,m+l\ 



§1 



(3.108) 
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Using (|3.108p and adding \dt j^iif^'^^f to both sides of (|3.107p . we arrive at 

■>mH 

(0) 



9t£o™+i+2?"+i=AK'+i,/",/™+i), (3.109) 



where 



We expand 



7J™i?" = l_/^',' + 7Vi(/™), 
where -/Vi(5) = 0(|5p + jgeP + IseeP)- We can conveniently rewrite 



= j ^ {/r+'-/r+'(i-/9'^+A^i(r")) +(.r-.r+')/r+'+/r+v^r'} 



/t Jee + Je ht +U j/t - It ^^U/ ij 

i/ +Jee-Jee )Jt ~ It ^^UJ Jj- 

Using the previous calculation we can rewrite A(u™+^,/™,/™+^): 

^(w ^ ,/ ,/ ^ LW + fee" Joe jft -ft )| 

Upon integrating p.l09|) in time and using p.llip . we obtain: 



(3.110) 



(3.111) 



sro) {*)+ v^o) {^)ds^£^om+ Aiu^^+\r,r+^). (3.112) 

Summing (|3.10ip and p.ll2p . we arrive at the fundamental energy identity: 

1/ J J >/ "J S"'" 

(3.113) 

3.2 Energy estimates 

In this section we prove the basic a-priori estimates for the regularized Stefan problem. In order to bound 
the higher-order space-time energy, we first introduce the necessary notations. We denote := (B'^{u^ ,f^), 
ej" := e'g'^iu^J^) and (Big := ^ + ^ for any j e N, where in the definitions ([2331) and (fr36| ft is substituted 
by . Analogously we define S)^'", Sj"^ and T){g. For any jgN we finally set ^^'^ — Y^q^o^'^'^q and we 
define & analogously. Summing the e-independent energy & and the e-dependent energy & , we define 
._^j.,v j^ ^j,v ^ In a fully analogous fashion, we introduce D-'''', S-''" and Tt^'^ . We follow the notational 
convention of dropping the superscript v when v=\. Let us set 



3,f 
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Lemma 3.2 Let C '■= ^ jirj> 0. Then there exist positive constants a, u, Ca, Ca, o,i^d t^ , such that 
for any 9,9i <9o, t<t'^ such that if Ea.i,{0) < 9/2 and 

i-i i-i I I 

k=l fc=l k=l 

KAt)+ fD^:Ar)dr<9, sup Wfnin^. <Ca{9i+ f 'S^n 

Jo 0<s<f Jo 



fc=l 



then 
i-i 



i-i 



EiiiPo/t"^'iu^+Eii^i/t"^'iii^^^-^™^'' EiiiPo/t"^'iii^+EiiiPi/"'''iii^^^-^™^'' (3.114) 



fc=i 



k=l 



k=l 



k=l 



E:^t\t)^ f D^t\r)dT<9 and sup \\f"^+YH^^ <Ca{Oi + ^ f ^ 

Jo 0<s<t e Jo 



771+1,1/ \ 



(3.115) 



Remark. For any jeN, any Sobolev exponent fc€NU{0} and a function v.H.^ -^M., we shall denote the 
Sobolev norm ||w||if'=(03-±) simply by | |w| Ijiffej-Qj) (recall that || • ||//'t(n±) is defined in (|1.2ip ). 

3.2.1 Proof of ()3.114p (zero-th and the first modes) and second claim of ()3.115p 

Let pi= X and P2 = y be the homogeneous polynomials of degree 1 uniquely associated with the spherical 
harmonics si=cos6' and S2 = sin0, respectively. In other words: Si = Tr|si (pi). Multiply (|3.82p by pioS™ 
and integrate over fi. Here 5™ ■.H.-^ft is a diffeomorphism with the property ^^jr™ — (0™)^^- 



R"'f^+'s, + ef"!+'s,d9^ I uT^'p,oS"'+ I Vu"'+'-V(p,oS' 



f"ni-\-l 



Hence 



(1+e) fr+'s^e 



r''fr+'s,d9+ I uT^'p,oS"'+ I Vu"'+'-V{p,oS 



,m+l 



,m+l 



(3.116) 



Let l<k<l and apply the differential operator d^k-i to both sides of (|3.116p . Applying the Leibniz product 
rule and the Cauchy-Schwarz inequality, we arrive at: 



fc-i 



<CE{ll<+V||L^(0".) + ||V<+i|U.(n™)}+CEll/,T'IU^Ell/r'lU- (3.117) 

j=0 j'=0 3 = 1 



Observe that, due to the assumed bound on f^D^^^^,, we have supQ<g<j 1 |^2i < Ca{9+ ^Ji-i )■ Hence, 
fe-i fc-i „ 

Eii/tT'iiL^<ii(Po+Pi)riiL^+Eii(iPo+i'i)/t"NU^+(o'/'<ccy'((^ 
l(c(Ji,)^/%i+cy^).v.4. 



The estimates p.ll7p and p.llSp give 



§1 



k — 1 k 

<cE{IK^^'lli=(n'") + ll^"™^'lli=(f^")}+^/^Ell/tT+'llL- 
In the next step we exploit the stability condition C > 0. Note that for any fc > 1 we have: 

(<+i)2- / (/"+i)2= / {Pui:+'fdx+\n\(Pou"i+'f- f (P/,T+')'-|§+Po/,T+')' = 



(3.118) 
(3.119) 



(Pu'^l^ydx^ / (p/; 



rm+l\2 



§1 



+l^2 



(3.120) 
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To bound the last term on the right-most side of (|3.120p . we shall use the conservation law (|3.97p . Namely, 
applying the differential operator dfk-i to (|3.97[) . we have 



dfk-i 



+ 1 fl'k 



dtk-iinrn 



(3.121) 



To extract the highest order contribution from the first term on LHS of (|3.12ip we observe the following 
differentiation rule for a function U defined on some moving domain 

dt I U:^ [ Ut+ I VU, (3.122) 

Jn(t) Jn{t} Jdn{t) 



where V stands for the normal velocity of the moving boundary dQ,{t). Using (|3.122p (fc— 1) times consecu- 
tively, it is easy to prove the formula 



fc-i 



where -Bfc(Fr™ ,u™+i) := E^Zj (/p,„ T^r^ K,'+^]t) . Hence, from UjJJm and (IXT^ . we obtain 

= dtk-.{r.fr+')-Bk{Vr.^,u"^+^) = -B]!^, (3.124) 



where /si dt^-i {r fr+') + (Fr^ , 



,m+l\ 



Hence, using (|3.124p . we obtain 



m+l 



In ' Jfi 

By the product rule and the Cauchy-Schwarz inequality: 

I „ I fc-l k 



fc=l "^^^ k=lj'=0 3 = 1 3 = 1 



L2. 



In an analogous fashion we show 



fc-i 



iSfci<cEii/t"iiL^Eii^ 



m+l 1 1 



fc-i 



Furthermore, from (|3.124p we obtain for any l<k<l 

iipo/,t+'iIl2<ciiu™+1|u.(o)+c / dt.-.{rfr+')+B. 

Now we may use p.l26p . p.l27p and the definition of _B™ to conclude 

I I „ I 

Ei^ri<E| / dt.-.{rfr^')\+j2\B,\<Cf,Y^ 

k=l k=l •'^^ 

Using this observation together with (|3.125p . we conclude that 

k=i ■'^ -^^^ k=i ■ - 

-fCM^(Eii/r'iii^+Eikr'ii?.MO' 



■;:+'\\L2+C^,J2\\u^+'\\HHn^ 

k=l k=l 



k=l 



k=l 



k=l 



(3.125) 



(3.126) 



(3.127) 



(3.128) 



(3.129) 
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Analogously we prove ^ Z — 1) 

E|(/<n'-(/ rr')" <^^(A + -^)E{lKr^lll,MO"-) + ll/r'lli4- (3-130) 

fc=l -^^ -^S' fc=l ^ ^ 

Note that (|3.126p . (|3.127p and (|3.128p in particular imply 

Y.\\^of:rYL^<cY,{\\uirYLHn'^^)+mf^^^^^^ (3.131) 

Having estimated the last term in (|3.120p . we now deal with the unpleasant negative term — /gi (P/^*^^)^. 
First fundamental observation is that the Wirtinger inequality implies: 

^et-p- [ \FL,t>^\^>Cw\\r2+uJeM\h{giy 
J 

Even though this would suffice to establish positivity of the energy, it features only the second and higher- 
order modes of uj on the right-hand side. However, the inequality (|3.117p gives the obvious bound on the 
i^-norm of the first modes of /™^^. Namely, using the smallness of 9 and v, together with (|3.13ip . we 
obtain from (|3.117p : 

I I 

T.\\^^f^^'\\U^n^^T.{\K^^'\\Un-)^ (3-132) 
fc=i fc=i 

With 6 small enough and Ca larger if necessary, from (|3.13ip and (|3.132p we immediately deduce the second 
claim of p.ll4p . An analogous proof, using (|3.130p instead of (|3.129p gives the first claim of p.ll4p . 

Proof of the second claim of (|3.115p . Note that the equation (|3.85p can be rewritten in the form: 

fr^'+efX' = \r\[uT^]t-fr^'r- (3.133) 

Let q<2l — A. Apply the differential operator dgq to the equation p.l33p . take the square, integrate over §^ 
and then integrate in time. We readily obtain 

r(ll/;,?'lli^+26|lO\|li.+e'||/,"i\||i.) <2 /* / [a,,(|g"|K+i]t)]%2 f [ [9,,(/r+V")]' 
Jo ^ 'Jo Js^ Jo JSi 

(3.134) 

Using the Leibniz rule, L^-L°° type estimates and the trace inequality, we easily conclude that there exists 
a constant Ci, such that 

f\\fr.tH\\h<Ci sup ||n|?,.-4 (3.135) 



On the other hand, since 9t||/^^^||L2 < ll/e^^^^lli^, using the smallness assumptions in Lemma 15^ we easily 
deduce 

e'|lOi|li.<e'||/,™i^(O)||i.+CiC^(0i+0)t Ts^+i. 

Jo 

Thus, using the above for any 0<q<2l — 4, assuming that | |/o| |^2i < -ft'^i for some constant K>0 and 
choosing 9,6i such that {9 + di)Ci,Kdi < 1 suitably small, we obtain 



Jo 



t 



(3.136) 



and this proves the second claim of p.llSp . If we apply the operator dgg with q<2l — 2 to the equa- 
tion (j3.133p . in the same way as above we get (with no e!) 

\\r+'\\W-<CA{0i+t Td^+i) (3.137) 

Note that the same proof also gives the bound 

\\r+'\\l2.-2<CA{e,+t' sup €™+i(.s)). (3.138) 

0<s<t 
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3.2.2 Positive definiteness 

Let 77>0 be a small number to be specified later. Using the definition (|3.99p of and combining (|3.120p . 
P.129p and (|3.132p . we have 



f^^^ '-Jn J ft Jn 

k 



After regrouping, we obtain 



fc=i 

k=i j=i 

k=l 

First choosing rj and then A, 0, /i and e small, we obtain 

I 

V^-+'>C''Y,{\\u^+YLHn) + \N^T^^'\\lH^^ (3.139) 
fc=i 

Let us denote 

I 

fc=i 



Using (|3.139p together with (|3.13ip and the definition p.99p of , we conclude that there exists a constant 
C > such that 

1 

C 

Similarly we introduce 



I?,i(u'"+\/™+1)<I?"+i<CI?||(m™+\/"+1). (3.140) 



/-I i-i 

£ll{uJ):^Y.{\\''t^\\lHn-) + \\^^t^\\lHn^-^^ (3-141) 
k=a fc=i 

A proof analogous to the proof of (|3.139p . which uses (|3.130p instead of (|3.129p leads to 

/-I 

^™+l>^™+l^^^^|||^,n+l||2^^^^^^||^^™^ (3.142) 
k=l 

Like above we conclude that there exists a constant C > such that 

II (^"+1,/™+!) < 5™+! < II (u™+\/™+i). (3.143) 



20 



3.2.3 Temporal energy estimates 

The last and the hardest estimate in the statement of Lemma 13.21 is the first claim of (|3.115p . We first 
estimate RHS of the zero-th order energy identity (|3.112p . We then collect a number of useful estimates in 
Lemma 13.31 and then use it to estimate RHS of the higher-order energy identity (|3.10ip . 

We start off by estimating RHS of the identity p.ll2p . Note that ^(m^+iJ^J^+I) is given by the 
formula (|3.11ip . The first term is a cross-term and is easily estimated as follows: 



t t t 

<X [ I?™+i + ^Ca(26Ii+ / D"+ / D 

Jo "'0 -'0 





m+1 



(3.144) 

where we used the second claim of (|3.115p and the assumption that t<l. As to the second term on RHS 
of (|3.11ip . note that from the assumptions of Lemma [3?2l we have 



Js^ 



Cot 



Js^ 



{frr+^ sup iiiVi(r")iii.<A/ 



1 Co^At 



A 



{01+ f D™)' 
Jo 



A o<s<t 

(3.145) 

In order to deal with the third term on RHS of p.ll2p , we first use the conservation law p.98p to expand 
that expression: 



dt 



.7^1 Jam ./si 



rn+l 



m+1 



§1 



,m+l 



§1 



m+1 



Using the inequality between the geometric and arithmetic mean and (|3.146p . we easily deduce 



< A 



Jgi 



H — — sup [\\u 

0<s<t 



m+l||2 



L2 ^ 



(/r+'f+A 



Jn" 



<A fv"'+^ + ^ sup ^"^+\s) + ^Ca{0i+ f^"'+ sup (£"+1)2(5)), 

Jo A o<s<t A Jo Jo o<s<t 



(3.147) 



where we used both p.l36p and (I3.138p . combined with the assumption t<l. The last term on RHS 
of (IXTT^ is /o /gi A™, where A™ is given by (|3.106p . The first two expressions in the definition of A"'- are 
the cross-terms. We first use the integration by parts and then reason similarly to (|3.144p : 



-Jen Ue2 -Je^) + Je^t U, - 



03 



< 



"'§1 



ftn+lt £711+1 I ptn+1 f -f rii i £m\\ 

h^t \J0 +■'03 -{le +Je3)) 



A(^ J^\\fo'it'\\h + ^i^l \\ii~^i)n\H-^+^l II(/-Pi)/™+^II1^3)<a(^ ^'"^0 

{01+ [ D™-f / 
Jo Jo 



CoCAt 



(3.148) 
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As to the third term on RHS of (|3.106p . we start by integrating by parts: 



1) 



<e 



rm rni + 1 / 
/O "'§1 ^ 15 



1) 



"'§1 



-^e sup (IK^;;^ 

c* ft 



li.o||/,".1li.)<A(/ + ^ sup (||P/"||?,.6|re||l,0 

A 0<s<i 







<A( / 2?™+i + e / sup \\r\\mCA{e,+ I 

Jo Jo ^ 0<s<t Jo 

(3.149) 

where we recall the definitions (|2.46|) and (|1.33p of T)'^ and V respectively. In the last estimate we used the 
smallness assumption on ||/™||jLf2! of Lemma [321 Finally, the last term on RHS of (|3.106p can be estimated 
analogously to p.l49p and we obtain: 



Jsi 



sup \\r\\H^CA{Oi+ /*!)"). (3.150) 

0<s<i Jo 



The following lemma states a number of inequalities that will be used in estimating the remaining error 
terms. 

Lemma 3.3 The following inequalities hold: 

(a) For any 0< fc< / — 1 the following inequalities hold: 

f ||[7™,atfc+i]u™+i||i2(r™)<C sup €™(s) / f sup £"+i(s) 

Jo Q<s<t JO JO 0<s<t 

/*||[7'",9t.>™+i||l,.(r™), t\\b"'A^]u::+'\\Ur^j<C sup £"(5) f ^"^+\ 

Jo Jo 0<s<t Jo 

(b) For any 0<k<l~l, the following inequality holds (df is defined in \2. 59(l ): 

f m.V^^+'\\l.(r-)<C sup (l + (S'"(s)) 

Jo 0<s<t Jo 



-im+l 



(c) Let ^gC°°(J,M) and J QM. an interval such that every derivative of ^ is uniformly bounded on J. Let 
<r <l. Then there exists a positive constant C such that 



< 



C{V£ 



m 1 ^/ifin \ 



and 



C 



\\dtr[mn^+UT?)\\\H^<^ws^' 



(3.151) 
(3.152) 



Proof. Recall the formula (|2.6ip : 



(-1 



m+l 



p—Oq-\-r—l—p 

Note that for any l<p<l~l | |5fi-i-p (F'Vi^ o 0™)| 1^2 < C^C™ and for any g,reNo such that q + r^l-p, 

we have 

||9,.5„,u™+i|U.o(r™)<C||V<p'+i||H<-P(r.")<cVs^, 

where we use the trace inequality and the fact that Z — p+ 1 < 2(Z— 1 — p) + 2 (observe (|1.37p ). If p = 0, then 
||5ti-i(F'?I^^o0™)||i2 <Cr"V^ and from the trace inequality, (fT35|l and the fact that q + r = l: 



I |9,.a„,u™+i I Uoc(r™) < c\ I |ff<+i(r..) < cVe^, 
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where we note that ||u™~''"'^||/fi+i(r'") < CV if l>2 (observe (|1.35p '). Combinmg the above inequahties, 
we deduce ^ ^ ^ 

/ ||[7™,atfc+i]u"+i||i2(r™) <C sup €™(s) / /" J)™ sup €™+i(s) 

Jo 0<s<t Jo Jo 0<s<t 

and this proves the first claim of part (a) of Lemma 13.31 The other two claims are proven similarly. As to 
the part (b), fix < fc < ? - 1. Then 



< 



^ r>m pt pt 

^'^p / IIC++i'lli^<c(i+ sup / V 

^0<s<t |ff"r* Jo 0<s<t Jo 



m+1 



and this finishes the proof of part (b). Note that we used the Leibniz rule in the first inequality, and 
then standard L°°-L^ type inequality, combined with the Sobolev inequality. Part (c) follows using similar 
techniques and the definition (|1.34p of € and D. This proof is in fact analogous to the proof of part (a) of 
Lemma 3.3 in [21]. □ 

We proceed now to estimate RHS of the energy identity p.lOip . 
Estimating O™. Recall that is defined by (|3.102p and 



In order to estimate the first term on RHS of (|3.153p . we apply the relation (|2.6ip to write 
We next note that from p.49p 



(3.153) 



(3.154) 



,m+ll + 



/o Jr 

< 



< 







where we used part (a) of Lemma 13.31 in the last line. Furthermore, using p.84p yields: 



(3.155) 



JT" 



JT" 



c 
< — 
- A 



/o Jsi 



where A™+^ is defined in the line after (|3.84p . Note that by parts (a) and (b) of Lemma [3731 



2) 



m+l 



(3.156) 
(3.157) 



Furthermore, recall the formula dtk ( 



+ K,2 iK,2 given by ((3J6l) ). Hence 



Js^i 



m+l||2 



l_f/i(r™) " 



Jsi 



Jsi 

Jo Jgi 



(3.158) 
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Note that by p.96p and integration by parts the last term above is bounded as follows: 



/o Jsi 
■- e 



Note that from 



(3.159) 



Similarly eP^^ J'"+i(^)^,_J|i2 < CV€^VP™+i. From ((3A59)) . (|3l60l) and the estimate after ((3A601) . 
together with part (a) of Lemma 13.31 we deduce 



C [ [7",a,.]u"+io</,™/^^2|5'"| <C sup /*||[7'",9,.>'"+i|| 

Jgi 0<s<t Jo JQ 



m+ll |2 

Hi(r" 



(3.161) 



<C sup / I?™+^+C sup €™ / P 

0<s<t Jo 0<s<t Jo 

From p.l56p - (|3.158p . (|3.16ip and part (b) of Lemma [5751 we finally conclude 



[7'",94.]u™+ia;.K+ii 



JT" 

-C sup / V"'+' + C sup r" / V 

0<s<t Jo 0<s<t Jo 



t rt 
m I -Am+1 I ^ rem I -j^m+l 



<CVe^ I +CeA/ V^+'+Ce^X V 

Jo Jo 



(3.162) 



Estimating /g/r™^"- Recall that P^' is defined by (|3.102|) and ((2J4)) : 

n"=(at([7'^at.]«"+l'±)-[7",9t]<+'■^)9„^.;:+^'±+9;.+,^.™+M7'^5,.]K'+']^+^^^ 

It is, however, easy to see that the expression in the first parenthesis above simplifies, i.e., 



m+l,± |2 



m+l 



Using this observation together with (|2.60p applied to [dnU^ti], we may rewrite PiU)'- 

P>m [ m cx 1 m+lo m+l,± . o* m+lr mo ir , it-? m+l.±|2 

J, =[7 ,dtk+i\u ^ dnUjfc +04^+11* ^ [7 ,dffcj[u,j JI + l/r™|Vu ^ ■ | . 

We start by estimating the first contribution on RHS of p.l63p . Note that 

r'c. 



(3.163) 



[7'",54.+i]7x'"+ia„<!+i-± 



JT" ^ 

<^ sup / 1)™+! + ^ sup €"+1 / I?" + CA / 

A 0<s<t Jo A o<s<t 



lL2(r" 



(3.164) 



where we used part (a) of Lemma 13.31 and the trace inequality in the second inequality. As to the second 
term on RHS of (l3l63ll . first recall that dj^+itt'^+i o0™ = - + 0^+1 i^T+i defined by fjSMh ). 
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First note that 
ft 



JS^ 



Jo ^ Jo 



c 



(3.165) 



where we used part (a) of Lemma 13.31 in the last estimate. Integrating by parts and reasoning in the same 
way as in (j3.165p . we obtain 



lo Js^ R^gn 







(J 

— ve™ / D 

Jo 



?n+l 



Note that by ^jSM^i 



(3.166) 



Jo JS^ 



< 



JS^ 



q=0 



JS^ 



(-^^) [7™, 9,. ] o 0™ I 



(3.167) 



The first term on RHS of p.l67p is easily estimated by the Cauchy-Schwarz inequality and part (a) of 
Lemma ([33]) (recall p^M]) ) iN*{f) defined in the hne after ^Ml): 



iV*(r),.+,[7™,a,.][<+i]lo</>-|g-| < / ||7V*(/™),.+i||i.+ / ||[7™,5,.]K+il+l|2 



Js^ ^ JO^ 

<C sup (||/'"(s)||^3+(£™(s)) / V'^ + cVe^ [ 2)"+^ 

0<s<t Jo Jo 



(3.168) 



As to the second term on RHS of (|3.167p . we first use the integration by parts, and then the standard L'^-L°° 



type estimates in analogy to (j3.160p . For any 0<q<k (recall d^ = dgdf): 



-)[7™,a,.]K+i]lo0"%"M 



< 



1 



Js^ 



1 



,0 jJ'^"^^'^'~'^'R^-W^ 



MW^\d,.][u":+']tor\9"'\) 



<C sup £'"( / V" 

0<s<t Jo 



£1 



(3.169) 

where we used part (a) of Lemma [3731 in the last estimate. The last term on RHS of p.l63p is estimated as 
follows: 



r / Fr™|V<+i-±p<C sup ||yr™|U=.(r™) < C^e^ / 

Jo JT"^ 0<s<t Jo Jo 







(3.170) 



Estimating /g/siQ"- Recall that is defined by (|XTU51) . where Q is given by with (jXTUOl) . The 

first four terms on RHS of (I2.75P are the cross-terms, first of which is estimated as follows: 



Jsi 



<A^'ll/;iv+i'lli. + y_^\ll/,T+^lli. + ll/,T|iiO- (3-171) 



The second cross-term is estimated analogously: 



Jsi 



C 



The third cross-term is estimated similarly: 



(3.172) 



Js^ Jo Jo 
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As to the fourth cross-term, in the same way we show 



Jgi 



C 



(3.174) 



The fifth term on RHS of p.75p is estimated as follows: 



"'gi 



0<s<t 



< sup \\ris)\\L^ / {ii/,T,+.^iii.+ii/,Tiii.}<c sup \\ris)\\H^{ / / p™) 



0<s<t 



(3.175) 

In the sixth term on RHS of l|2.75p with p.lOOp we first integrate by parts and then estimate in the same 
way as in (|3.175p : 



Jg^tkJi 



"'§1 



■-1 



< 



1 



Jgi 



■-1 



<c\\ 



{ I If^ti I lis + 1 l/g^fc I lia } 



<C sup {\\r{s)\\HS + ^¥Hs)){ / / v^), 

0<s<t 



where we observe that 



^)e\\L^<C\\n\w^,^. 



(3.176) 



Note that the seventh term of (|2.75p with p.lOOp is simply and since F=l and the eighth term on 
RHS of (|2.75p is estimated in the same way as in ()3.176p . After integrating by parts, we obtain 



1 



j?m rm-\-l / 
j^J'^'^'^SH^+'^'R-^lg- 



<C sup ||/"(s)||ff.( / / 

0<s<t 



(3.177) 



Note that next-to-last term of /gi Q^, given by (|3.103p and ([275]) (with (I3.100p and ^^), involves 



the term h'^^hf^+eh'^^, given by (I5T961) . Note that 



Jo Jo 



JS^ 



Through the standard L°°-L^ type estimate and Sobolev inequality, we obtain 



(3.178) 



where we used Lemma 13.31 part (c), in the second estimate. Observing that u™ = 9*fcU 
[7™,9jfc]u™+^ we estimate the first term on RHS of p.l78p as follows: 



m+l „ 



lO JSi 



hT.i^T-nn < / ||/i^il|LHKIlL=(r".) 

t t (3.179) 

Jo Jo Jo 

where part (a) of Lemma [3731 has been used in the last inequality. As to the second term on RHS of (I3.178p . 
in the same way 



Jsi 



<C sup €™ / I?'"+^ + C sup r" / V"'+\ 



Jgi 



0<s<t Jo 



0<s<t Jo 



(3.180) 
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where we used the inequahty p.l6ip in the last estimate. Thus, combining p.l78p - p.lSOp . we obtain 



Jq Jsi Jo Jo 



-,m+l / rT\m+l\ 



(3.181) 



The last term of /gi Q™, given by (|3.103p and (|2.75p . involves the term G™, which is defined by p.95p 
with (|3.100p . Using the standard L°°-L^ type estimate, Sobolev inequality and part (c) of Lemma [3T3l it is 

easy to show that 

||G^||l= <C(||Pi/"||l~ (3.182) 
Thus, by the Cauchy-Schwarz inequality, we obtain 



Jsi 



Gm pm+1 jDrn 
k /ffc + l ^ 



<q|GriU2||/;r+i|U.<c( sup ||Pi/"||L=o + ye^)( fv"'+'). (3.i83) 

0<s<t Jo Jo 



Furthermore, integrating by parts and applying the Cauchy-Schwarz inequality we obtain (recall d^j^-^ - 



Gma4 rm+1 
k '-'k+lJ 



J& 



<e 



Jsi 



L2- 



(3.184) 



Recalling the definition (|3.95p of G™, we have 



q=0 



(3.185) 



Using the Leibniz rule, L°°-L^ type estimates, it follows that {N*{f) is defined in the line after p.96p ): 



Further, for any < < fc — 1 



(3.186) 



1 



Similarly, for any < q < fc — 1 

k — l 



1 



fc-i 



.-^\L^<^reY.\\^lr\\L^ 



1 



Dm \fjm I / dt^ 
q=0 ^ 1^ I q=0 

fe/2 k-l 

<gE(||/™||h3 + Vp^)v^+G E (x/7||n|H3 + v^)Vs 

q=0 q=k/2 



(3.188) 



Hence, from p.l85p - p.lSSp . we obtain 

yfe\ \{Gk)0\\L^<C{\\r\\m + \/e^) {VV^ + 7^ + V^') + GP™ + GP'^ 
Combining p.l83p - (|3.189p . we conclude the bound on the last term of /J /gi Q™: 

/* / Gr(i?"/™+i+6/;;j+^0 <G( sup \\r\\^, + ^){ f vT+ fvT+') 

Jo Jsi 0<s<t Jo Jo 



(3.189) 



(3.190) 



Estimating /o/si'^T- Recall that 5^ is given by l|3.103p and ((2J6l) (with (|3.100p ). The first four terms 
on RHS of (|2.76p are the cross-terms: 

lo JS} _ , _ Jo Jo Jo ^ Jo V 



<X I V"' + ^ I + sup £ 



C 



?m+l 



0<s<t 



(3.191) 
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where we used the relation (|1.3p and the trace inequahty in the last inequality. Similarly, for the second 
term, we have 

C 



Uetfc + i J gtk + l)J Qf-k + l 



< 



Jn 



Q_ 

c 



\f> 



m |2 



I fm+l 1 2 



+ -A 



V Jo Jn 



|V2^m+l|2 



(3.192) 



As to the third cross-term on RHS of (|2.76p . we have 



<A / P™ + -^(ry||w™+i||2^, + ^||u™+i||2^,)<A / + ^ / + sup £ 



(3.193) 



0<s<t 



To estimate the fourth cross-term on RHS of (|2.76p . we apply the same estimates as in p.l93l) : 



Jgi 



rm+L / prn-\-L rm \ 

J g3fk + l\J g3fk + l J gSfk+l ) 



91 

eX 



C 

1 sup 



(3.194) 



In full analogy to the estimates (|3.175p - (|3.177p . we can estimate the fifth, sixth, seventh and eighth term 
on RHS of (|2.76p (note that the seventh term is actually 0, since F=l). We obtain: 



Jtk+iJtk + l J 



JSi 



rm fva+1 ( 
Jg2fk + lJ^k + l [ 



^9" 



Jsi 



rm rrn+1 ( 

Jg2fk + ljgi^k + l \ 



<C sup \\r{s)\\H^[ I V\ 

Q<S<t 



,m+l 



(3.195) 



The next-to-last term on RHS of Jg^ S*™ involves the expression G™ which is given by (|3.95p and hence 



fe-i 



fc-i 



)t — ~ / Je^ti^ 



q=0 



HfoH->{^;;:^)tk-,+i+dtk+iN*{D 



Separating the cases q<k/2 and q > fc/2 and taking and _L°°-norms accordingly, we conclude 



\\{GT)i\\h<cmir\\i^+ sup €^(5)) s)^ 

0<s<t 



(3.196) 



Hence 



where we used the trace inequality, part (a) of Lemma 13.31 and the estimate p.l96p . 
Let us introduce a small parameter 62 such that 



sup {E{s) + <t{s))+ I {V{s)+'£){s)]ds<d2. 

0<s<t Jo 



(3.197) 



(3.198) 



We shall later express the desired constant (from the statement of Lemma [3T2|) with the help of 02- Using 
the identity (|3ll3l) and estimates (l3T44ll - (13^50]) . (13^55]) . ((3T621) . ([3164]) ([STTO]) . ([STTT]) - (ISTTT]) . 
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(|3.18ip . p.l90p and (|3.19ip - p.l76p . we arrive at: 

Jo_ _ Jo Jo , , 

+(^+-^)f'"+i+^(fr+^r+')+^^(A+ sup ii/"iiH^+\/e^+€™)( / / 

eryA A o<s<t Jo Jo 



m+l\ 



(CA+^ + ^ + C7fJ/')^*D'"+i+C sup 



0<s<t 



(3.199) 



+ (CA+^ + ^ + C0^/V^) f^"^+' + Ce, sup €^ 
A eA A 7o o<s<t 

where we used the smallness assumption above. Let us set 



m+l 



C f Cf 

lC{62,t):^ ^{6^+62) + —e2 + C{\ + 62 + ^2)02, (3.200) 
A A 

which shortens the notation for the constant term appearing on RHS of (|3.199p . 
3.2.4 Reduction and representations of J-pVssVn and /pUfW„ 

Just to provide intuition, assume for the moment that (w,/) is the solution to the Stefan problem p.lOp - 
(jl.l2p . The energies € and D arise by applying the space-time differential operators to the heat equation p.ip . 
For instance, after applying dxdt to the equation p.ip . multiplying by —Adxdtu and integrating over Q, 
we run into the trouble of extracting the boundary energy contribution. The difficulty is contained in the 
restriction of the spatial differentiation operator d^i to the moving surface F. Namely d^i =Tids+nidn on 
F, whereby \ni\ is of order 1 and hence the normal component of d^i can not be assumed small. Thus, the 
restriction of the operator d'^ to F inherits non-trivial high order normal derivatives of u. This presents a 
difficulty in extracting the energy contribution, which is overcome through the reduction procedure explained 
in this subsection. The main idea is to recover the spatial derivatives via the equation in the bulk and the 
established temporal energy estimate in Subsection 13.2.31 We first decompose the Laplace operator on F 
as Au\Y = Usg + KUn + Unn- Since Ut = Au in the bulk, we can represent m„„ in terms of the Ut (which is 
controlled by the temporal energy) and Uss (which helps us extract the boundary energy contribution). We 
now rigorously justify the above heuristics. 

In order to facilitate our arguments, we shall denote the multi-index fi of order 2p, by an or- 
dered 2p-tuple (ii,Z2,...,i2p), where Zfee{l,2} for k = l,...2p. Let c = l—p—l and set v = d^u^^^ = 
dx'idxi2,---d^i2pdfU"^~^^. After applying the differential operator d'^df to the equation (jl.ip we obtain 

Vt-Av^O. (3.201) 

Multiplying the above equation by —At; and integrating over il™, we get J^„^{Av)'^ — J^,„vtAv = 0. On the 
other hand, using the Stokes theorem and denoting dxiV = Vi, we obtain 

(At;)2= / V^wiV^-y- / ViVi.inj+ / Aot„. (3.202) 
Jn™ Jr™u9n Jr^uao 

Recall that dxi —Tids+riidn on the moving boundary. In the following we employ the Einstein summation 
convention and sum over the repeated indices: 

/ ViVijTlj^ / TiVsVijnj+ / niVnVijTlj^ / Vsn{TidsVj +nidnVj)nj + / VnVnn 

Jr™ijaf2 Jr'"u9n Jr-^uao Jr"^uao Jr™u9n 

VsdsVjnj+ Vn{Av- KVn-Vss)^ Vs{ds{Vjnj) - Vjds{nj)) + VnAv 

•uaa Jv^udn Jr™u9f2 Jr^uon 

•u9n Jr^^uan Jv^^uan Jv^^uan Jr^^uan 



/r™u9n Jr'"uaa Jr™u9n Jr™u9n Jr"'uaa 
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In the second equality we used the identity UiVijUj —Vnn- In the third equahty we use r,;ri = |Tp = l and 
Aw^Wss + KWn + Wnn. In the fifth equahty we use the fact that dsUj = KTj together with VjTj —Vs- We also 
substitute Vt for Av due to the equation p.20ip . In the last equality we use integration by parts. Further 
note that 

Hence, going back to p.202p and the remark before it, we obtain 

V\:V^v+ldt [ \Vv\^+ [ n{vl+vl)^ ( vtVn-2 [ VssVn-^ [ Vr^\Vv\^ 
2 jj^m Jr'"u9o Jr'"u90 Jr^uoa 2jp,„ 

Upon integrating in time over the interval [0,t], we arrive at the fundamental identity 



\Vvit)\'+ / \V'v\'+ / n{vi+vi) 

Jo Jr'"uSf2 Jo Jr^uan ^■ 



(3.203) 



Convenient representation of Vss, Vn and vt. Motivated by the identity (|3.203p . our goal is to 
understand better the boundary integral terms /pmygj^Wss^'n and /pmuan^t^i- To do so, we first proceed 
toward a description of terms Vgs, and Vt in terms of tangential and time derivatives, as given by the 
formulas IKTT^i . (IX^ and (1^^ . 

Let C denote either F™ or dfl. Note that for i = l,2, the operator dj.i on C can be written as d^i = 
Tids + riidn- Hence, for a given multi-index ^— (ii,...,i2p) of order 2p, we have 

2p 2p 2p-j 

A:=0 JC{ii.....i2p} jeJ jeJ" j=l a=0 

J|=fc 

Here Faj{T,n) is given as a sum of the expressions of the form C9gm(ii) (t^j Vni^), . . . ,i9^m(i2p) (Ti^^ Vni^^), 

where m{ik) > and = J ((?sm(i) {riVrii) equals either dgm{i)Ti or dgm(i)ni). Each of the summands 

in the second term on RHS of (13.204^ is of lower order, and by this we mean that the number of derivatives 
hitting Mjc equals at most 2p — 1. As to the first term of on RHS of p.204p . we carry out the crucial technical 
step, that we will refer to as the reduction procedure. Given b>2 and cG N we note the following: 

d''^Ut^^d''-^{Aut^~Kd,,ut^~dssUt^)^d''-\t^+i-dssd'i-^Ut^-Kd''-^ut^. (3.205) 

Carefully using (j3.205p . by induction we deduce for b even: 

6/2 /L/o\ 

9„.^.,c=^(-l)™( M92™„,.+,^,_+ ^ Mlp^^{K)d:df,u,.^.; (3.206) 

3=0 V / a+l3+2-t<b 

and for b odd: 

L'/'J f\bl2\\ 

d^.Ut.= Y,{-ir[^ 'J\df^dr.Ut.+ yW2,-„.+ Mlp^^iK)d:d^Ut.^-,. (3.207) 

TO=0 ^ a+p+2j<b 

For a given triple (a,/3,7), the function ^ ^{n) is a sum of the expressions of the form Cdgii k. . .dgi^ k, 
where «i,...,«m>0 and ^|ife|=Q;' for some a' such that a' + a< L^J- This is a counting argument and is 
easily seen by counting the occurrences of k in the repeated application of the formula (|3.205p . To illustrate 
this, assume 6 = 4 and note that in the second iteration of (|3.205p . the term dssd'^^Uf generates the contri- 
bution dss{nd^-i^^ut>:) and from the product rule this expression always has the form 9" ndfd^^Uf , where 
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a + a < 2 = [f J . We can argue inductively to conclude the claim. Using the relations (|3.206p and (|3.207p . 
we may write for k even 

ds2p-kdnkUt<:—''^{~iy^(^^jdg2p-k+2mUfa+b/2~rr,+ ^ Af^ ' ^ ^ (t, 9" 9^ c+t ; (3.208) 

m=0 V / a+/3+27<2p 

and for k odd: 

Lfe/2J 

ds2p-kdnkUf== ^ (-1)^ 
m=0 

Here, M^''^ ^{T,n) is defined analogously to M^p^{k), whereby the total number a' of ds derivatives in 
the summands of the form d^ii k. . .dgim k satisfies a' + a<2p—k+[^\. Hence, using p.204p . (|3.208p 
and (|3.209p . we obtain: 



YLfc/2J 

V m 



ds2p-k+2^dnUt.+ iw2i-,^+ Ml;^pn^T,n)d':d^Ut.+-,. (3.209) 

a+/3+27<2p 



fccvcn|,/|=fc i^-' m=0 ^ ^ 



\k/2\ 



- E E n E (-1)^ 

feodd|J|=fc i^J m=0 



'g2p-k + 2mU^C'\-k/2- 



Q+/3+27<2p 



where 



2p 2p 2p-j 

3 = 1 a=0 



fc=o Jc{ii,...,i2p} jeJ jeJ" 
\j\=k 



With the use of formula p.210p we write out the expression for Vgs on T^' 



(3.210) 



k/2 



fccvcnljl— /c m— 



9ss5^Utc = ^ ^ rtj-Tj/ (-1)™ ^^^^^^ ]ds2p-k + 2m + 2Uta+k/2 



Lfe/2J 



+ E E n E 

/codd|J|=/c jeJ \ ™ / 



T?l = 

fc/2 [fc/2j 

E E n '^i^j'E°"+E E n "j^i' E 

feodd|J|=fc je-f m=0 



fcovcn|,7|=fc i^J m=0 



(3.211) 



where we define F^:— dss{^ci+i3+2-i<2pFa,ii.,-i{T,n)d'^df^Uic+-y^ , and by product rule, we can write it in the 
form 

a+l3+2i<2p+2 

Furthermore, 



(k/2 

am = (— 1)*™ 9s2p-fc + 2m + 2U(e+fc/2-m, ^'m = (^l)'^ 



Lfc/2J 
m 



's2p-fc + 2m + 2 9„Mjc+Lfc/2J 



(3.212) 
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For Vn on F"* we obtain: 

k/2 



fccvon|J|=fe j<^J m=0 ^ ^ 



'„Ujc+fc/2-r 



- "jT?' 5Z j {9s2p-fc + 2mUic+Lfc/2J-m + l — 9s2p--fc + 2m + 2M(c+Lfc/2j-m} +Fp (3.213) 

/£odd|J|=fc m=0 \ m J 

k/2 [fe/2j 

H H n "i^J'Hc,„+^ Jl njTjv ^ (d„,-e„) + i^^, 

/£oven|J|=fe JEJ m=0 /codd|J|=/£ JE-' m=0 



where 



fk/2\ /[fc/2j\ 

C„i = (— 1)"M 932p-fc + 2m9„Ujc + fc/2-m , dm = (— l)'" 952p - fc + 2m U(c+ L fc/ 2J - ™ + 1 , 

\m J \ m J 



Cm — (— 1)™( L ^ j 932p-fc + 2m + 2U(c+Lfc/2J 



(3.214) 



Furthermore, is defined by -F/j = X]a+/3+27<2p+i ^"./3-7(''''"')^?^n^t'=+^ ■ Similarly, using (|3.210p . we obtain 
an expression for vt on F™: 



= ^ ^ JJ^ "-jTj' ^ (^1)™ ( ^ )9s2p-fc + 2mUic + fc/2-,Ti+l 

fccvcn|J|=A: i^' m=0 \^ J 

Lfc/2J 



n %Tj' ^ (-l)"f jas2p-!= + 2,„a„Utc+Lfe/2J-™+l+ ^ Fa.J3,^{T,n)d"d!^Utc + y+l 



fcodd|,7|=fc je-' m=0 ^ ^ Q+,3+27<2p 

fc/2 [fc/2j 

=^ E E n E^™+ E E n "^'^^y E ^™+^/^' 

fceven|J|=fc JE-' m=0 fcodd|J|=fc JE-' m=0 

(3.215) 

where 

= (-l)"P^^V.^''-'=+^'" B™-(-l)"f^''^^JV.^''-''+^"^""t=+LV2j-^+i (3.216) 

\ m J \ m J 

and is of the form = J2c+f3+2'y<2pPa,i3,',{T,n)dfd^Ut.+-,+i. With formulas ((33TT|) . ((33T31) and ([ST^TS]) 
in hand, we can evaluate the the integrals Jp„ UssW„ and Jp„ utt;„. 

Evaluating the expression /p,„WssW„. The goal is to arrive at the identity (I3.223p . which helps us 
extract the boundary energy contribution and isolate the 'error' terms. Using (|3.21ip and (j3.213p . we may 
write 



y2 j VssVn^ dssd^'ut'=dnd^'ut<= 



2 fc/2 Lfc/2J 

= E / ( E E n ^j-^j' E E E n "j^j' E ^ 

ir = i, fccvcn|J|=A: ie-' m=0 fcodd|J|=A: ie.l m=0 

r=l,...,2p j'ejc I I -.^ja 

fc/2 Lfc/2J „ _ . 

E E n "j^i'E'=™+E E n '^j^j' E('^""^™))+ E / "^5-^^+ / w„^^p- XI i 

fcovcn|J|=fc je-' m=0 fcodd|J|=fc jeJ m=0 Ut|=2p |Ai|=2p |Ai|=2p 

(3.217) 

Let us simplify the first integral on RHS of (|3.204p . We observe that if J and J are two subsets of the index 
set / = {ii, . . . ,Z2p} such that J ^ J , then there obviously exists at least one index i^, such that is&JC\J'^ or 
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is G Jn J"^. Without loss of generality assume is £ JC\ J^. In this case, we note immediately 

2 2 
XI n n "^rTr' ^^n,,n, ]J P^Ty J| n^T^/ = 0, 

15=1 jeJ, j' eJ" reJ.r'eJ" «5=i jeJXisJ'eJ" re,hr'eJ''\is 

since X)?^=i"'i5'''i5 = ?T.-r = 0. Because of this, the first term on RHS of p.217p reads: 

2 k/2 k/2 2 Lfc/2J 1^/2] 



/ E EE n "'^'E°™E^™+/ E EE IT ^ME^-Ek--^™) 

ir = i fcovon |./| = fc 7"eJ,7'G./^ m=0 m=0 »r = i fc odd l./N* jG J, •i'SJ'^ m=0 m=0 

re{i 2p} ,/c/ re{i,...,2p} jc-r 

/ecvcn ^ ^ -^^m^O m^O /c odd ^ ^ m^O m^O 



(3.218) 

where we used X]j=i '^J — 5I]j=i j = ^ together with the fact that the number of fc-subsets of a set of size 2p 
is In the first term on the RHS of (j3.218p . we single out the Jf^ak/2Ck/2 contribution, as we shall use it 
to extract the leading order energy contribution (see ()3.212p and (I3.214|) ): 



E Pfc) /E«™^c„.= j2 Pf) / afc/2Cfc/2+ E Pfc)E(-"fc/2Cfc/2+E«™E 

fceven ^ ^ ^ m— m— even ^ ^ ^ even ^ ^ C m— m— C 

E i^k) J ^s''-P^^Ut-=ds2pdnUt'=+ ^ ^^^^ j (-afc/2Cfc/2+ X! "™ E 





fc/2 fc/2 



(3.219) 



k even ^ ^ ""^ fc even ^ / L. m=0 r?j.=0 

Similarly, (see (|3.212|) and (I3.214P ) 

Lfc/2J Lfc/2J 

\k/2\eik/2\ 



i:Cr)ii:''...i:(*..--)-i:(?)i'>L 

codd ^ m=0 m=0 fcodd ^ ''^ 

^ E ( ) / E E '^'"^ E ( f ) / (-^Lfc/2jeLfc/2j + E ^™ E ' 

fcodd ^ ^ r?j.=0 m=0 fcodd ^ ' m=0 m=0 



Lfc/2J Lfc/2J 



feeven ^ ^ fcodd ^ ^ t?x=0 m=0 

E ( fc ) / (-^Lfe/2jeLfcV2j + E ^"^ E 



(3.220) 



Note that 



E ( hi I ds2p+2Uf=ds2pdnUf= — { hi I ds2p+iUf=ds2p+ldnUtc ~~2'^P j ds2p+lUf=ds2p+ldnUt':, 

fceven^^'^-^C fcodd ^''^•^C Jc 



where we integrate by parts in the first integral on LHS and then use the binomial expansion formula. Hence, 
using this observation, we conclude from (|3.217p - p.220p 

VssVn = -2'^^ ds2p+lUt'=ds2p+ldnUt-=+ Ip+ ^ VssF^+ ^ jvnF^- ^ F^F^, (3.221) 

where 

fc/2 fc/2 , , [fc/2j [fc/2j 



'■= E \k) E (~"^/2Cfc/2 + E E + E \k) / E ^™ E 

fceven ^ ^ C m=0 m=0 fc odd ^ ^ C „i=a m=t) 

^ E ( f) / (-^Lfc/2jeLfe/2j + E ^"^ E ^'")- 

fcodd ^ ^ Tn=0 m=0 



(3.222) 
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We use the trace formula (|2.60p to further simpUfy the first integral on RHS of p.22ip : 

dg2p+iUf=ds2p+idnUtc = / 9s2p+i (^dfau-[j"^,dt<:]u)ds2p+i (9t*cU„- [7™,9tc]M„) 
' " . . f f 

ds2p+ldfcUds2p+ldfaUn— / ds2p + l[y"^ ,dt<:]uds2p+ldnUf: — / ds2p + lUt':ds2p+l['y"^ ,df]ur-, 

Jc Jc 



ds2p+i [7™ , 9tc ] uds2p+i [7™ ,df=]un. 

c 

Plugging the above in p.22ip , we obtain the following important formula: 

22p / d g2p+id^,:uds2p+idfcUn + 2^^ / dg2p+i[Y" ,dti:]udg2p+idn 
, „_vc Jc Jc 



Im1=2p 

+2^P I ds2p+iUt.dg2p+i[j"\dt.]u„-2^P I d,2p+i[j"\dt^]ud,2p+i[j"'A^]un (3.223) 



Jc JC 

|m|=2p |Ai|=2p"^'^ |/i|=2p"^'^ 

Evaluating the expression J^^VtVn. Our goal is to prove the formula (|3.227p . Using the represen- 
tations (|3.213p and (|3.215p of and Vf respectively, just like in p.217p . (|3.218p and p.219p . we arrive 
at: 

E / ^tvn^ E fc / E^™E^™+ E I) / E ^™ E E / 

|m|=2p fecvon ^ "'C^^g fcodd ^ ^ ^ m=0 m=0 \fi\=2p '^ 

+ E E / E / Ak/2Ck/2+ E / (-^fc/2Cfc/2+E^™Ec™) 

|At|=2p |/i|=2p fcovon^ "^'^ fc even ^ ^ "^"^ m=0 m=0 

+ j£BraY. idm-e^)+ E / + E / ^"^^ " E / ^^^M- 

(3.224) 

we only concentrate on the case when k is even, whereas the odd k-s are treated as an error term. Note that 

E (t) /,/^/2^''/2 = 22p-i I dg2pu^+,'d,,2pd^ul^i+\ (3.225) 



(^m ^m)- 



Let us abbreviate 

Jv-= nil) (-^'^/2c./2+E^™E^™) + ^ r J2BmJ2^ 

feevcn ^ / JC „j^Q ,„^Q fcodd ^ / -^C m,=0 

Using (|3.224p and p.225p . we arrive at the following identity: 

E j Vtv,,^2^P-^ j ds2pu^+lds2pdnuT+^+ I Jp+ fy^F^+Y^ L^F^^ f^t^^t^- 

lMl=2p-^C Jc Jc |^|=2p-^C M=2p-^<^ \t^\=2p''^ 

(3.226) 

Just like in (I3.223p . we use the trace formula (j2.60p to simplify the first integral on RHS of ()3.226p and we 
obtain, in analogy to ()3.223p . an important identity: 



Y I vtvn^2'p-^ 1 9,2p9;.+iu"+i9,2pa;.u;7+i-22p-i I a,2p[7",9tc+i]u'"+ia,2p9„<+i 

-22P-1 / d,2pU^+ldMTAAu^+^ + 2''P'' I a,2p[7'",9te+i]u"+i9,2p[7",9te]<+i (3.227) 



- fjp+Y I E / E / ^M^M- 

|/x|=2p-^C |;.|=2p-^C |M|=2p-^C 
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3.2.5 Space-time energy estimates 

We proceed to estimate RHS of p.223p . term by term. We start with the last three integrals on RHS 

ofdinni): 

Estimating terms J^VsgF^, J^^VnF^ and J^Ff^Ff^. Recall the definition (|3.213p of and the formula 
Ff_i = J2a+0+2f<2p+i-^o',t^,i(''''''^)'^s''^n''^t''+-y- Let us, thus, fixatriple (a,/?, 7) of nonnegative integers satisfying 
a + /3 + 27<2p. 



i<i 



VssFa,f:i,j{T,n)d"d^Utc^ 



< 



ll"s||ffi/2(c)||^'a./3,7(^,«)9"9fM"i^||i/i/2(c) 



< 



+ /3 



<Cr} 
<Ct] 



V 

c\ 



t<=+i llH»+i;i+i(n™) " 



^ I WVu. 



•^'7 Jo 



m+l||2 



p ft r 



(3.228) 

in the first inequality we use the ^/^-version of the Cauchy-Schwarz inequality and in the second inequality 
we use simply the Cauchy-Schwarz inequality. In the third estimate we apply the trace inequality and in the 
fourth estimate we use the interpolation inequality between Sobolev spaces. In order for the fifth inequality 
to be valid, we have to make sure that (recall the definition (11.37^ 1 



p ft 



(3.229) 



for our choice of the triple (a,/3,7). Note however, that c-l-7 = Z — 1 — (p — 7) and since a + /3 + 27 < 2p, we 
immediately conclude that a + 13 +1 <2{p--f) + 1. This implies that ||Vm™:^;^||^„+3+i(j^„^ < CDpLt^^ Since 
0<7<P, the inequality (|3.229p follows. Thus, from (|3.228n we infer 



lMl=2p 



"T" 



^ss F^ 



(3.230) 



Recall that = X]a+/3+27<2p+2-^">.3^7(''''")^?^n'"t'^+^- virtue of the trace inequality, it is easy to see from 
the previous formula that 



The proof of this estimate is analogous to the proof of (|3.229p . Therefore 



JO Jc Jo V Jo q=o-'° ^ •'^ ■^VJo 



(3.231) 



(3.232) 

where we used the estimate p.23ip in the second inequality and the trace inequality combined with the 
interpolation inequality between Sobolev spaces. Using the Cauchy-Schwarz inequality, the estimate p.23ip 
and the trace inequality, it is easy to obtain (recall that F^ is defined in the line after (|3.214p '): 



ft f P ft ft 

/ F,F, <XJ2 ^r' + T/ 

Jo Jc g=0-^0 



(3.233) 



Estimating term J^jJi^Ip- Recall the definition p.222p of /p and (|3.212p and (|3.214p . Note that for C = 917, 
— because in any of the products of the form amCj, bmdj and bmSj at least one term has the form 
9s"9„Uj6 and is thus equal to on dfl. Let us now assume that C = F™. Let < fc < 2p be an even number 
and fix 0<m,j< fc/2 such that m or j is strictly smaller that k/2. Assume without loss of generality that 
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j<k/2. Let us denote a:—p~ k/2 + m and f3:^p—k/2+j. Note that a<p and /3<p— 1. With these 
notations the expressions for a„i and Cj take the form a„i = C9s2c,+2u^_i_„ and Cj — CdgifidnU^i-i-fij where 



C is some positive constant. We can thus estimate 



/o Jv^ 



ds2a+-2U^i_i_^dg-2f)dnU^i_i_fi 



JT" 



< C\ |9s2° + lU^-l-c I lifl/2 I |l9s2/3a„U^_i_3 I 



<A / 2) 



-'4 



^2<, + i(f^,„)-ry 



C 



t!-l~/3||^23 + l(Q™) 



(3.234) 

In the first inequahty we used the H^^^ version of the Cauchy-Schwarz inequahty. In the last inequahty we 
made use of the fact that /3<p — 1. In the case that k is odd, in analogous manner we can conclude 



"T" 



^ 9=0-^0 



(3.235) 



for a pair of indices 0<mj'< such that m or j is strictly smaller than As to the expressions of 
the form Jq JY^mbmSj, we have to allow for the possibility that both m and j may be equal to [|J. Let us 
abbreviate a:—p— [^\ +m, and P:—p— [|j +j. Note that a,(3<p. With these notations the expressions 
for bm and Cj take the form bm = Cdg2a,+iu^i_^_^ , and Cj — C9^2(fi-i)+i9„it^_i_(^_i) , where C is some positive 
constant. Hence, using the same idea as in p.234p . we deduce 



/o Jr" 



Jo ^ Jo Jo „^nJO 



p-1 rt 



Here, we make use of the fact that f3<p which implies /3—l<p—l. Combining (|3.234[) - (|3.236p . we conclude 
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We now proceed to estimate the second, third and the fourth term on RHS of p.223p . 

Estimating term J^ds2p+i[y,dt<:]u"^'^^dg2p+i a„u^+^ From the definition (PIBD]) of [7,9tc], we realize 
that if C = dfl then ds2p+i [j ,dt<:]u'^'^^ ds2p+idnU^~^^ = 0. Hence, the expression we want to estimate is given 
by: 
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Let us fix some < d < c— 1 = I —p — 2. Then 
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(3.238) 



where we used the i/^/^-version of the Cauchy-Schwarz inequality. Note that, by the trace inequality, we 
can estimate 
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where the trace and the interpolation inequahties are used respectively. By taking L'^ and L°° norms 
respectively and the trace inequality, we obtain (recall (|1.34[) ) 

Ja 0<s<t q^i-'o 

After writing d = l — 1 — {I — I — d) and using the definition (|1.37p of Dq, we need to establish 2p+l + q + r< 
2(i — 1 — d) + 1 in order to bound the first term on the right-most side of (|3.239p by AJ)J^'|i^. Since q + r — c~d 
and c = /— p— litisa simple calculation to see that 2p-|-l-|-g-|-r<2(^ — 1 — (i) + lis equivalent to d < c, which 
obviously holds. Using the inequalities between p.238p and (|3.239p . we conclude 
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In a fully analogous manner we can prove 
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Evaluating j^ds2p+idiaudg2p+idlaUn via boundary conditions. Finally, we want to evaluate the first 
integral on RHS of p.223p using the regularization (|3.85p It is important to note that the pull-back of the 
operator ds (on F™, denoted by 9r™) takes the form Sr™ = J^^^s on the sphere Then the first integral 
on RHS of (p:^ takes the form (recaU i/'" = (.g",i?™) where H is given by (0): 
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Let h:§^ ^M. be a smooth function and let the operator dh be defined by dh — hdg ■ In addition to this let 
L be another smooth function, L : ^ M. It is easy to see that 
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(3.243) 



where Pk-q{h) — '^j^c^hg^i . . .hg^^k , where the summation runs over some multi-indices v — (vi, . . . jVk), sat- 
isfying ^'^^ii'i — k — q, i'i>0. In this notation, we realize that 9r™ — ^|g™|-i ^^^d we will make use of the 
formula p.243p with h=\g"^\-\ Recall the formula H"' = -f"- - ^^^^^ + N* {f""). Hence, using Leibniz' 
rule, we get 
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where is defined by dHU]). Hence from p.244p and (|3.243p . with h=\g™\-\ we get 
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where 

^2p-fl,c 

Also note that 
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where is defined by dSH]). Thus, from (|3.247p . (|3.243p (with /i=|.g™|"^) and the Leibniz' rule, we 
obtain 
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Using formulas p.244|) and (|3.247p . we shall first evaluate the first integral on the rightmost side of (|3.242l) : 
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(3.250) 

Now we may apply the identities (|2.78p and (|2.79p to the integral on the right- most side of p.250p . We set 
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f^+if., i^=|ff"|2p+i, ^ = i?", ff^ff", G' = G2p+i^= and /i = /jSp+L'^+i. Thus, from (1^751) 



and (|2.79p . we obtain: 
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where Q is defined by (|2.75p . For the sake of simplicity, let us denote Q2p+ic-~ 
Q(/e"p+itc,/e'?pVitc,^",|ff'"P^+^)- Recah the notation :^ fga^b for any a,beE. Before we proceed 
with the estimates for /q /§i Q™p+i c state the following inequality that will be repeatedly used. Namely 
for any pair {p,c) = (pj — 1 — p) , we have 
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The estimate p.252p follows by applying the operator d'^'^ to the regularized jump equation (|3.85p . We then 
take squares, integrate and use the trace inequality. We can further refine the estimate (I3.252P using the 
interpolation on the term (recall (|1.35p ) to obtain 
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(3.253) 



Estimating /o/giQ2p+ic- The first four terms on RHS of (|2.75p are the cross-terms. The generic 
constant used in the estimates for these terms will be denoted by C to indicate that they are cross- 
terms. In the case of the first two cross-terms, we only analyze the case when c>0, because the term 
/o /r™ c>s2i-iu'"+^9g2i-iy"+^ is estimated separately. For the first term we have 
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(3.254) 

As to the second term on RHS of (|2.75p . we have 
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(3.255) 

where we argued analogously to the estimate p.254p and in addition to that, in the last estimate we used 
the inequality (|3.253p with ri = e^. The third and the fourth cross-terms are estimated in a fully analogous 
way, making use of (j3.253p (with 77 = e^) and the trace inequality: 
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(3.256) 

As to the fifth and the sixth term on RHS of (|2.75p . we first sum them and rewrite in a convenient way: 
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Notice that by the Sobolev inequality and the boundedness of ||/™||ij2 
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Then, keeping in mind that c>0, we easily estimate the second integral on RHS of p.257p : 
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where we recall that the bound on the first moments of dc+if"^'^^ : | |5c+iPi/'""'""^| <C£>™"'"^ holds due 
to (I3A391) . As to the first integral on RHS of KWll . note that 92p+i(/™ + /™) = gsp+i^p^^ jm ^^p^^ jm) 
(recall p.20p ). Thus, we may write 
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Just like in (I3.259p . we easily bound the second integral on RHS of ()3.260|) : 
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The term - is estimated analogously to (|3.258p . As to the first integral on RHS of (|3.260p . 

we first extract a full time derivative and a cross-term: 
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(3.262) 

Observe that the first term on the right-most side of p.262p gives exactly the last energy contribution in the 
definition (jl.35p of <£g with q = p. Furthermore, 
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where we exploited the assumption c> in the estimate \ la^P+^Pa+J^+i] |2 ^ < C /g The third term 

on the right-most side of (|3.262p is a cross-term and the estimates are analogous to p.255p : 
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The fourth term on the right-most side of p.262p is estimated analogously to p.259p : 



(3.264) 
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(3.265) 



The last term on the right-most side of (|3.262p is easily estimated after a routine application of integration 
by parts so that the number of 9 derivatives hitting 9cF2+/™ is exactly 2p-|-2. We obtain 
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where the estimate analogous to p.258p was used to get the bound || |^„|4p+j — l||^i,oo <C||/™||//3. Sum- 
ming (|3.259p - (|3.266p and using (|3.257p . we conclude the bound on the sum of the fifth and the sixth term 
in the expression for /q /gi Q2p+i,c- 

^ .t I .t .t (3.267) 
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In order to estimate the expression Jp„ 5s2p+i9f*cu"'^^c)s2p+i9(* F'"^^ in the case c = 0, we shaU not 
evaluate it using the boundary conditions for tt'"+^. Instead, we proceed in the following way (note that 
c = imphes p = l — l): 



(3.268) 

where the trace inequality and the interpolation are used in the first inequahty (note that the potentially 
large constant ^ stands in front of the temporal dissipation term) . We now estimate the seventh and the 

eighth term in the expression for /q /gi (52f)+i c (where, recall, Q is given by ()2.75p ). Here, however, the 
possibility that c = is permitted and in some instances we have to proceed differently in our estimates if 
c — 0. In order to estimate the seventh and the eighth term in the expression for Jgi QTp+i c (where, recall, 
Q is given by (|2.75p ). we employ similar strategy as in the estimates leading to p.267p . In analogy to p.257p 
we have: 
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Upon integrating by parts twice, the second term on RHS of (|3.269p is estimated as follows: 
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(3.270) 

To estimate the first integral on RHS of (|3.269p . in analogy to p.270p . we first note that 
<Ce sup \\n\U\d'/+'P2+n\H^+Cte /Vc"+Ve2+'llL < C sup \\n\m^p+Ct Ts^+i 

0<s<t Jo 0<s<t Jo 

(3.271) 
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As to the first integral on RHS of p.269p . in analogy to (|3.262p . we note that 
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Observe that the first term on RHS of p.272p gives exactly the last energy contribution in the definition (|2.47p 
of 2:^. Routine estimates give the following bound: 
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(3.273) 

Note that we made use of the estimate p.252p (with 77==e), to get the bound e /g ||5^+tVe"3+^||i2 < p£"+^ + 
f^D^'p"^'^. for anyp = 0,...,/-l. From ([3:2691) - KTm . we finally obtain 
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(3.274) 

The last two terms in the expression for /q /gi Q^p+i c involve the expressions G^^"*"^'^ and defined 
by (|3.246p and p.249p respectively. The following inequality holds 
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(3.275) 

where in the case c = 0, we only estimate the e-dependent terms on LHS. These estimates follow by a routine 
use of and L°" norms, together with the Cauchy-Schwarz inequality. 

Reasoning in the way similar to the derivation of estimates (|3.230p . (|3.232|) . (|3.233p . (|3.237p . (|3.240p 
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and (|3.24ip . we can deduce the bound on the last seven terms on RHS of p.227p : 
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Estimating RHS of (|3.227p . We now extract the leading order contribution from the first term on RHS 
of (|3.227p . Note that 
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We now apply the identity (|2.8ip to evaluate the integral on the right-most side of p.277p . We set x = feipf ' 



f^tl, -F=|ff"pP, ^/' = i?'", 3 = 5™, G = G2P^'= and /i = /i2p.c+i. Thus, using (jXWTIl . we obtain 
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(3.278) 

where 5' is defined by SUm . For the sake of simplicity, let us denote ^^^^ — '5'(/e™ptc,/e'5p+o,i?'", |5"Tp+^)- 
Estimating /q /gi c- The first four terms in the expression for /g /gi S"^ ^ are cross-terms (given by 
the first four terms on RHS of (|2.76p ). We proceed analogously to the estimates (I3.19ip - p.l94p : 
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where we used the inequality (|3.253p in the last estimate. The second, third and the fourth term in the 
expression for Jg^ ^ are estimated in fully analogous manner. We obtain 
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As to the fifth term, we have 



10 Jsi 
rt 



Ct 



jm.+ l 



Jgi 



Dm 



< 



i|gm|4p 
im+l 



(3.280) 



(3.281) 



where we used the estimate analogous to p.258l) to bound 1 1 | gm | 4p ~ 1||l=°- In a fully analogous manner we 
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estimate the sixth, seventh and the eighth term in the expression for /q /gi ^, where S is given by (|2.76p : 
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The last two terms in the expression for /o /gi5^,c involve the quantities G^P+i'^ and /i^p.^+i^ defined 
through p.246p and (|3.249p respectively. Similarly to (|3.275p . the following inequality holds 



f 



f 
J t 



Q2p + 2fa 



10 JSi 



|(^m|2p /j>m|gm|2p+l ^ 



, rm+l r/m fm+1 



|gm|2p4 



15" pP 



^2p,c+l(-f2p+l,c+l 



Jsi 



<c[VW^ + \\r\\w){ s)™)+c(%/e^+||niHO( / 2) 



(3.283) 



3.2.6 Proof of Lemma [372] 

Combining the estimates (13:2321) . (13:2381) . (13:2411) . (13:2541) . (13:2551) . (13:2561) . (13:2631) - (13:2641) . (13:2791) 
- p.283p together with the identity (I3.203|) , we conclude for any < p < / - 1 , < t < T < 1 



/ ®rp+'<*S.,p(0) + C77 / + — / S);,"+i + C(A + \/€^) V / + — + 

Jo ' Jo V Jo <;=p+l"'0 











Jo 




Jo 





+tE ['k^'+c sup ||/™||ff4(€™ + e^'+i)+c(t+i)%/e^(€™ + £: 



rn+l\ 
P ) 



A 







A 



01 +— e;"+i+(7< sup ||r"||H^[e;;" + g™+i] 



A X-qe ' o<s<i 



CCAt^ 



+c sup iini^.(-(g™ , . . , , 

o<s<t e Jo ^^Jo e o<s<t 

+Ce sup + 

0<s<t 

(3.284) 

where we employ the convention X^a'" — if a,6eN and a>5 and :=0. We choose A, r] small, and 

allow the constants C to become larger (depending on the choice of A and rj). Further, recall the smallness 
assumption supQ<g<j (£(s) + £(s)) + [T>{s)+D{s)] ds < 02 introduced in the paragraph after (|3.197p . We 
conclude that there exist constants Ci, C2 and C2 such that for any Q<p<l — 1: 



^l_Cit-^) / I)™+i) + (l-(72i-C2t)(£: 

^ Jo 



m+1 
P 



1)^+1) 



q=p+l 



q=0 



cxt 



<£e,p(0) + z^' V / D"+i + Cy / + ^ sup £^+^ + C* / P^+i, 

e o<s<t Jo 



(3.285) 



whereby > is a small, but fixed number and 



A A 



+ sup ||/''"||ff2(S™ 

0<s<t 



-9^e, + CU£l-,^+Ce sup ||nil,4f" + Ce||Pinii.. 
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Furthermore, let Di ■.= 1-Cit-^, L»2 := 1 - (^2^ - C2i. For each 0<p</-l, muhiply (|3.285p by i^p and 
sum over all p = 0, 1 , . . . , ^ — 1 . We obtain 

p 



p=0 p=0 p=0 p=0 p=0 

Z-1 i-l „t i-l p-1 „f i-l ,/V\-i /-t 

+^'y,vP V / S)™+i+cVi.fV / S)™+i+Vz.p(^ sup + / 



From here we infer 

i-i i-i „t i-i p-i i-i 



-iin+l 



p=0 g=l "^^ p=l 9=0 q=p+l ° 

+i?2^i^^ r2);^+i)<^i^^'ee,p(o)+^i^p^^+^i^p(^ sup 

n -^0 n . n n ^ 0<s<t Jo 



p=0 p=0 p=0 p=0 

Let us set rfo = ^i ^ X]q=i '^'^ for 1 < p < Z — 1 set dp — Div^ — v'' X]q=o ^"^ ^ Sg=p+i • From the previous 
inequality we immediately obtain 

^0 Jo Jo ^ ^ ^ 0<s<t Jo 

(3.286) 

where (Jq :=X]p=o'^^'^e,p(0) ^^id i/"^ :=X]p=o'^^'^''- -'^'^^ 0<Q!<1 be a real number to be specified later. 
Multiplying the inequality p.286p by a and adding it to p.l99p . we arrive at 

^m+i^ f d^+^^Jd^ sup €"+i^'' + ydp / S)™+i)+ai?2( sup / s^+i.-) 

Jo ^ 0<s<i Jo ^ 0<s<t Jo 

<£,(O)+^(02,<) + ^(f™+^ + (r"+i)2) + (^ + -^)£"+i + ^£™+i+C(Af?2 + y^) Tl^r' 

+ (CA+^ + ^ + CT^/') rS)™+i+C02 sup (£™+i+af<5- + J- + y,.pf^ sup Tp^+i 
A eA Jo o<s<t ^ ^ ^ e o<s<t Jo 

(3.287) 

Absorbing the factors of f /oX>^+^ and /o^r^^ from RHS into LHS and recaUing that + 

/o*S''<2 + /o*D<i^"('"^^(e'' + /o*S''), we obtain 



i-i i-i 

-im+l 



(i_^_^_^_^y^p)£™+i+(i_c(A+02+v^)-«c*y^'') / i?r 

A 77A eyyA 7o 

+ (ai?i-C02^-('-^)) sup e-^+^^'^ + yfadp-CX-^-^-Cel^'-^) (3-288) 
o<s<t ^ ^ A eA A / Jo 

+aZ?2( sup /*i)™+i''')<f,(0)+/C(^2,i)+a(5^+aJ" + ^J^"'''"'^(£"+''")' 

0<s<t JO ^ 

Now set a = min(^, =^7—, — ,1) and choose A, u, 9i, 60, v and then t — f^ so that 

' 20c* X]p=o'^'' : r-j II 
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Let 6q = 02V^ ^ and assume £^{0) + a5'^ <9[)/2 (and hence £^e(0) + q;(5o < ^'2/2), we obtain from ^ 
sup / sup /" 3™+!^"^) < 16*2 + ^(0 sup (S^+i''')^ 

0<s<t Jo 0<s<t Jo 9 9 o<s<t 



From the previous inequahty we conclude that there exists a (possibly smaller) such that the first claim 
of (|3.f f 5p holds. This finishes the proof of Lemma [ 



3.3 Regularized Stefan problem 

The principal goal of this subsection is the following local existence theorem: 

Theorem 3.4 Let C>0. For any sufficiently small L>0 there exists t'^>0 depending on L and e and 
constants < q;,i^< 1 such that if 

then there exists a unique solution (u',/') to the regularized Stefan problem defined on the time interval 
[0,i'^[, satisfying the bound 

sup E^^,{u\f'){t)+ f D^,,{u\f'){T)dT<L 

0<t<t' Jo 

and Ea.^{u'^ ,f'^){-) is continuous on [0,t'^[. 

3.3.1 Local existence for the regularized Stefan problem on [0,t'^[ 

Having obtained the uniform bounds on the sequence {u^ ,R"^), our goal in this section is to prove that this 
sequence is actually Cauchy in the energy space. In order to do that, we need to compare two consecutive 
iterates in the energy norm. However, the expression ^^(u™^^ — — i?™), does not a-priori make 

sense, because the functions u™"*"^ and are defined on different domains il™ and respectively. 
To overcome this difficulty, we shall pull-back the functions u*' from Of^^ onto f75j^(3.g gj, where we recall 
Q,s^[xQfl)=^\Si{xQ^Q). To apply this change of variables, we first define the map x = 0™(i,x), where 
0™(f,x) = (a;o,0) +7r™(i,x)(x— (xo,0)) and tt™ is a smooth scalar-valued function with the property 

^ ' \l, |x-(xo,0)|-l>2d ^ ' 

By abuse of notation, we shall interchangeably use the fact that R"^^^{ |xl(xo'o)| ) — for a unique 6 

such that (cos6', sin6') = |^:[;^°'°]| . For any x e r™-i , we have x = (xq , 0) -f i?™- ^ (9) ^^Z[ll'l]i , by the definition 
of r™-i. Hence 



X 



(a^o,0) x-(a;o,0) . , 1 nm-i./jx ^~ (^o,0) 



e™(i,x) = (a;o,O) + 7r"(t,i?'»-^(0)| p-^)i?"-^(^)- pL^ = (a;o,0) -— ^i?™-^(0) 

|x-(a;o,0)| |x-(a;o,0)| i?™ ^[9) |x-(a;o,0)| 

-K,0) + ^^^e5i(xo,0). 
|x-(a;o,0)| 

Thus 8™ does the job for us: 8™ : il™^^ — > r25j^(3;g o)- Note that this map is natural from the point of view 
of our geometric assumption that the evolving interfaces F™^^ are, in fact, 'graphs' over the unit sphere 
S'i(a;o,0) centered at (2:0,0). For notational simplicity, from now on we shall only work with the case xo = 0, 
i.e. S'i(a;o,0) = S^ and we shall also drop the boldface notation and denote the vectors x and x by a; and x 
respectively. 

The inverse map x{x) to the above change of variables (|3.289p is given by 

x = p"^(x)x. 
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To convince ourselves that the function : flgi is well defined, we observe that has to satisfy the 
relation n™{t,p'^{t,x)x)p^{x) — 1= 0. The existence of such a p™ can be established by the implicit function 
theorem, applied to the equation F(t,s,a;) = 0, where F:MxRxK"^Kis given by F{t,s,v) =TT"^{t,sv)s—l. 
Namely, it,s,v) = (Vtt™ • w)s + tt" > since tt" w 1 and Vtt'" = 0{f^^^) is small. We define : f^gi ^ M 
by setting 

u"\x) ■.= u"\t,p"'{t,x)x). (3.290) 

The heat operator dt — A on the domain fJ™^^ will transform into a more complicated operator in the new 
coordinates. The following lemma describes this description in the new coordinates: 

Lemma 3.5 The push forward of the operator dt — A with respect to the map 0™ : Sl™^^ — + §^ reads: 

{dt-A)*u^ut-AsU- a'^jUs.^s,^ - , 

where 

a™ := ((7r")2 - 1)*^ + 27r"a;%^ IVtt^P, b^' ■= 27r™ + Att^x'' - ttTx'. (3.291) 

Furthermore, 

= (Kr+']t ) ° r (3.292) 

Proof. Note that (I3.289|) implies (e; is the i-th unit vector, i — 1,2) 

x^^^TT"'e, + TT^,x. (3.293) 

From tt"^{x) — l/p™(x), we obtain tt™ = — (Vp™-a;^i)/(p'")^, which in turn, combined with (|3.293[) . after an 
elementary calculation implies 

7rg = - , (3.294) 

where ^t^™ = p"* — Vp™ • x. After further differentiating (|3.294p with repspect to x^ and using the rela- 
tion (|3.293p . we arrive at 

^ -Ap™ 2|Vp™|25'™-2|Vp'"|2Vp'"-x + p'"|Vp™|2(*™)-ixJxV^._j, 

Similarly, it is straightforward to see that ttJ" = — p™/(p™)^. In order to evaluate the Laplacian in new 
coordinates, by p.290p we first write 

u'"(t,a;)=u'"(i,7r"'(i,a;)a;). (3.295) 
Applying A^ to the left-hand side, we obtain 

Au™(t,x) = (7r")2A2w"+(27r™x^'7r™4-a;*x^'|V7r™|2)t2™_^ + (27r;'|-KA7r™a;*)M^ (3.296) 
Applying dt to the left hand side of l|3.295p . we obtain 

m m 

Ut {t,x) = Ut -l-TTj Vxw -a: = 7/4 - ■^^^;;r)2 ■^ = "t -^^JiU^,. (3.297) 

Since {dt - A)u™ = 0, we conclude from ((3:296)) and ([3:297)) 

(dt-As)u =a,^u^^^^+b^ 
where and &^ are given by ()3.29ip . Next, using ()3.294p . it is easy to see that 

% (p™)2 ^(-pm)2^m + (■pm-)2(-^m)2- 
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We now turn to the proof of (|3.292p . Note that 

m+l 

From here, we infer the formula 

Denote ^ = (cos0, sin0) and let n and r stand for the unit normal and unit tangent on respectively. From 
the above formula, we obtain 

1 

[Vu-^y_or[e)^-;;:^^[Vu-^r-{0--^^^ (3.298) 

Note that 

since = (due to the fact that m™+1'+|si =u'"+i^"|si) and ^-71=17112 = 1. Observe further that 

p'"+i(^) = i?™(0). Furthermore, since Vp™+i(C) • C = ^ep^+V • ^ = 0, we also have *™+i(^) =p™+i(^) = 
R™'{9). Form these observations and from (|3.298p we obtain the formula 



unit normal on respectively. Hence 



It is straightforward to see that 77r™ o0™ = t^^77— -M^t, where r and n stand for the unit tangent and the 



( [V7i™+i]+-r7r™) o0™= f— [<'+i]+77-^i^[<"+i]+r) • f ^n-^^r) 



Iff"! |ff™|(i?'")2^ " (i?™)2^ " 



and this proves p.292p . □ 
Since (^[Vu'"+^]± ■ rip^n ^ o 0'" = ^-^j^l'i^"' + e ~|-^''„t| , from (I3.292p . we immediately obtain 

j-m+l/nm^S fm+l/'nm\2 fm+lrfm fm+1 . fm+l r>m j-ni+1 , / nm\2 

r-m+ii+^ Jt ) /e^f \^ ) ^ It ^ hH I / /t ^ hn \( \^ ) 



(3.299) 



where 



/•m+l nm rm+l ( rtm\2 I ^"i I 

S(/") = S™ := f ^V-^ + (^T^ - 1) = (1 - + 



15"! 15"! ' ^ {R^Y" " 



For any fceN, let ■.— u^'^^ — u^ and cr''"+^ =/''^^ — /'' • By subtracting two consecutive equations in the 
iteration process, we obtain: 

?;™+i~Aw™+i = /°, on f7™ (3.300) 

m 

i;™+i = -a™ r^ + G° on §\ (3.301) 

?;™+i = on dil, (3.302) 

m+l nm rr"^+^ 

= + C on §\ (3.303) 
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Here 



Gi,=N*{r)~N*{r~'), 

rym ryrii—l i i 
7 r7n f ^ \ _ c f^^^ ( \ 'Z^'ii^ I ':7m— 1 

Here, the quadratic nonlinearity N*{-) is defined in the hne after p.96p . Upon applying the differential 
operator dtk to the equations p.300p . (|3.301|) and (|3.303p and singling out the leading order terms, we arrive 
at: 

v^4l-^v^^'^fL (3.304) 
vT' = -<^-J^\+G'^ §^ (3-305) 

[dnvrr-=-<n^^^~e^+K, on §\ (3.306) 

fL^^t^f:^ (3.307) 
fc-i ^ 

= - ^ o-^t, ( I I ) J,-, + G°„, . (3.308) 



where 



c=-E<^^H^)*.-,-^E'^«"^^+^(n^)*-.+^*^^J^- (3-309) 

In order to obtain the energy identities for the problem (|3.304p - p.306p . we may use the identities from 
Section O Respecting the notations from (|2.62p - (|2.66p from Section [21 for any j G N we set 

T^f'^,F=l,G = g'^,h = hr,U^vl+\u^al+\x^ , and V = W . (3.310) 

Note that r = §^ and hence, since the domain is not time-dependent, we have = and the terms and 
V in equations (|2.63p and (|2.65p are simply equal to 0. In other words, the expressions OiJA) and -P(W) 
defined in (|2.73p and (|2.74p are equal to 0. We define the energy: 

i-\ i-i 
:=E^.(«'\'^'), :=E^^("''^')' 

k=0 fe=0 

where and are defined by (|2.67p and (|2.68p respectively. Furthermore, using (|2.7ip and (|3.310p and 
keeping in mind that O = P = 0, we arrive at 

e™+HO+ / d"'+'^{s)ds^ [ [ {g™ + s"}, (3.311) 



JS,^ 

where = Ei=oQ(^™''^"^''^'")' = E1:1o 5(a7,',a™+\i?™). Here, Q and S are defined by 
and (|2.76p respectively. Just like in Section [3.21 for any 0<j^< 1, we introduce the higher order energy: 

eJ>:=€''(t;^CT^'), V^" ■.= t^''{v\(j^), O^-'^ :=S)'^(w^cr^), :=I)'''^(w^cr^'); (3.312) 

here jGN, 0<fc<^^l. Like in the previous sections, we drop the index v if v = l. Before we explain the 
energy estimates, we first examine the form of G"; and featuring RHS of equations (|3.300p . (|3.30ip 
and (I3.303P respectively. Note that f°-^ can be rewritten in the form 
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We observe that a'ij'^^ — a™ can be written as Aij (cr'", cr™) for some smooth function with bounded derivatives 
Aij. Similarly we can write b"^^^~b™ in the form Bi{a™' ^a'^g^aY^)^ for some smooth function Bi. In 
other words /°j takes the form 

Jni^A,j{a ,ag )u^j +aijV,.j + B,{a ,ag ,(Jgg,a^ 

In order to write in a similar form, note that 

^m — 1 1 



= - (1 - tIt^) + rat (tIt^ - 



and hence 



|gm| ■'e''tV|gm| \g-m-l\) ^" (^m )2 ^ ^ ^ n J - V (^™)2 (i?m-l)2. 



(3.313) 

Similar structure is also shared by G°„. The conclusion is that f°^, G°j and have a quadratic structure, 
always containing one term that can be bounded by the norms of v^^^ — v^ and a^^^ =a^ , and the other 
term bounded by the norms of and (and hence a-priori bounded due to Lemma l3.2p . In order to bound 
the first modes of crj^"'^^, we shall use relation p.303p . If we rewrite (|3.303p in the form 

at ^ +ecrg^J =-\g \Z + \g \h^-f ct^ ^ , 

upon applying dfi , multiplying the whole equation by cr^^ti integrating over S"'^ , we obtain 

^":t^YL^-+e\WXl,^\\h<{\\i\9n[Vn^^^^ 



5 
< - 
- 4 



m+l||2 ^ 

(3.314) 



We note that ||(|ff™|[C+^]o<?!)™)t,||2, < ^e^+i + AUVw^+^ll^i^f^). It is further easy to see that 
lff"|/i^;^)t^lli. + ll(/™^r+')t^lli.<C^^Elrio(lk"lli^ + ll<^^'lliO- Thus, using 11311, we obtain the 



bound 

J2{ sup ||Pi<.'+i||i.+ f\\Piai:t^\\h}<^t sup e-^+i + A f\^+' 

t^n 0<s<t Jo 0<s<t Jo 

'-o<s<t Jo -* 



(3.315) 



+ce 

k=o 

Motivated by p.315p . we introduce the following temporal energy quantities: 

e^-:=e-^+|]||Piaf,||i.; d^' := ||Pia^,,J|i.. (3.316) 

k=0 k=0 

Remark. Observe that the energy norm is not equivalent to £^ , because the first modes of are also 
controlled by in other words > C_b| IPicr-' 1 1^2 for some constant Ce- An analogous remark holds for . 

The next lemma is the analogue of Lemma 13.21 where we collect preparatory statements, necessary to 
establish (I3.318p . Note that in the following the constant Ce depends on e and we shall make use of the 
e-regularization (|3.85p . 

Lemma 3.6 (a) For any jCzN, the following conservation law holds: 



^Jw Jgi JSi Jsi 
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(h) The energy norms and are equivalent to E\\+\\S'i(t\\'^j^2 and T)\\ respectively. 

(c) There exist constants K* and < i^* < 1 such that for any j G N the following inequality for space-time 
energies holds: 

sup + f 1}^+^'''' <C,i0 + t){ sup e^-'"' + f <}^'''')+K* f d^+\T)dT. (3.317) 

0<s<t Jo 0<s<t Jo Jo 



(d) For a possibly smaller there exists a C such that for any t<t'- the following inequality for the temporal 
energies holds: 



sup e™+i' 

0<s<t 



/ d"'+\T)dT<C,{9 + t)( sup e"'{t)+ d"'{T)dT). (3.318) 

Jo ^0<s<t Jo ' 



Sketch of the proof. The proof of part (a) is identical to the proof of (I3.97P and follows by simply integrat- 
ing ()3.300p over SI™. The proof of part (b)is analogous to the proof of (j3.114p presented in Subsection l3.2.1l 
It relies on part (a) of the lemma. 

Proof of part (c) . To prove (|3.317p , we proceed in the same way as in Subsection 13.2.51 For any pair 
of indices (/i,c) such that + 2c < 2Z — 2, we apply the differential operator 9^ to the equation (|3.300p . 
multiply with — A9^'w™+^ and integrate over Jlgi. We then carry out the same reduction procedure as in 
Subsection 13.2.51 (which in fact becomes simpler because the boundary is time-independent). The es- 
timates carry over analogously and the only terms where additional care is needed are the formally new 
expressions that involve G°j and h°^. To illustrate this point, note that for instance, the analogue of the 
expression p.242p in this case is simply given by 

/ / / / 52f+i(-a™--^f^)52j'+i(-^V-r — ^-f^) 

io 4 c L „ J- c V Rm\gm\> c \ |^™| j^^j^ ^^^^^^^ 

Jo Jgi Jo JSi 

Observe that the first integral on RHS above will render an energy contribution to e™+^, as it is explained 
in the energy identities from Section [2l However, the other three terms are formally new. Let us prove the 
energy estimate for the second term on RHS of (|3.319p : /q /gi First, we integrate by 

parts and use the Cauchy-Schwarz inequality: 

f f dl^+^-^+'dl^+'hl^ = f I dl^deev"^+'dl^hl^<[\\dl^deev^+'\\l.+ [\\dl^hU\l^. (3.320) 
Jo Jsi Jo Jsi Jo Jo 

In order to estimate | 1^2 (recall that is given in the third line after (|3.303p ). we use the inequal- 
ity (|3.252p (to estimate e/^J, the quadratic structure of /i^ as given in p.313p and the a-priori estimate 
+ /q* < for j = TO, TO + 1 . We obtain 

I \dl^h°J\l-- < Ce™ (S™ + + Cr"5™+i + CD"'e". (3.321) 

On the other hand, using the boundary condition (|3.301D for w^+i and the definition of e^, we obtain 
\\d^Pdeev"'+^\\l2 < f e". Using this, the estimates ([332Tt and ((3320)) . we obtain 



lo Jsi 



<C{-+ f 1?™ + ^ sup sup e^-fC sup €™ / 

e Jo e o<s<t o<s<t o<s<t Jo 

<C(i/,a)(-+0 + ^) sup e" + C7L / 
e e o<s<t Jo 



Using the definition of (given in the second equation after p.303p ) and the the definition of e^ , we easily 
obtain | |9^'''''^G5^| |^2 < ^^t"^ ■ Using the integration by parts, the trace inequality and the previous estimate. 
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ft 

/o Js\ 



we bound the third term on RHS of (|3.319p : 

t t 

<?/ l|5.^^+^G^IIi^+A/ ||a,^^K+^]i||i.(,.)<^ sup e" + Arc,"+^ 
^ Jo Jo ' ' Ae o<s<t Jo 

Similarly we can deal with the last integral on RHS of (|3.319p . Note that we can afford to have constants 
in our estimates that contain e in the denominator. This fact allowed us to use the jump relation (j3.303p 
to estimate the highest order terms above. This is in fact analogous to the proof of uniqueness in Section 5 
from [21J. Having obtained the energy estimates, we can conclude (|3.317p in a way analogous to the derivation 
of (|3.286p starting from (|3.284p . Here we also exploit the fact that crJ (0) v-' (0) = for any jeN. 
Proof of part (d). The inequality p.318p follows by estimating RHS of the identity (|3.31ip . using (|3.315p 
and Lemma 13.61 These estimates are, indeed, analogous to the estimates from Subsection 13.2.31 The only 
novelty are the formally new terms dtkf°^, dtkG°^^ and dfkh°^^ appearing in the definitions p.307p . p.308p 
and (|3.309p of //^, and h'^ respectively, but they can be dealt with analogously to the proof of part (c) 
of Lemma 13.61 □ 

Using p.318p together with part (c) of Lemma 13.61 yields a contraction bound on the total energy: 



sup e'"+\s)+ / d"'+\T)dT+ sup e™+i^'''+ / 5" 

0<s<t Jo 0<s<t Jo 

<a( sup e"(t)+ / d"'{T)dT]+C{9 + t)( sup e'"''^*+ / d"""') + K* a( sup e"(t)+ / (i"(T)dr) 
<A*( sup e"(s)+ / d"\T)dT+ sup e'"^'''+ / 0" 

^0<s<t Jo 0<s<t Jo 



^0<s<t Jo ^ ^0<s<t Jo ^ ^0<s<t Jo 

rt ft 



(3.322) 

where A* := k{l + K*) + C{9 + t). Note that A* < 1, if we choose 9 small enough and possibly a smaller so 
that A< i^j^^. Having proved ((33221) . we conclude that the sequence (?;'^,o''^)fegN is a Cauchy sequence 
in the energy space defined for some small enough. Passing to the limit we recover the solution to the 
regularized Stefan problem on the time interval [0,t'^[. 

3.3.2 Uniqueness 

We shall use Gronwall's lemma to prove uniqueness of the solutions to the regularized Stefan problem: 

Lemma 3.7 In the class of functions u,f satisfying supp^^^jj + /g Da,i,{u,f){s)ds<L (where L 

may be chosen smaller if necessary), the solution {u,f) to the regularized Stefan problem is unique. 

Proof. Let us assume that there exists another solution (u,/) satisfying the same initial condition 
(m(x,0),/(6',0)) = {uo(x),fo{9)) and the same smallness bound supg<g<(£ Ea,i,{u,f) + /J D a,v{u, f){s) ds < L. 
We shall formulate the Stefan problem on the fixed domain by pulling both functions u and it onto the 
fixed domain $15^(3,^ 9). Without loss of generality we assume 175^(3,^ 0) = Jlgi . To this end, we shall apply 
the change of variables described in Subsection 13.3.11 We define <d{t,x) ■.— Tr{t,x)x and Q{t,x) ■.= 7r{t,x)x, 
where tt and tt are defined as in (|3.289p by substituting i?™^^ with (1 + /) and (1 + /) respectively. We set 
U:=uoQ and U = uoQ. The corresponding coefficients a,y , hi, dij, hi are defined analogously to a™, 5™ 
(given by (|3.29f p ). Subtracting the two sets of equations (that we get for (C/,/) and ([/,/)), we obtain a 



new problem for (y, /i) :=([/-[/,/ - /) (note that i? = f + / and |.g| = i?^ + /2) : 

Vt-/W^A°, (3.323) 

V = -h--^+B°, (3.324) 

^-^I^^.^^^ + C'" on §\ (3.325) 
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where 

A° = (a^j - aij)Uij + dijVij + {h -bi)Ui + biVi, 
B" = A'-(/)-A-(/>/..(^-5^). 

Here N*{f) is defined in line after (|3.96p and and are defined by replacing /™ by / and / in 

the definition of in the line after p.299p . respectively. We use the notation from (|2.62p - (|2.66p and 

Lemma 12.21 from Section [2] to derive the accompanying energy identities. To this end, we set 

J'^A°,F=l,G = B°,h^C°,U = V,LJ = x = h Siud^^R. 

Same remark as in the line after (|3.310[) applies here; since r = S^ is time-independent, we conclude that 
the expressions 0{U) and P{U) defined in ()2.73p and ()2.74|) are equal to 0. In formal comparison to the 
problem p.300|) - (|3.303p . / takes the role of /™+^, and / the role of f"\ Additionally, the cross-terms 
vanish since Ljj^x = h- We are naturally led to the following energy quantities: 

where £e and are defined by (|2.67[) and (|2.68[) . respectively. In addition to this, we define Q* = Q{h,h,R) 
and analogously 5*. Using the identity (j2.71|) . we obtain 

£*{t)+ f V*{s)ds= f f {V + Vt)A+ f f {Q* + S*}. (3.326) 
Jo Jo Ji\i Jo Js^ 

In addition to the above energies £* and V* , we define 

2)* = I|VF||?,.(^)-)-/ {\het\'-\h,\'}+e [ {\fost\'-\feH?} 

Note that the quantities <B* and J)* are identical to £ and D in the case when 1 = 1. Our goal is to prove 
that for some < a* < 1 the following inequality holds 



a*e*{t)+ f {V*{s)+a*'S*{3)}ds<C f {£*{s) + e*{s)}ds + C{VL + a*) f V*{s)ds, (3.327) 
Jo Jo Jo 



which would enable us to absorb the multiple of J* V* on the right-hand side into the left-hand side and then 
use the Gronwall's inequality to conclude that £* (t) -f €* {t) = for any t > 0. It is essential that the constant 
C in the above estimate does not depend on e so that the smallness bound on L remains independent of e. 
That the identity (|3.327|) indeed holds, follows analogously to the energy estimates from Subsections 13.2.51 
and 13 . 2 . 3l applied to the right-hand side of (|3.326p . Here we strongly exploit the uniform bounds on Ea,u{u,f) 
and Ea,v{u^ f). A major difference from the existence part of the proof is the absence of cross-terms in the 
energy identities (since 0' = % in the notation of Section [2|). In addition to that, we work in a lower order 
energy space (1 = 1) and we can thus use the above uniform estimates to bound the term by C^/L in 

L°°-norm. Knowing that the e-dependence comes only from the estimates of the cross-terms (cf. (|3.254p - 
(I3.256P . (|3.279p . (|3.280p . (|3.192p - (|3.194p ). we conclude that the constants on the right-hand side of (|3.327p 
do nof depend one. Choosing L and a* suitably smaU ()3.327p implies f *(i) -h (t) < C/^ {£*(s) €*(s)} ds, 
which, by Gronwall's inequality leads to £*{t) + £*(t) = 0. This finishes the proof of the uniqueness statement 
of Theorem [H31 □ 
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3.3.3 Continuity for the regularized Stefan problem 

The continuity statement is important because it will allow us to extend the solution beyond the the time 
interval [0,i'^[. 

Lemma 3.8 The total energy j^(u^ f)(t) + /q*-Dq ^{u,f)(s)ds is continuous as a function oft on [0,t'^[. 

Sketch of the proof. Integrating over [s,t] instead of [0,t] while arriving at the identity p.ll3p in Sec- 
tion [XI] we obtain 

£'^'+\t)-rj'+\s)+ j v'j'+^{T)dT^ j A{u'''+\f''\f"'+^)+ j J {o"+p"'}+y J Iq^+s""}. 

" ^ ' ^ (3.328) 

We claim that we may pass to the limit as to ^ oo in each of the terms in (I3.328|) to obtain: 

£,{t)-£,{s) + j'v,{T)dT^ J' A + J' J^[0 + P}+ J' J^^ {q + s}, (3.329) 



where 



i-i 



A=Aiujj), o^J^^M, P^Y.PM.Q^Y.'^^'^i-^^i-^^) s=Y,s{ut>^Jt>'.R)- 



fe=i 



k=l 



k=0 



The passage to the limit is justified for the following reason: we may apply the differential operator dtk to 
the equality (|3.295p to express the derivatives of in terms of derivatives of and /™. This enables 

us to express the norms of d^kU™'^^ on fJ™ in terms of the space-time norms of and /™ on figi and 

S^, respectively. Hence, using the strong convergence of the sequence («■' ,/■'), we may pass to the limit to 
conclude (|3.329p . In order to bound RHS of (|3.329|) . by the same argument, we pass to the limit as to— >oo 
in the estimate p.l99p . where we keep in mind that we work on the time interval [s,t]. Note that all the 
estimates coming from the cross-terms will vanish, and those are precisely the ones involving the constant 
C. As a consequence we obtain 



A + 



O + P^ 



Q + S 



<CCo{t-s) sup (£(r) + €(T) + e'i+ / S) 



S<T<t 



+C{\+ sup \\¥if\\l,+cV€) / I?, + C(Vs + A+(£(o))i/2) / S), 

0<T<t Js Jo 

Using (|3.329p . p.330p and the triangle inequality, we obtain for any < s < i < i*^: 
\£,{t)-£,{s)\<C f V,{T)dT (l + C{X+ sup ||Pi/||2,+Ve)) 

Js ^ 0<T<t ' 

+CCo{t-s) sup {£{T) + €iT) + 9i+ [ D) + C(V2 + A+(5(o))i/2) / S)e 



(3.330) 



>0 as s-^t. 



since supq<^<je f e(T) and (£e are bounded on [0,i'^[. This proves the continuity of £e on the small time interval 
[Ojf^l. In a completely analogous fashion, by passing to the limit in (|3.284p . it follows that (Be is continuous 
on the interval [0,i'[. □ 



3.3.4 Momentum conservation and smallness of the first modes 

Recall that the smallness condition on , which is repeatedly used to close the estimates, is guaranteed 

through a t and e-dependent estimate (|3.136p . It is thus not applicable if we want to prove that the solution 
exists on a uniform time interval independent of e, nor is it useful for extending the solution globally 
because the bound grows as t becomes large. In order to overcome these difficulties, we shall prove a crucial 
'conservation-of- momentum' law. Let M<L/2, where L is given by Theorem 13.41 and let (u,f) be the 
unique solution of the regularized Stefan problem with the initial conditions satisfying i?a,i/(wo,/o) < M. 
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Lemma 3.9 Let Paix,y) = x + ^^p^ and pb{x,y) — y+^^^. Then, on the time interval of existence of the 
solution to the regularized Stefan problem, the following conservation laws hold: 

dt[upa = dt[ Fa{R,0)dO + e[ {fg,,Rcoae + Rlfg.t— ^cos^ + " (3 331) 



an 



d 



dt f up,^dt [ Ft{R,9)de + e f {foHRsine + RlfeH-j, .^^S'f^p ■ ^^2 )^^' (3-332) 

where Fa and Fi, are defined by 11.25]) and il.26\) respectively. 

Proof. Note that both pa and pi, are harmonic away from the origin r = and enjoy the property 

(^nPa\g^^dnPb\g^ = 0- (3.333) 

Let further ■.= Q\Bs{0). Muhiply (jl.ip by Pa and integrate over fls. Using the Green's second identity, 
harmonicity of Pa and (|3.333|) , we obtain 

ufPa - / dnU+PadS+ / W+9„Pq = 0, 

a+ Jr Jr 



/ M* Pa 4 


- J dnU PadS- J 







UnPa + I udnPa = 0- 

SsiO) JSsiO) 

Summing the above equations and using the jump condition p.3p . we obtain 

dt / Upa^ VpadS+ / UnPa~ / udnPa- 

J Us Jr Jssia) Jss(o) 

Note that dnPa\ss{o) = drPa\ss{o) ^ cos9 - Rl^ . Hence 

udnPadS^ [ u{S,e){coae-Rl^^)Sde^ f uiS,e)Scoa9-Rl [ u{6,e)^^d9. 



Sa(0) JSi -'SI Jsi 

Note that the first integral on right- most side obviously converges to as ^— >0. As to the second integral, 
observe that 

f u{S,d)^dd= f "('^^^)-"(0^0) eosgd^^ / M,(O,O)cos0d0 = O asJ^O. 
Similarly Jg^^^^UnPa—^0 as (5— >0. Finally, passing to the limit as (5^0, we get 

dt f upa^ f VpadS, (3.334) 



where the integral upa is interpreted in the sense of the existing limit as (5 ^ 0, i.e. upa = lim^^o J^^ upa 
(keep in mind that Pa~ 7 at the origin and this singularity is integrable). Upon abbreviating R^^ = {Rcos9 + 
xo)"^ + {Rs'm9y and recalling the definition (|1.25p of Fa, we get 

/ Vpa^ / (-p-p+e— — )(i?cos6' + a;o + i?: ^ )\g\d9 

Jr Jgi l5l l5l Rxo 

= 1 ftdRFa{R,9) + e f {f0HRcos9 + RlfeH^^^^^P^)d9 (3.335) 
^dtf {Fa{R,9)-Fa{l,9)}+e f (/,4,i?cos0 + i?^/,4,-^'°'^ + '''' ^ 



i?2 



)d9. 



and this, together with (|3.334p proves the first claim of the lemma. The second claim follows in a fully 
analogous way. □ 
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Remark. One could conceivably try to obtain a hierarchy of conservation laws for higher momenta, by 
multiplying the equation (jl.lO[) by appropriate functions that restrict to a combination of cosines and sines 
of integer multiples of onS^. For instance, in order to control second momenta, a natural choice for the 
test function p is p{x,y) — y'^ — + R^ ^ ~f since p is harmonic, dnp\dn = and p|r = —i?^cos20 — 
However, this strategy fails to render new conservation laws, since p contains a singularity of order 
which fails to be integrable. 

We shall exploit the conservation laws given by (|1.22|) and Lemma 13.91 to prove that the natural energy 
£(0) (given by ()1.29p ) is in fact positive and controls all the modes of /. To this end, the following lemma is 
necessary: 

Lemma 3.10 There exists a constant C such that the following inequalities hold: 

|ai| + |foi|<c(|mo| + |mJ + |mb|+q|Vu|U2(f,±)+C^(|afc| + |&fc|) + Ce\/i( f WftWl^f'^) (3.336) 

k>2 •'^ 

|ao|<c(|mo| + |ma| + |mfc| + ||w|U2(n±) + ||Vu|U2(o±)+^(|afc| + |6fc|) + et / \\ft\\h). (3.337) 

^ fe>2 ^" ^ 

Proof. Recall that mi — J^uoPi + Jgi Fi{Ro,0) for i — a,b and toq is given by (ll.27p . Observe that 

Fa{i + f,e)-Fa{i,e)de= f dfFa{i,9)de+ f o{f). (3.338) 



§1 



From the definition p.25p of Fa we have dfFa{l,9) = coa6 + xq + Rl (cose+ao)^+'sin9^ • Using the decomposition 
f = J2k=o'^kCOsk6 + bksmk9 and abbreviating Ga{R,0) :— fe4fRcos9 + Rlfgit , we conclude 

f / 2n n n.2 cos^ 6* + COS 6* N ,„ , /" . „ ^ , 6* COS 6* + xq slu 6* ,^ 

oi / ( cos^ 6* + a;o cos 61 + .^^dO + hi sm6{cos6 + xo) + Rl—— ^—^d0 

Jgi l + 2xocos0 + a;5' Jgi l + 2a;ocosy + a;Q 

= / Fail + f,e)-Fail,e)dO- [ 0{f)- [ dfFailMao + y2akCOske + bkSmke) + e [ Ga{R,e)d9. 

(3.339) 

Observe that the second integral on LHS above vanishes since smecose+xosme j function. On the 

other hand, the following claim holds: 

Claim: 

I (cos^e + xoCose + Rl^^^;^±^^^^)de^{l + Rl)7r. (3.340) 
Jsi *l + 2a;ocos6l + a;2'' v *^ v ; 

Proof of the claim: Obviously Jgi (cos^^? + xocos0)(i^? = 7r. To prove the claim it suffices to show that 

<'^'^ cos^0 + a;ocos6 



Q l + a;Q + 2xoCOS0 



de^TT. 



We first compute 

2 



1 /■2-{cos0 + ^} - '-^cost 



de=^ I ^ A , ^" '-^de 



l + Xo + 2xoCOS^ 2X0 7o i±^_|_cos0 

1 1,2 ^,2 r2n cos9+^-^ 



_/ {ll^ + cos9}d9-^^ ^d9 

2X0 Jo 2X0 i+£i+cos0 



2x„ 

{l+x2}2^ {1+X2}27r {l+x2}2 /•2- 1 



4xo 2xq 4xq Jo 1 + Xo + 2xoCos( 

4xo 4xq Jo 1 + Xo + 2xocos6I 



-d9 
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By a standard integration formula, 







^ — de = arctan{ ^ tan -}r ^ = ^. 

We note that in the stable vase, |xo| < V2— 1 < 1. We hence deduce 

Jo l + a;^ + 2a;ocos6' 4a;g(l-a;^) l^^^o l~^o 

On the other hand 



9 ndO^ I ^^LZZZ:^ — 2_! — ^—dO 

H-a;o + 2xoCOS0 Jo 1 + 2:0 + 2a;o cost' 

l + xl r^'' 1 ,^ l + xl 27r 2a:g7r 





^TT-^-r^ / -de = TT 

From here, we obtain 



2 Jo l + a;o + 2a;ocos6' " " 2 1 — 2:0 ^^^l' 



cos^e' + xocose' {l + a;^}7r 2x^7r 

9 7. — 9 9 ~ 

l + a;o + 2a;ocos6' I — ^^q 1 — 

and this finishes the proof of the claim. □ 

Using the fact that J^^ Fa{l + f,0) - Fa{l,0)de = ma- J^upa, we deduce from p.339p and (|3.340p 

{1 + Rl)7rai^ma- j upa~ [ 0{p)~ [ 9/i^a(l,e)(ao + Vafccosfc0 + 6fcsinfc0) + e / GaiR,e)d9. 

(3.341) 

In a fully analogous way, we obtain 

{1 + Rl)7rbi:^mb- [ upb~ [ 0{f)- [ 9/Fb(l,6l)(ao + VafcCOsfc6l + 6fcSinfc6i) + e /" Gb{R,e)d0. 

(3.342) 

Note that for i = a,b, we have 

/ 0(/')| + | / a/F,(l,0)(ao + ^afcCosfc0 + &fcSinfc0)|<C(^a2. + &2 + |ao|+^|a^| + |&^|) 

•^^^ "^^^ fc>2 fe>0 k>2 

<CVMi\ai\ + \bi\) + C{l + VM)J2iM + \bk\). 

k^l 

In order to bound the second terms on RHS of p.34ip and p.342p respectively, we use the polar coordinates: 

[ upa = [ [ u{r,e){r+—)coserdrde= [ [ u{r, 6) [r^ + RI) cos 6 drdO 
Jn Ja Jgi ^ Jo Jgi 



= - I I U0{r,e){r'^ +Rl)smedrde = - [ [ {-u.j:rsm9 + Uyrcos0){r'^ + R^^)siii9drde 
Jo Js^ ^ Jo Js^ 

<{Rl + Rl)l I {\u4 + \uy\)rdrde<C\\S/u\\L2(n), 
Jo Js^ 



and analogously for J^upb- Using the integration by parts, boundedness of on [0,T] and the Cauchy- 

Schwarz inequality, it is easy to see that 



/ / GaiR,0)d9 +e f f Gb{R,6)de <CeVt{ [ 
Jo Js^ Jo Js^ Jo 



Using p.341[) . (|3.342p and the previous three inequalities, together with the smallness of M, we conclude 
that 

\ai\ + \b,\<Gi\ma\ + \mb\) + C\\Vu\\L2in±)+C\ao\+Cy^{\ak\ + \bk\) + CeVt{ f WftWh)'^^ . (3.343) 

I.->9 -^0 
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Note further that the conservation-of-mass law (|1.22p impHes 

^— ( f = !I^-— [ -—I 



where we recall (|1.27p . From here and the smallness of Oq, we easily infer 

|ao| < C\mo\+C\\u\\L2^n±)+CY,i4 + bl). (3.344) 

k>l 

Using in we obtain I^MB- The second claim (15357)) follows from ([53^^]) and dSSM]). □ 

We use Lemma [3. 101 to prove a necessary smallness estimate on ||/||//2i. 

Lemma 3.11 There exists a constant C such that the following inequality holds: 

\\f\\H2i <C{ml + ml + ml) + CM + CetM. (3.345) 

Proof. The boundary condition (|3.83p can be rewritten in the following way: 

-f99-.f^uo(f>-N{f), (3.346) 

where A''(/) is the quadratic nonlinearity defined by (|1.9|) . Apply the differential operator to both 

sides of p.346p . multiply with f 921-1 and integrate over to obtain 

{fg2,f-{fe2,-if=[ df-^{uocjy^)fg2^-i+ [ df-\N{f))fe2.-. (3.347) 

From (j3.347p . using the Leibniz rule, the Cauchy-Schwarz inequality, the trace inequality and the smallness 
of II/II//21, we can conclude (see (|1.20l) ') 

\\^2+fe-A\h<C\\u\\ + C\\^if\\h- 

By the Sobolev inequality, we conclude 

||F2+/||?j2,<Ce: + C'||Pi/||2,<Ce: + C£ + C(m^ + m2+mg) + CM||P2+/||i2+Cet f UWl^, 

Jo 

where we used the estimate (j3.336l) in the second inequality. This estimate implies 

\\P2+f\\U<C{ml + ml + ml) + C{<^ + £) + Cet f \\ft\\h. (3.348) 

Jo 

From p.348p . p.336p and (|3.337p . we easily obtain 

11/111^2, < C{ml+ml + ml) + C{€ + £)<C{ml + ml+ml) + C{M/{aj^'-^)+M) + CetM 
< C{ml+ml + ml) + CM + CetM. 

as claimed in the lemma. □ 

Estimates for £(o) revisited. Let M<L/2, where L is given by Theorem 13.41 and let (u,/) be the unique 
solution of the regularized Stefan problem with the initial conditions satisfying Ea{uQ,fo) <M. Now recall 
the definition (fL29)) and the expansion (HH) |.g| - 1 / + /|/2 + ^'(/), where "^{f) = 0{\f\^ +\fe\^). We have 

£(0)i^J) = \l /"+ / {\9\-l) = \ I I n+[ f+ I fl- I ^if) 
^ Jo Jfi Jsi ^Jn Jn Jgi Jgi Jgi 

u) 
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Note that, again due to (|1.22p . 



/)'-(-o-/ y)(/--/ /), 



(3.349) 



and it is hence of third order. Using the estimates p.336p . p.337p and the equahty p.349p . it is easy to see 
that 



2|§ 



The last ineqiiahty imphes that there exists a constant C* such that 



+CVM(ml + ml + ml + | \u\ ^.^^^^ + \\Vu\ + Uo ^ ^f) + j^UWl^) ■ 



^(0)>a(||^/|li.(n±)+ / {fl-Vf})-C.{ml + ml + ml)-C.^I\\^u\\l.^^±)~C.^et \\ft\\ 



§1 



(3.350) 

Thus, for small M this imphes that 5(o) +C*(TOo + "^a + "^6)+'^'(i) + /o^(i) ^ positive definite quantity. 
After summing the zero-th and the first level energy identities and adding C, (to§+to^ + to^) to both sides, 
we obtain 



C.{ml+ml + ml)+£(^){t)+£l^^{t) + \\Vuml2(n^)+ 1 1 Vu| + / + / / {fet^ ft] 

Jo Jo Jo Jsi 

(3.351) 

whereby Oi, Pi, Qi and Si are defined by dropping the index m in the definitions (|3.102p and (|3.103p . 
Estimates analogous to p.l49p and p.l50p imply the bound on the e-dependent terms on RHS of (|3.35ip . 



fe^feH{j^^-l) + ^ I {Nif) + -{\g\-')efg)fgH <Co{\ / V' + \e j V) 



R 



+ ^ sup m\HA^))+C,te sup (||/||?,3||Pi/||i.] 



(3.352) 



3.3.5 Uniform interval of existence: t^>T for small e 

Our objective is to prove that the solution (u'^,/'^) exists on a time interval independent of e. 

Lemma 3.12 There exists aT>0 independent of e such that the solution (u'^,/'^) exists on the time interval 
[0,T] and Ea^u is continuous on [0,r[. 

Proof. As in Subsection 13.3.31 we can pass to the limit as to -^oo in the estimate p.284p . In doing so, we 
use the bound (|3.345p to ensure the smallness of H/^Hh'- Namely, 



ll/"llH'<ll/ll/f' + lir"-/llif'<2ll/llff'<C(TO^ + < + "i6) + CA/ + CeiM, 



(3.353) 



for m large enough. Note that all the estimates involving the constant C simply vanish, and for any 
< p < ? — 1 , we recover: 

sup €p+ / Sp + (1-C2i)( sup / 'hp)<5p + v^ / Dp+i+C^ / Dq + eJi + C* / V, (3.354) 

0<s<t Jo 0<s<t Jo Jo „— nJO Jo 



9=0- 



where Ji := Ctf ™ +Csupo<3<t 1 1/™| 1^4^ " + C| IPi/"] l^^. Multiplying (1?:^ by and summing over all 
p = 0,l,...,Z— 1, just like in the derivation of l|3.286p . we can infer 

sup C' + ^dp / Dp + (1-C2t)( sup e'^+ / S^'")<€:(0) + e^i^'?J'i+^i.«C* / (3.355) 

0<s<t -'O 0<s<t Jo 9=0 -'O 
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whereby we observe that none of the constants in the inequaUty (j3.355p contains e in the denominator. 
Secondly, we pass to the hmit as m—^oo in the estimate (I3.199p . Thereby we use the identity (|3.35ip 
together with the energy estimate (|3.352p . Note that aU the terms containing the constant C vanish and we 
arrive at: 

sup £{s) + C^{ml+ml + ml)+ sup £'{s)+ / V,<£,{0) + C^{ml+ml + ml) + sup £fo){s) 

0<s<t Q<s<t Jo A o<s<t 

+Cote sup (||/||1^3||Pi/||i2)+C(A+||/||H4 +\/m + M) / V,+C{X+\\f\\H^ +Vm) f D 

0<s<t Jo Jo 

(3.356) 

Multiplying p.355p by a and summing it with (|3.356p . we obtain 

i-i pt 

sup £{s) + C\{ml + ml + ml)+a sup + (1 - C(A+ 1 1/| |^f4 +\/M + M) - a V / V 

0<s<t 0<s<t ^^jj Jo 

+ V(adp-C(A+||/||ff4+x/M)) / 2)p + (l-Cot) sup £'{s) + {l^C{\+\\f\\H^ +Ve + 9)) / 

p^O •'^ 0<s<t Jo 

+a{l-C2t){ sup f i)''''')<£,{0) + C*{ml+ml + ml) + a(B'',{Q) + eJ2, 

Q<s<t Jq 

(3.357) 

where J2 « E,=o '^'^i + ^ot supo<s<t (||/||^3||Pi/|||2)- Using the bound (|3.353p to ensure the smallness 
of ^^'^ choosing A and M sufficiently small, we obtain the following energy bound for any t<t'^ 

from d^Wn, : 

sup £{s) + C4ml + ml + ml) + a sup (B" (s) + - eCT M) { [ V + a [ S)^) + (l-Coi) sup ^'(s) + 7 / V' 

0<S<t 0<S<t 4 Jq Jq 0<S<t 4Jq 

+a(l-C2<)( sup / i)'^'')<£{0) + C,{ml + ml + ml)+£'{0) + a€''iO) + ae^''{0) + eJ2. 

0<s<t Jo 

(3.358) 

Now extend [0,i^[ to a maximal interval [0,i^[ such that i?Q,^(s) + j^Docw is continuous on [0,P[ and the 
inequality (j3.358p holds for any < t < P. Let L<9, where 9 is given by Lemma [3T2l Set 

T:=sup|0<t<F: sup Ec,,.{u\f){s) + C,{ml + ml + ml)+ f D„,,(uS/^)(r)dT < l|. 

t 0<s<t Jo 

Note that T>t'^>0. Furthermore, let us choose eo,9>0,T small enough such that CoT,C2T<|, and 
eoCTM<i. Note that £:,(0) + C,(mg + m2+TO2)_^Q,g>'(0)<L/2. From p.358p . we conclude that for any 
t<T and any e < Co the following inequality holds: 



/■* 4 

Jo 5 



(3.359) 



However, if < T, then T <T and the inequality (|3.359p holds for any t<T. This contradicts the definition 
of T together with the continuity of 5e and hence t'^>T. □ 



3.4 Proofs of Stability Theorems 
Proof of Theorem 11.11 

Since the inequality ()3.358p holds on a time interval [0,r], where T does not depend on e, we may pass to 
the limit as e^O. The sequence (w^,/*^) converges to the solution (u,/) of the original Stefan problem. In 
addition to that eCTM^O and eJ2^0 as e^O, since J2 is bounded. Due to the choice of (uo'/o)' have 
£'(0), £'^'"(0) — > 0, as e— >0. Together with the weak lower semicontinuity we obtain the energy bound: 



sup (£{s) + C,{ml + ml + ml) + ae''{s))+l( f V{T)dT + af T)''{T)dT 



0<S<t 

< £-(0) + {ml + + ) + (0) , 



(3.360) 
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for any 0<t<T, where T — min( 2^)- The uniqueness proof is analogous to the proof of uniqueness in 
Theorem l3.4[ where we formally set e = 0. Given the initial data (uqi/o) such that the smallness assumptions 
of Theorem 11.11 are satisfied, we can construct a solution on time interval [0,T[. From (|3.360p . at t — T/2, 
we can conclude f (r/2) + (m^ +mf ) +a(£(r/2) < £"(0) +C*(mg + 771^ +mg) +a£(0). The idea is to solve 
the Stefan problem with new initial data {u^{x,0),f^{9,0) = {u{x,-j),f{9,^)). The problem allows for a 
unique solution that exists at least on a time interval of length [0,r[. Due to the uniqueness we have 
{u^ {x ,t) , {9 ,t)) = {u{x ,t + -j) , f {9 ,t + ^)) and we have thus extended the solution {u,f) to the interval 
[0,^[. Iterating this procedure, we conclude that the solution exists for all t>0. Furthermore we conclude 
that (|3.360p holds for any t>0. In order to prove (|1.38p . we apply the same idea. Let first M* < ^ and 
assume -E'a,;y(uo,/o) + C*(mQ + m^ + m^) < A/*. By the above considerations there exists a global solution 
(m,/) with the initial value (uo,/o) and it satisfies the global bound supQ<g<oQ£'a,i/(s) + C*(r7iQ + TO^ + TOj) + 
/q°° Da.u{T)dT < Fix any s > 0. We solve the Stefan problem with the new initial data (u^(a::,0),/^(6',0)) = 
{u(x,s), f{9,s)). The problem allows for the unique solution since 

M 

Ec,,,{ulJ^) + C4ml + ml + ml) = Ea^,{uJ){s) + C4ml + ml + ml)<—. 

In addition to this we have the global bound supQ^^^^Ea,i,{u^,f^){t) + C^:{mQ + rri^ + ml) + 
Da,^{u^ ,f^){T)dT < M (again by the same consideration as above). We are thus in the uniqueness 
regime and we conclude {u^,f^){t) = {u,f){t + s) for any t>0. We may now use the estimate p.360p to 
obtain pTSS)) . 



Proof of Theorem [TT2l 

Our first objective is to find the candidate point {xo + xo,yo) close to (a;o,0) and the limiting radius i?, so 
that the solution {v{t),r{t)) would converge to {1 / R, S ji{xo + XQ,yo)) based on the 'mass'- and 'momentum 
conservation'. To this end, for any pair (xcj/o) define 



f''{rcos9 + xo+xo + Rl-, J^»^^+^o + ^o )dr, 

' ^ {rcos9 + xo + xoj +{rsmO + yo)'' 



J (rcost' + xo + a;o)"'-l-(rsmt' + j/o) 

Choose R (close to 1) to be the solution of J^-^ + ^ /gi R^ = /^(l + ^0) + 5 /gi (1 + /o)^- 
Lemma 3.13 There exists a pair {xo,yo) close to (0,0) such that 

F^'>>y«{R,9)d9^ma+ I Fa{l,9)d9; [ F^"'"" {R,9)d9 = 171^+ [ Fb{l,9)d9. (3.361) 

Proof. The solvability of the above system of equations follows easily from the implicit function theorem. 
We have to check that at the point {xo,yQ) — (0,0) the Jacobian 



J:=9,„( y^^Ff-'"'(i?,0)d0)a,„(y^^F,^-«"°(i?,^)rf^) -d,,[J^^F^"-y°iR,9)d9)dy„[J^^F^«'y«iR,9)d9 
is different from zero. To do this, we first easily check that at (xQ,yo) — (0,0) 

° " J ^ *r^ + xl + 2rxQCOs9 * {r^ +xl + 2rxocos9)^ ^ ^R2 + + 2x0 co 

From here, by a standard computation, we conclude that at {xo^yo) — (0,0) 

/ d.,F:o:^W)d9^ [ {R + Ri§^^^4i§^)d9 

Jsi _ Jsi _ R'^+xi + 2xocos9' 

- RlR r / (R + xo)cos^ . -| 

= 27ri?+— V<^6' + 2arctan( J ^ ^) \ 

2x1 I ^Rsml^-xnsmi' i 



„n sin0 
R^ 



Xq 



2-KR + = 2nR. 
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By a similar calculation, we conclude that at (xo,yo) — (OjO) 

Hence J = Att^R'^ y^O. This finishes the proof of the lemma. □ 
We parametrize F as a 'graph' over Sj^{xo+Xo,yo) - the unit sphere centered at {xo+xo,yo). We set 
r = {Rcosd + xo + xo, Rs'mO + yo) , R=l + f and v = l/R + u. With respect to the new parametrization of F, 
the conservation laws (|1.24p take the form 

dt [ u + dt [ {Rf+^}; dt [ upa + dt [ F^°^y°{R,e) = 0; dt [ upt + dt [ F^°'^° {R,e) = 0. (3.362) 
Jo Jgi ^ Jo Jsi Jo Jsi 

After integrating the above conservation laws in time, we note that our choice of the center (xq +xo,yo) and 
the radius R of the shifted circle, together with the definitions (|1.27p and (|1.28p of and rUb, implies 



upa+ {F;^°-y"{R,e)-F^«^y«ii,e)de}^o; upb+ {F^°'y°{R,e)~F^"'y"ii,9)de}^0; 

Jgi Jo Jgi 



fu+[ {Rf + ll}^0 
Jo Jgi ^ 



(3.363) 



The coefficients in the Fourier expansion of / will be denoted by and bk, i.e. / = ^^QafeCOsfc0 + 6fcSinfc0. 
We also introduce the map 0: [O,oo[x§i-^ [0,cx)[xF, (f,(^t, 6) :^ {R{t, e)cose + xo + xo, R{t,6)sm9 + yo). With 
respect to the parametrization 0, we redefine the zero-th order energy £(o)- 

£io)iu,f) = U {{l/R + uf-l/R]+ ( {^{R + fr + fi-R}. 
^ Jo Jgi 

The only difference from the definition (|1.29p is in subtracting 1/R instead of 1 in the first integral, and R 
instead of 1 in the second integral above. The higher order energies and are defined precisely as 
before. 

Lemma 3.14 The temporal energy £(u,f) is positive definite. 

Proof. Using the new conservation laws (|3.363p . by a small perturbation of p.338p and (|3.339p . as in 
Lemma |3.10[ we can bound the first and zero-th order momenta of /: 



k>2 k>2 



|5i| + |6i|<c(||VS|U2(a±) + C^(|afc| + |6fc|)), |ao|<c(||u|U2(o±) + ||VS|U2(o±)+^(|5fc| + |6fc| 

(3.364) 



With the above bounds we obtain an estimate analogous to (|3.350p : 

f(o)(",/)>Ci(||u||i.(f,±)+ / {/|-P/2})-C7i^/M||Vii||i.(f,±), 



where from we easily infer the claim of the lemma. □ 
The following lemma states the analogue of the stability estimate (ll.38p . where the energy quantities are 
expressed as functions of u and /. 

Lemma 3.15 The following inequality holds: 

1 



E^^,{uJ){t) + -j D^^,{uJ){T)dT<Ea.,,{uJ)[s), t>s>0. 

Proof. The proof of this lemma is completely analogous to the proof of ()1.38p . Formally, the only difference 
is the absence of the term C* (mg H-m^ + m^) in the statement of the inequality. The reason for this is 
Lemma I3.14[ which (again, due to our definition of u and /) enables us to conclude positivity of £{u,f) 
without adding any constant to it. □ 
From now on, we shall abuse the notation and denote u, i?, /, 9 and <j) respectively by u, R, /, 9 and (f>. The 
main ingredient in proving the exponential decay is to control the instant energy in terms of the dissipation, 
as stated in the following lemma: 
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Lemma 3.16 There exists a constant C>0 such that Ea ,^<CDa^i,. 



Proof. Note that the only term, which a-priori can not be bounded by a constant multiple of D, is precisely 
£(0) . Recall the expansion p.9p : uo(j>= —f — fgg + N{f) . For any k>2, let us multiply both sides of the first 
relation in (|1.9p by cosk9 and smk6 respectively, and integrate over S^. Observing that Jgi{— f — fee) cos k6 = 

i^^a.,weget 



{k^-iW\Js 



uo (hcoskO - 



[k^-imj^ 



N{f)cosk9 



(3.365) 



Note that 



uoAcoskO 



§1 



where we use the Sobolev inequality in the last estimate above. The last inequality, together with (|3.365p . 
immediately implies 

|afc|<C\/P + C / N{f)cosk9 , k>2 (3.366) 

and fully analogously, we arrive at 

\bk\<CVv + C [ N{f)smk9 

In order to estimate Oq, we first invoke the decomposition 

u°(t>\g\^-f- fee + q{f), 



, k>2. 



(3.367) 



(3.368) 

where q{f) stands for the nonlinear remainder with the leading quadratic order. Assuming without loss of 
generality that _R= 1 in the last relation in p.363p . we integrating (|3.368[) over S^, we find 

(3.369) 



u^- f+ qif) 



§1 



Multiplying the conservation law J^u + J^j^{f + /2} — by jij and (|3.369p by and subtracting the two 
equations, we obtain 



Note that C - 



_! i_ 

mi in 



-qn and hence \(- 



tL) \ <C(Af*)^/^, which, for M* small enough im- 



plies that 1 1^ — 1^ I > C/2 > 0. Hence, upon dividing p.370p by Ki :— jgrj ( j^jT ~ W\) ' ^® conclude 



— f- 

Ki\\n\ 



1 



2^1 1171 



I V |n 

i^iiri 



qif)- 



(3.371) 



From the mean value theorem, we deduce 

1 

m 



<C||VM||i2(o±). 



Thus, from p.37ip and the previous inequality, we deduce 

|ao|<q|V7.|U.(n±)+C||/||i. 
Summing (jg^T^ . ^MM, (P^^ and ^M7]), we obtain 



|ao| 



-Y,{W'^\ + \bk\)<C\\Wu\\L2(n±)+C\\f\\l,+CY,{\ / N{f)cosk0\+C\ / N{f)smke\} 

k=l k=2 •'^^ ''^^ 



(3.372) 



(3.373) 



<c||Vzi|U.(o±)+q|/|li.+q|iv(/)||L^<q|Vu|U.(o±)+c||/||^ 
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Using the smallness of from (|3.373p we obtain 



oo 

|ao| +^ (Iflfcl + \bk\) < C\ |Vw| U2(n±) < CVv. (3.374) 



fe=i 

Finally, from the previous inequality and the conservation law (|1.22p . we immediately deduce 



u 



<CVV. (3.375) 



The inequalities p.374p and (|3.375p imply the bound £(o) 1^ CT), which in turn implies the estimate claimed 
in the lemma. □ 

Using the same argument as in [21] to prove decay, we use the inequality p.38|) and Lemma [3. 161 to conclude 
that there exist constants fci , /c2 > such that 

Sa,.(M,/)(t)<fcie-^^*, i>0. 

This finishes the proof of Theorem 1 1.21 



4 The unstable case C < 
4.1 A-priori estimates 

In order to bound the higher-order space-time energy, we first introduce the necessary notations. We denote 
TVq :=971,(m^/J), m{ :=OT^(mJ,/J), for any j eN, whereby in the definition ((233)) . VL is substituted by VtK 
Analogously we define 0^^ and For any j e N wc finally set = Y.l=i and we define OJl-' analogously. 
Summing the e- independent and the e-dependent energies and Tl^ , we introduce dJll :=9Jl-' . In a 
fully analogous fashion, we introduce W , W and ^{ . Let us introduce 

In analogy to Lemma [3T2l we state the a-priori estimates for the unstable case C<0. 

Lemma 4.1 Let <^<0. Then there exist positive constants (3, Cb, o.'^d t'^ , such that for any <(,o, 
t<T'' such that if Mf3{0) <^/2 and 

M^{t)+ fN^{T)dT<^, sup ||nii.<CB(ei+ r^™) 

Jo Q<s<t Jo 



then 



MT+\t)+ fN^+\T)dT<^ and sup < ^^(6 +i T W"+i). (4.376) 

Jo 0<s<t Jo 

Remark. Note that the norm ||/e||//2i-i is already contained in the definitions of A4 and VJl and hence, 
at variance with the formulation of Lemma 13. 2[ the second claim of Lemma 14.11 is only concerned with the 
L^-norm of 

Remark. We claim that for any Q<k<l~l the following inequality holds: 

f^m.'mh+nn'it^uih+^'m^^^^ f%i\iu- (4-377) 

Jo ^ ^ V Jo Jo 

This claim follows by first applying the differential operator dtk to the regularized jump equation (|3.85|) . 
taking the squares and integrating over . By the standard trace inequality and the assumptions of smallness 
on and it is straightforward to get 

q=0 ' q=0 
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We now exploit the fact that = Aw'"+i to conclude that E^^o Il'^^l^"^^lli2(n™) < C'||Vu'"+i||^,fc(^„). 
Using the standard interpolation inequality and the definition of we obtain (|4.377p . 
Proof of Lemma 14.11 The second claim of (|4.376l) is proved in the same way as the inequality (|3.136p is 
deduced from the jump relation p.l33p . The proof of the first claim of the lemma is analogous to the proof 
of the first claim of (j3.115p in Lemma 13.21 The energy identities from Section [2] are nearly the same, except 
for the additional quadratic terms that appear in the definitions (|2.75p and (|2.76p of Q" and /S", respectively. 
We first analyze the temporal energies: the identity (|2.72p implies 

7W,(s)+ f Af,{s)ds = MM+ I I {0™'" + ^'"'"}+ [ [ {Q^'^ + S"'"}. (4.378) 

Here 0™'« EL=o = Eto^™'", 0"'" = Ei=oQr'" and 5"'" = EL=o^r'". whereby 

0'^":=0''{dt.u"'+^), Pf''":=P"(9t.u™+i) (4.379) 

and 

g^'":=g"(i?™,i?™+i,i?™), ^["'":=5"(i?;^,i?™+i,i?"). (4.380) 

Estimating RHS of (|4.378p is analogous to the estimates from Subsection ^231 Namley, O™-" = O™, P"'" = 
P™ and 



Qm.u^Qm^y^f f {(Po+Pi)/™+l(Po+Pl)/,T++i'+ePl/e^,t'Pl/,^+ii} 

= + ^ f f {((Po+Pi)/,T++.')' + e(Pi/,'?,tiO'}. 



We estimate the second term on RHS of the first equation in (j4.38ip in the following manner: 

IE r / {(iPo+Pi)/r'(iPo+Pi)/,T.+.^+^Pi/,'^,rPi/,^^i.}|< fmt'wh+^wfrn'^.wh} 

Jo ^ ^ Jo Jo V Jo 



(4.382) 

where, in the second inequality we used the estimate (|4.377p and the second inequality in (|4.376p to bound 
loiWf'^^^Wh+'^Wfei^^Wh} by the last two terms on the right-most side of (14382]) . The second term on 
RHS of the second equation in (I4.38ip is again estimated using the estimate (|4.377p : 



e/7 {(iPi/"^i')'+^(iPi/«5,tij'} <r?r^r 

u^nJO JS^ Jo 



.1 , c r 



k=0 



V Jo 



m+1 1 ^ / I l„,m+l 1 12 



(4.383) 



Using the identity (|4.378p and the estimates analogous to p.l44p - (|3.150p . (|3.155p . (|3.162p . (I3.164P (I3.170p . 
(|XT7T|1 - ([?T77|1 . (jXWIl . (jXTgO)) and (jXTM]) - ^J7^, together with the additional estimates 
and (I4.383p . we arrive at the following inequality: 

/ AC™+'<Xe(0) + -^(Al'"+Al'"+i) + (-^ + -^)A^"+i 

Jo ( ( ^^^^ J J 

+C7(A+ sup \\r\\H^+Vm^){ [ AC"+ / AC"+')+^(6+i / W"+') + ^ / (4.384) 
0<s<t _ Jo Jo ^ Jo ^ Jo 

+ {C\+^ + ^ + C^M^) f W+^+C f N"' sup OT"+i(s) + ^ / 
A eA Jo Jq o<s<t A Jo 



,m+l||2 



We now turn our attention to the estimates for the space-time energies and W . Note that the esti- 
mates p.232p . p.238p . (|3.24ip carry over to the unstable case, where the energy terms D™, 2?™+^ 
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are replaced by 01™, fR^+i and 7V'"+\ respectively. The identity (|3.25ip reads 



"T" 



+1 /^{|/,"iVP-|P2+/;i++VP} ^= / 



(4.385) 



Js^ 



where is defined by (|2.75p . For the sake of simplicity, let us denote 

0"(/e™ +itc,/^?++\tc,i?'Mff"|'^+') and note that 

Am.ti _/nm , f f m fm+1 m fm+1 i , /" f W fm+1 p rm+1 

^2p+l,C ~ V2p+l,C ^ / / ^l./fl2p + 2tcJi 1 jg2p + 2tc+l / / 1 J g2p + 4tc 1 ./ g2p + 4(c+ 1 • 

Jo JSi Jo Jgi 

The estimates analogous to (|32S11), (ElSSl), (E^Si), (02231) - (EMI) imply the bound on g^;,_^i_^: 



'2p+l,c 



<c{t+Vm^+x+ sup linuOf / K+ I ^"+') + (? + tt) 

0<s<t ^Jo Jo / A Ae 0<s<t 

(4.386) 



Furthermore, 

i 



5 /.m+1 p i-rn+1 i c f [ W f™+l p 

IJ g2p + 2fc'^ IJ g2p + 2fc+l ~r t ; Jri7g2p + 4tcJl lJ92p + 4jc+l 

/o Jgi Jo Jsi 



<c 



/o Jsi 



< A j^* I iPi/™ V I ii2 + f ^* I iPi I ii2 < + ^ 

(4.387) 

where we used the estimate (|4.377p in the third inequality above. Observe that the identity p.278p reads 



a,2p9;.,,^'"+ia,2pa*K+i]±= / / {|/,^+\,e+4p-|P2+/,^+i+ip} 

^0 t ^ (4.388) 

I jm+l |2 |m frn+l |2\ / / qu / fm fm+1 r>m \ m\2p\ 

|/g2p+3(c+i| F2+ Jg2p+2tc+i I f— 11'-' Ue^pf'Je^pt"^-^ 'If I ^' 
Jsi Jo Jsi 



where S"" is defined by (^TSD . For the sake of simplicity, let us denote S*^'" := S''{rjip^^J^+l,R"\ |g™pP+i) 
and note that 



'5'2^;c— '5'2p.c+ / / |Pi/^i\tc+iP + e / / |Pi/( 
Jo Jsi Jo Jsi 

The estimates analogous to p.279p - p. 2831) imply the bound on 



m+l 

■p 



\s^pJ<c{Vti^+\+ sup ||/™|U2)( / {m^+m^)+ [ {m^+^+m^+')) + ^ [ m^+' + -^m. 

Q<s<t ^Jo Jo ^ Jo ^AE 

(4.389) 

Furthermore, analogously to (|4.387p . using the estimate (|4.377p we obtain 



Jsi Jo Jsi 



(4.390) 

From the argument in the paragraph preceding (|4.385p and the estimates (|4.386p - (I4.390p . we obtain the 
energy inequality analogous to (|3.285p : 



^ Jo Jo ' , , , 

l-t v-i i-t i-t (4.391) 

<97lep(0) + !^' / W'+^ + CY^ / W„'+^+KP + — sup M"'+^ + C* / AA'"+i 

^0 ~:^iJO e 0<s<t Jo 
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where is a small parameter and is given by: 

icp^c{t+Vm^+x+ sup ||/™||i.) / {K+K) + {^ + Ti) "^T- 

o<s<t Jq '^e o<s<t 

Using the inequalities (|4.384p and (|4.39ip in the same way as it is done in Subsection l3.2.6[ we conclude the 
proof of Lemma UTJ □ 

4.2 Local existence 

The goal of this subsection is to prove the analogue of Theorem 13.41 stating the local existence, uniqueness 
and continuity of solutions to the regularized Stefan problem over a time interval that does not depend on e. 
In addition to that, we state the fundamental energy estimate that will be used to establish the instability 
in Subsection 14.41 



Theorem 4.2 Let (<0. For any sufficiently small L>0 there exists t>0 depending on L and constants 
Q < <1 such that if 

then there exists a unique solution (u*^,/*^) to the regularized Stefan problem defined on the time interval 
[0,t[, satisfying the bound 



sup Ml^{u\rf){t)+ / Nl,{u\f%s)ds<L 

0<t<T Jo 

and ^{u"^ , f'^)(-) is continuous on [0,r[. 

Proof. The proof of the local existence on an e-dependent time interval [0,t^[ together with the proof of 
continuity is fully analogous to the proofs in Subsections 13.3.11 13.3.21 and 13.3.31 In order to prove that the 
solution exists on a short time interval independent of e, we proceed analogously to the Subsection 13.3.51 
We first pass to the limit as m ^ oo in the inequality l|4.39ip . All the terms involving C, C in the estimates 
vanish in the limit, because they are used to bound the cross-terms. Thus, in turn we obtain an estimate 
with constants that do not contain e in the denominator: 

mp+ f mp+{i-C2t){m;+ f %)<6p+iy' f mp+i+cj^ f m,+iCi+c* f f \\iu,f)\\l., 

JQ Jo Jo q=0-^^ "^0 '7^ Jo 

(4.392) 

where ||(u,/)||2,:=||«||2,^^^ + ||/||2^ ^nd ICi=C{t + Vm + X + snpo<s<t\\f\\L-) lo + Note that the 
key novelty with respect to the energy estimates in the stable case, is the presence of the (potentially) 
large multiple of the L^-norm of (u,/) on RHS of (|4.392p . Multiplying (j4.392p by and summing over 
p = 0, — 1, just like in p.355p . we obtain 

sup m'^+J^dp %, + {!- C2t){ sup m'-''+ / ^'^''^) <9Jl^(0) + e^z^'?/Ci+^i/«C* / Af 
o<s<t Jo o<s<t Jo Jo (4.393) 



77A 



1 pt 



IIK/)II 



L2, 



We now pass to the limit as m— >oo in (I4.384p . Note that the term '^^^ (^1 +t /q^™'*'^) used to estimate 
/o Il/™^^lli2 can, in the limit, be replaced by the bound /„ Il/Ili2 < C/J IKm,/)!!^^. We obtain 



Me+ I K<MM+c{>^^ sup \\f\\H^ + Vmi) f j\fc + c{\+VM) f m+c f N sup 9JI(s) + y 

Jo 0<s<t Jo Jo Jo 0<s<t X 

(4.394) 

Combining the estimates (j4.393p and (|4.394p in the same way as the estimates (|3.355p and (j3.356p are used 
in Subsection 13. 3. 5[ we conclude that there exists a t > such that > r for any e. □ 
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4.3 Variational framework and linear growth 

Let (/)o : §^ — + 5i(xo,0), (j)o{0) — {xo + cos9,sm9) be the parametrization of Si{xo,0). In order to find the grow- 
ing mode Xq of the Stefan problem with surface tension, we turn our attention to the associated eigenvalue 
problem 

Xv-Av = onr2\S'i(xo,0), (4.395) 

vo<t>o = -f-fee on§i, (4.396) 

[vn]to^Q = -\f onS\ (4.397) 

t;„ = ondVL; [v]t^Q onS'i(xo,0). (4.398) 

For any A:eNU{0} set 

H^{VL^) = ^u:n^m.;ul^i(n±)^H^{c\{^^))], (4.399) 
where fl"*" := f7\cl(i3i(a:o,0)) and fi^ := Bi{xq,Q). We define the linear operator C 

C{vJ):={Av-[d,Mt), (4.400) 

where the domain of C is given by 

D{C) {(i;,/) e H^{n^) x i/^/2(§i); = on 9^!; [v]t^Q on ^i(xo,0); i;o0o = -/- /ee on S^}. (4.401) 

The following lemma states a number of important properties of the operator C and the proof can be found 
in [IH]: 

Lemma 4.3 There are countably many eigenvalues Ao,Ai,A2..., they are all of finite multiplicity and the 
associated eigenvectors are smooth. There is exactly one simple positive eigenvalue Xq. Furthermore, Ai = 
and its eigenspace has dimension three. All other eigenvalues are negative. 

Let us define a bilinear operator (•,•) on D{C). Let {v,f), {w,g)^D{C). Set 

{{v,f),{w,g)):^ / vw- / f{I + d0g)g^ / vw+ / {fege-fg) 



Lemma 4.4 The operator C is a symmetric operator on D{C), with respect to (•,•). Furthermore, the 
eigenvectors associated to any two different eigenvalues are orthogonal with respect to (•,•). 

Proof. Where no confusion can arise, we abuse the notation and write v = —f— fgg on instead of w o (/jg = 
—f — fee- Analogous remark holds for [vn\^ = —Xf on S^. A straightforward computation yields: 

{C{v,f),{w,g))^{{Av,-[dnv]t),{w,g))^ ( Avw+ ( [d,,v]t{I + dee)g 



Awv+ / [i9„v]lw— / / [w„]lo0ow= / Awv— j [dnw]^v 

n Jgi Js^ Jn 



Awv- / ~[dnw]tiI + dee)f = {ivJ),iAw-[dnw]t))^{ivJ),Ciw,g)). 
n Jsi 

Note that in the second equality, we used the Green's identity to rewrite J^Avw and we also used the 
boundary condition w~—{I + d0g)g to rewrite J^i[dnv]'^{I + dgg)g . Note also that the third equality implies 

{C{vJ),iw,g))=- f VvVw. (4.402) 
Jn 

If A and fi are two different eigenvalues (w.l.o.g. ^^0), with associated eigenvectors (w,/) and {w,g) 
respectively, then 

{C{vJ),{w,g))^X{iv,f),iw,g))^-{{vJ),Ciw,g)) 
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and hence, due to the symmetry of (•,•), we have {C{v , f) , {w , g)) =0. □ 
To every eigenvalue A, with multiphcity m{i), we associate the eigenvectors e^^i, . . . ,6^ where we denote 
Si,k = ivi,k: ^i.k) for 1 < ^ < m(i) . The null-space S of the operator C is spanned by three eigenvectors ei i = 
(1,-1), ei_2 = (0,cos6') and 61,3 = (O,sin0) associated to the eigenvalue Ai=0. Let us define the functional 
space I?! by: 

:= { (f, /) e (n^) X (§1 ) ; vn = on 9^!; [v]t = on Si{xo,0); vocj)o = -f - fee on } (4.403) 
By 5^ C -Di we denote the orthogonal complement of S in Di with respect to (•,•). In other words, 

S^ = {ivJ)eDi\ jj + j^J^^^ Fi/ = 0}. 



Lemma 4.5 The following variational characterization of the positive eigenvalue Aq holds: 

— mm 



Proof. Without loss of generality, we assume xq — 0. For (w,/) G let us set 

(4.404) 

Let t = inf(„j)g5iX(w,/). 

Boundedness from below: let (w, /) G 5^ . Since /q + /gi / = and Jgi v — —J^if, we deduce that JqV — JgiV. 
By the mean value theorem | fny /fj'^ ~ Ifrj/si ''^l — ^1 ^^"^ ^-'^^^ ICI | ^ C"! I V'f'l lL2(n±), where 
Il'lli2(n±) is defined by (I4.399p . Since \(^\ >0, we conclude 

(ivj),{v,f))^ I \Pv\'+ I (/,2-|P/n + C( / ^')'>C( / t')'>-yl|Vz;||i.(o±), 
JQ. J& Jn Jn s 

and hence l> —00. 

Negativity: l<0. Assume for the moment xq = 0, i.e. S'i(a;o;0) Let 

{d 0<r<l-e and l + e<r<i?*, 
d+i_i(r-(l-e)), l-e<r<l, 
t)d,,(2-r), l<r<l + e. 

Note that Vd^e{l) = l and we define fd.e — —^- Hence IPi/c;,£ = for any d,e. For a fixed d note that J^Vd,e = 
RldTT + 0{e) and thus the condition jQVd,e + J§ifd,e = reads R^dir —2t: — 0{e) and hence (i=2/i?^-|-0(e). 
On the other hand, it is easy to see that /^^ |Pvd,eP = 0(e). Note however 

I{vd.,e,fd.,e)^<:i f ^d.^Y+l |Pt'rf..|'=Ci?*dV2+O(e)<0 

Jn Jn 

for small e. To conclude that t<0, we slightly perturb the function Vd^c to a smooth function v^^g so that 
the |Vw„egP is well defined. (If xq t^O, we define the function u^^g inside S'i+e(xo,0) by simply translating 
the function Vncg by xq. On Si,;, (0) \ 5'i+e(a:o,0) set v^°g = d. It is easy to see that X{v^°g, — 1) <0.) 
Existence of the minimizer. Let (w„, /„) be an infimizing sequence for the above variational problem. Without 
loss of generality assume |Vw„p = 1 (otherwise rescale). As in the proof of boundedness, 

I / v„\^ <^\\Vvn\\L^n±) (4.405) 
Jn ^ 

and hence from here and the Poincare inequahty | |wri||L2(o±) < |Vw„| |i2(Q±-) = C. Due to (|4.405|) . Pi/„ = 
and the condition v\E>^ = — f — fee , we deduce that there exists a positive constant C such that ||/n||//2 <C 
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for any n. By Banach-Alaoglu, we deduce that there exist weak hmits v* and /* of {vn} and {/„} in H^{n) 
and H'^{S^) respectively. Hence {{vn,fn),{vn,fn))^{{v*,f*),{v*,f*)) as n— >oo. Furthermore, by weak 
lower semi-continuity of || • | jj^i , we deduce 1 =liminf„^oo | IVz^nl |L2(fj±-) > | [Vi;*! |i2(-Q±). Since 1(vn,fn) < 
for n large enough, we deduce 

;, = lmimf2(u„,/„) >limmf — ^ =I{v ,/ ) 

and hence {v*,f*) is a minimizer. Note that ||Vw*||£,2(s-j±') = 1. Otherwise assume ||Vw*||^2(s^±-) < 1. Then 
((w*,/*),(w*,/*))=i||Vi;*||22(j.2±pi^liminf„^oo(K,/n),K, /«)> = ((«*, /*),(«*,/*)), which is a contra- 
diction (note that we used the fact t<0). 

Euler- Lagrange equation. Setting i = — , the Eulcr-Lagrange equation for the variational problem of min- 
imizing I over allows to conclude that the minimizer (w*,/*) is the weak solution of the eigenvalue 
problem £{v,f) = Xo{v,f), {v,f)£Di. In other words, for any test function (^s e C°°(ri^): 

I {Xov*ip-\7v* ■\7ip}-Xo I f^p = 0. (4.406) 

By standard elliptic regularity, it is easy to see that v* is smooth in the interior of il+ and fl^ . We want to 
prove that v* is up to the boundary i.e. v *eif2(f7±) (recah (l4399ll ). 

Claim. Let 71 G H^{Q^) and 72 £ H^{S^) be given. Let v £ H^{n^) be the weak solution of the problem 

Au = 7i in Q, 

[dnv]t^j2 on §\ 
dnV = on 

Then v£H^{VL^). 

Proof of the Claim. Written in polar coordinates, this means that for any test function (^GC°°(ri^), 
the following holds: 

( [ {-Vriprv-^^drde}- [ 72(^d6l= / / -fitpdrde. (4.407) 

Let Vgf (r, 6>) = v(r.0+h)-v(r,e) r^O. Let ip£C°°{Q) be supported on the complement of a small ball 
Bs{Q). From (|4.407p . we immediately obtain 



/ / {-S/'^Vr'frr-^^^^^drde}^ [ ^'Sj2<fde^f [ V^7i(/?rdrci6'. 

Js Jgi ^ JSi ^(5 Jgi 

Let s>0 be a fixed small real number. Let (p£H^{n^) be a test function such that ip{r,6) — —S/gV for 
(r,6l)eOi :=[l-s,l + s] x§i and supp(v3) G x Let us also denote O2 -.^ {[6,1- s]U[l + s,R^]} xS^. 

Then 

^ -^s^ (4.408) 

-/ V'^-fiV^iprdrd0+ {V'^vriprr+-^^^drde} 
J01UO2 J02 ^ 

Note that 

/ V^72V^i;d0 = / V-''V^72fd0 <||72||//2|kl|L2(si)<||72||H2||z;||^i(n±), 
where we used the trace inequality in the last estimate. Furthermore, 

/ V^,'7lV^V'■'^'^^^^'<l|V^'7l||L2(OlUO2)l|V00||L2(OlUO2)<<^ll7l||i^H^^±)(ll"llffU^^±) + ll'/'llffMO2))• 
JO1UO2 
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Using the interior regularity, it easy to bound the remaining term on RHS of (|4.408p : 



{V^VrLPrr+-^ drdO) < C||t;||//2(o^)||v3||^i(o^) < C||t;||j^i(o±)||<^||^i(n±) 

Hence, we may pass to the limit as h-^oo in (|4.408p to conclude Vr,vg£ Hg(Q^). Since Av^ = {drr + ^ + 
we deduce that 

drrV^ =7ilci(n+) — (— H — Y')''^^ distributional sense. 

Since the RHS of the above expression actually belongs to L^(cl(r2+)), we deduce € L^(cl(r2+)) in distri- 
butional sense and thus is up to the boundary We argue similarly for v~ to conclude the proof 
of the Claim. 

We now apply the above Claim to the function v* £ H^{n^). Upon setting 71 = Xqv* and 72 = — Aq/*, we 
conclude v* £ H'^{n^). Hence Aq is a positive eigenvalue for the linearized Stefan problem, with an associated 
eigenvector (u*,/*). Note that w*|si € iJ^/^(S^) and the Dirichlet boundary condition v*\si = — f* — fgg 
implies /* G i/^/^^gi)^ gy ^ classical bootstrap argument, we conclude that (u*,/*) € C°°{n^) x C°°(§i). □ 



Lemma 4.6 Let eo.i = (wo,i,zo,i) be an eigenvector associated to the growing mode Aq, for the Stefan problem 
linearized around 5*1(0,0) =§^. Then vo,i is spherically symmetric. 

Proof. Let vo^i{r,9) — J2'kLo-Pk{r,d) be the Fourier decomposition of vq^i. Due to the linear independence 
of I cosfc6',sinfc6'}^, we deduce that each Pk solves the eigenvalue problem (|1.15p - (|1.18p with A = Ao >0. If 
fc>l, it is easy to see th at ((Ffe,Pfc(l,0)), (Pfc,Pfe(l,0))) =^F| + J^ddePkM)^ - (PkM)^ >0- On the 
other hand, using (I4.402p and the previous observation 

- f \^Pk\^^{C{Pk,Pk{l,0)),{Pk,Pk{im^^o{{Pk,Pk{l,0)),{Pk,Pk{l,0)))>0. 
Jn 

This is only possible if Pfc = for all fc> 1. From here we immediately deduce that (wo,i:^o,i) — (^o(^)ic) for 
some function Pq and a constant c= — Po(l)- 1^ 
Now, we shall turn 5^ into a normed space, by introducing a new norm || ■ ||/ for any fixed />0. The 
new scalar product (•,•)/ is defined by: 

{{vJ),iw,g)), = {{vJ),iw,g))+{^ + I) f VvVw, ivJ),iw,g)eS^. (4.409) 

^0 Jn 

From Lemma 14.51 we immedaitely see that ((«,/), (f ,/))/> /| |Vu| ||2(q±') > 0. Furthermore, if 
{{v,f),{v,f))j = 0, then ||Vw||i2(f2±) =0 and hence v= const, which imphes {v,f) — /ci(l, — 1) + fc2(O,cos0) + 
A:3(O,sin0) = fcici^ + ^2^1,2 + ^361,3, for some reals ki,k2,k3. Since {v,f )£ S^, it follows («,/) = (0,0). Thus, 
is an inner product on and it defines the norm || • ||/ := \/ (•,•)/• 

Lemma 4.7 There exists a constant C >0 such that for any {v,f)£ 5^, I kllL2(n) + 1 1/| 11,2 < C||(f ,/)| |/. 

Proof. As in the proof of Lemma |4.4[ we abuse notation and write vo(j)Q — v. Let {v,f)£S^. As in 
the proof of boundedness in the previous lemma, we conclude |<^|Uj^f | < C||Vw||2^2(q±). From this and 
the Poincarc inequality, we deduce ||w||i2(Q') < C||Vi;||/,2(f2±-). Furthermore, from v\S^ — —f — fgg, we ob- 
tain /gi/e'-|P/P<C||f|U2(si)+C|/gi/|2. Since Pi/ = and /| = < C||i;||l2(si), from the trace 

inequality we deduce 

ll/lli^=^/«-|p/l'+^(^/)'<c^llHli^(s.)<ciHil,i(o±)<c||vt;|ii2(o±). 

Since ||(w,/)||/ > /||Vu||^2(q±-), we conclude the claim of the lemma. □ 
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Lemma 4.8 The set B — LI^QLIk<m{i) {^i.k} forms an orthonormal basis of with respect to (•,•)/• 

Proof. Without loss of generality, we assume xq = 0. Let A and A* be two different eigenvalues with associated 
eigenvectors e = (Wjz) and e* = Lemma H^ imphes (e,e*)=0. From (|4.402p . we deduce J^Vv-Vv* =0 

and hence {e,e*)j = 0. Let us assume span[S] and let denote the orthogonal complement of span[S] 
with respect to (■,■)/• Let us set 

J:— inf 2(v,f), 
where I is defined by (|4.404p . We distinguish two cases: J < and J > 0. 

Case 1: J <0. Note that J>i, where t is defined in the line after (|4.404p . The same proof as in Lemma |4T5| 
yields the existence of a smooth minimizer {v,f)^B^ such that I{v,f) = J. From the symmetry of C, we 
deduce 

{^vJ),{v^^kJr,k))j = {{v,f),C{v^^k,fi,k))j = {{vJ)At{Vt,kJt,k))j^O, 

Since {vJ)eB^. Rence £{vJ)gB-^ and if we set J = —j, the Euler-Lagrange variational equation implies 
that A is an eigenvalue of C with an associated eigenvector {v,f)G B'^ , which is a contradiction. 
Case 2: J> 0. In this case, we analyze the problem of minimizing 

1 L^^v? 



I{vJ) {{vj),{vj)) 

over the set B^. Since J>0, it is clear that jiirj) ^'^^ ("^^f) S'S^. Set K — mi(^^ j^^g± jj^-Yj ^^"^ 
{vn,fn) be an infimizing sequence. After rescaling assume {{vmfn)Avn,fn)) = '^- Just like in the proof of 
Lemma [4.51 we conclude that there exists a C>0 independent of n such that for any neN | I'^nl |ffi(n±) + 
||/ti||h2 < C. We may pass to a weak limit (w,/) by Banach-Alaoglu. From the compactness of the embedding 
^ and H^^H^, we conclude {{vn,fn),{vn,fn))^{{v,f),{v,f)') as n^oo. Furthermore, by weak 
lower semi-continuity, liminf„_>oo ||Vt'„||i2(Q±-) > ||Vu||i2(s^±-). Thus 

if^liminf — -, — > ■ 



and thus (u,/) is the minimizer. If we set K = —\^ the Euler-Lagrange variational equation implies that A 
is an eigenvalue of C with an associated eigenvector (w,/) (iB^, which is a contradiction. This finishes the 
proof of the lemma. □ 
For a given yG-Di, let y ~ys + , where ys is the projection onto S and t/^ projection onto S^, with 
respect to (•,•). We define 

\\y\\'--=\\ys\\l^n).mB^) + \\y^\\l (4-4io) 

Due to Lemma W77\ we immediately see that 

||2/||L=(a)xL^(si)<C||y||. (4.411) 



Lemma 4.9 Let (u^Oj/o) = X]i^oSi<fc<m(i) "^^-fc^^fc- Then there exists a constant Cc, such that 

\\e^'{woJo)\\<Cce^°'\\{woJo)\\- 

Proof. Note that 

||e^*(wo,/o)|P = II «i,fcei,fc+ «»*=e^fcll^ 

l<fe<3 i>0,j:#l l<fc<Tn(i) 

= 11 Y «i,fcei.fclli2(o)xL^(si)+e'^"*l«o,i|'+^e2^'* 

l<k<3 i>2 l<k<m{i) 



i,k 

and the claim follows. □ 
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4.4 Proof of Theorem 11.31 

The abstract formulation. We shall use the bootstrap method developed in [in] to prove the instability 

of the stationary surfaces in the case C < 0. The following lemma is a simple modification of Lemma 1 in [15], 
which suffices to establish the nonlinear instability of steady spheres for the Stefan problem. 

Lemma 4.10 Assume that X is a Banach space with the associated norm \\-\\x o,nd C:D<zX-^X is a 
linear operator. Assume a norm || • || with the property || ■ ||x < C*!! ■ || o^nd such that e^'-' generates a strongly 
continuous semigroup satisfying 

\\e'^\\x^x<Cce^°' (4.412) 

for some Cc and Aq > (l|-||x^Ji: stands for the operator norm). Assume a nonlinear operator U -.D^X, 
a norm \\\-\\\ on D and a constant C'u , such that 

\my)\\<Cu\\\y\\\' (4.413) 

for all y^D and \ \\y\\ \ <oo. Assume that there exists a small constant a such that for any solution y[t) to 
the equation y' — Ly + U{y), with \\\y{t)\\\<a , there exist constants Ci and C2 and a small constant ^, such 
that the following energy estimate holds: 

\Mt)\\?<Ci\\\y^\\\^ + C2 f\\y{s)\\'ds. (4.414) 

Consider a family of initial data y^ [Q) — Syo with\\yQ\\ — \ anc? |||j/o| II < 00 and let 9q be a sufficiently small 
(fixed) number. Then there exists some constant C>0 such that if 

0<t<r:^ilog^, 



\y{t)-6e'''yo\\<c{\\\yo\\\' + l)s' 



2„2Aot 



In particular, if there exists a constant Cp such that \\Se^*yo\\>CpSe''^°'^ , then there exists en escape time 
<T such that 

\\yiT%>Ko>0, (4.415) 
where Kq depends explicitly on C'c, Cjj , C2, Cp, yo, Xq and is independent of S. 

Proof of nonlinear instability. Let M**<^, where L is given by Theorem 14.21 and let {u,f) be the 
unique solution of the regularized Stefan problem with the initial data satisfying ^.(uoi/o) < A/**. We 
multiply (|4.393P by P and sum it with (|4.394p . In analogy to (|3.357p and (|3.358p . we choose A, M**, v small 
and use the smallness of ||/||^; to arrive at 

sup M^,^(s) + i / Xii,^{s)ds^{\-Ct) sup M/3,^(s) + i / /V0,^(s)ds 

0<s<T ^ Jo 0<s<T ^ Jo 

<M^,40) + M/3,,(0) + e/C2 + C / \\{u,f)\\\.As, 

Jo 

where the constant /C2 is defined analogously to J2 as in the line after (j3.357p . Passing to the limit as 0, 
we obtain a solution (w,/) to the original Stefan problem. Note that M/3_jy(0)— +0 as e^O, by the choice of 
the initial condition. Also, e/C2— >0 as e— >0 because /C2 is bounded. As a result, we obtain the estimate 

sup M^,,(u,/)(s) + i / NpA^J){s)ds<M^AuoJo) + C f \\{u,f)\\l,ds (4.416) 



0<s<r 



as stated in Theorem ll.3l To apply Lemma [4 .101 and finish the proof of Theorem ll.31 we have to formulate the 
Stefan problem on a fixed domain. To this end we apply the change of variables described in Subsection l3.3.11 
we set x — TT{t,x)x, where tt is defined by dropping the super-indices in the definition ([3.289p . Without loss 
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of generality, we furthermore assume xq — 0. Hence, setting v:=uot: leads to an equation for the function 

Vt — Av~aijU^i^j +biU^i, (4.417) 
where and bi are defined by dropping the index m in (|3.29ip . From (|3.292p . we deduce 



ftR' 



on §\ 



(4.418) 



The Dirichlet boundary condition simply reads: 



-.f-fee + N{f) onE\ 



(4.419) 



In order to set the problem (|4.417p - (|4.419p in the framework of Lemma I4.10[ we have to consider the 
new unknown u;(a;) := u(x) — -/V(/( j|j)). This way we insure that the nonlinear boundary condition (I4.419P 
transforms into a linear condition for w. We obtain 



wt- Aw^TZ{w,f), inilgi, 
w^-f-fee onS\ 



_ = 



Here 



Wn = on do.. 

n{wj) = a,,w,.,.+M,.-(9t-A)iV(/( — )) 

= a,, (z« + 7V(/(-)))^.^, +h.,{w + N{f{-)))^^ - (9t - A)iV(/(-)). 



(4.420) 
(4.421) 

(4.422) 

(4.423) 



To apply Lemmaimniset X := ^^(f^) x L2(§1) and D = Di as defined in (I4.403p . The norm \\-\\ is defined 
by (|4.410p and we also define 

Observe that if {w,f)GX is a solution to the Stefan problem (I4.420p - (|4.423p . then the norm |||'||| is 
equivavlent to the norms defined by i.. This follows easily from the boundedness of the Jacobian of the 
coordinate transformation x — 7T{t,x)x and the fact that the time derivatives of w and / can be recursively 

1 2 , 

expressed in terms of purely spatial derivatives via the relations wt = Aw + F and ft^ — [w„] _ . From this 
remark, the bound (|4.41ip and the inequality (|4.416p . we immediately conclude that there exist constants 
Ci and C2, such that the inequality (|4.414p is satisfied. Furthermore, we define 



[/(«;,/) = (7^(«;,/),-(^-l)K]^). 



(4.424) 



Note that the time derivative of / that appears in the definition of TZ{w,f) can be alternatively expressed 
as a purely spatial differential operator via (I4.422p . With this notation, the problem ()4.420p - (I4.423P takes 
the form dt{w,f) — C{w,f) + U{w,f). It is straightforward to verify ||[/(u;,/)|| <C|||(w,/)||p and thus, the 
assumption (|4.413p from Lemma 14.101 is satisfied whereas Lemma 14.91 guarantees (|4.412p . We may apply 
Lemma 14.101 to conclude that p.4ip holds and this finishes the proof of Theorem 11.31 □ 
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